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.

A Cartesian tensor of order $n$ is a multilinear geometric object whose components under an orthogonal change of axes obey the tensor component transformation law,An isotropic tensor has the same components under every proper rotation. For the Kronecker delta,Thus ${\delta }_{ij}$ is isotropic, even under improper orthogonal transformations.

To prove multiplicativity of the determinant from the given formula, first consider $C_{pqr}={ϵ}_{ijk}{A}_{ip}{A}_{jq}{A}_{kr}$. Exchanging any two of $p,q,r$ changes its sign, so it is proportional to the Levi-Civita symbol. Its $123$ component is the determinant, as follows by expanding the given double-epsilon expression. ThereforeSubstitute $(AB{)}_{il}=A_{ip}{B}_{pl}$ into that same determinant formula:HenceNo invertibility assumption was used, so this identity includes singular matrices.

Now differentiate the double-epsilon formula for a differentiable $A(t)$. Its three differentiated terms are equal after renaming the paired dummy indices: the simultaneous swaps in both epsilon symbols cancel their signs. Thuswhere $C_{il}$ is the cofactor of entry $A_{il}$. For an invertible matrix, the adjugate matrix identity gives $C_{il}=(\det A)(A^{-1}{)}_{li}$. Consequently the Jacobi determinant derivative formula isThe order of the factors inside the trace follows from the reversed cofactor indices; no commutativity of $A$ and $\dot{A}$ is assumed. The cofactor version is also Jacobi's formula and remains valid when $A$ is singular, although the inverse version then cannot be used.

For integrable $f,g$, define their convolution by $(f\ast g)(x)={\int }_{\mathbb{R}}f(y)g(x-y)dy$. Absolute integrability of the double integral permits Fubini's theorem. With $s=x-y$,The convolution theorem, together with uniqueness of the Fourier transform, therefore identifies the inverse transform of that product with $f\ast g$. With the stated convention the inverse carries $1/(2\pi )$, so there is no extra factor in the convolution itself.

Now take $f=e^{-|x|}$. For $x\ge 0$, split the convolution at $y=0$ and $y=x$:The convolution of two even functions is even, so the self-convolution of a bilateral exponential givesIn particular the centre value is one; both pieces extend continuously there.

Splitting the Fourier integral at zero givesHence $\tilde{h}(k)=4/(1+k^2{)}^2$. The Fourier inversion theorem yieldsThe integral is absolutely convergent. At $x=0$ it gives ${\int }_{\mathbb{R}}(1+k^2{)}^{-2}dk=\pi /2$, providing a normalization check.

The $k_0^{2}uv$ terms cancel in $u{L}_{H}v-v{L}_{H}u$. For a ball $B_R$, the product rule and divergence theorem give Green's second identity in the formThus the whole-space integral vanishes when this boundary term tends to zero. For example, uniform bounds $u,v=O(R^{-1-ϵ})$ and $\nabla u,\nabla v=O(R^{-1-ϵ})$, with $ϵ>0$, make it $O(R^{-2ϵ})$.

Pointwise vanishing alone, as stated in the source, is insufficient. Smooth fields agreeing with $u=\sin (k_{0}r)/r$ and $v=\cos (k_{0}r)/r$ outside a ball both tend to zero, but their distant boundary integral is $-4\pi k_0$. One may smooth them inside the ball without changing that flux. The requested zero identity therefore needs sufficient decay or boundary/radiation conditions that actually eliminate the boundary term.

For the generalized plane wave eigenfunctions, their orthogonality is distributional:Here $\int e^{iqx}dx=2\pi \delta (q)$ is the Fourier representation of the Dirac delta function. These eigenfunctions are not ordinary square-integrable functions or decaying test fields.

Expand $\Phi =\int A_k{\psi }_k{d}^{3}k$ and define $V_k=\int {\psi }_k^{\ast }(r^{\prime })V(r^{\prime })d^3{r}^{\prime }$. Since $L_H{\psi }_k=(k_0^2-k^2){\psi }_k$, projection of the source equation gives $(k_0^2-k^2)A_k=V_k$. Formally dividing and substituting back givesAt the shell $k=k_0$, a pole prescription is required; homogeneous solutions can otherwise be added. The formula is an inverse of the Helmholtz equation with that convention still to be selected.

Let $H$ denote the kernel in this source formula, so $L_{H}H=+{\delta }^{(3)}$. Put $l=|r-r^{\prime }|$ and align the Fourier polar axis with $r-r^{\prime }$. Angular integration givesConsequently the radial Fourier representation of the Helmholtz Green function givesThe last Green function written in the PDF has the opposite denominator. It is obtained by setting $G=-H$, the common convention $L_{H}G=-{\delta }^{(3)}$:Thus the two displayed source and Green formulas are consistent only when this minus sign is retained; $G$ cannot also be the positive-delta inverse kernel.

Finally the integral is not an ordinary improper integral across $k=k_0$. A Cauchy principal value yields $G_{\mathrm{PV}}(l)=\cos (k_{0}l)/(4\pi l)$. The outgoing prescription $k^2-k_0^2-i0$ yields the outgoing Green function for the three-dimensional Helmholtz equation, $G_{+}(l)=e^{i{k}_{0}l}/(4\pi l)$; the opposite sign gives the incoming wave. The Sokhotski–Plemelj formula accounts for their difference. The source does not specify which prescription is intended, so the boxed radial expression must be read with one chosen consistently.

Use separation of variables with $u=W(x)T(t)$ and choose $W^{\prime \prime }/W=T^{\prime }/T=k$. Then $T=e^{kt}$ and $W^{\prime }(0)=-W(0)$, $W^{\prime }(1)=-W(1)$.

If $k=0$, write $W=A+Bx$. The first condition gives $B=-A$, and the second gives $A+2B=0$. Together they force $A=B=0$, so this case contributes only the zero solution.

If $k={\alpha }^2>0$, write $W=A\cosh (\alpha x)+B\sinh (\alpha x)$. The first condition gives $B=-A/\alpha$. The second becomesA nonzero solution therefore requires $\alpha =1$, giving $W\propto e^{-x}$ and time factor $e^t$.

If $k=-{\beta }^2<0$, write $W=A\cos (\beta x)+B\sin (\beta x)$. Now $B=-A/\beta$, and the remaining condition is $-A(\beta +{\beta }^{-1})\sin \beta =0$. Thus $\beta =n\pi$, $n\ge 1$. Choosing a convenient amplitude gives $W_n=n\pi \cos (n\pi x)-\sin (n\pi x)$. The complete Sturm-Liouville eigenfunction expansion is consequentlyThese are self-adjoint separated Robin boundary conditions: the boundary form $[f^{\ast }{g}^{\prime }-f^{\prime \ast }g{]}_0^1$ vanishes because both functions satisfy $y^{\prime }=-y$ at each end. The regular interval problem has a complete orthogonal eigenbasis; the calculation includes its one negative eigenvalue for $-D^2$, not only its positive ones.

To project the initial data, calculate the squared norms:The second identity uses the half-unit integrals of ${\cos }^2(n\pi x)$ and ${\sin }^2(n\pi x)$ and the zero integral of their product. By orthogonal projection,For square-integrable $f$, the initial expansion is understood in that norm. Classical agreement of boundary derivatives at $t=0$ requires compatible smoother data.

This is Robin diffusion with a growing mode, in the reflected version of its positive-coordinate-derivative example. Growth is not a contradiction: the outward derivative is $-u_x$ at the left endpoint, so these conditions are not dissipative outward-normal conditions at both ends. Indeedand the boundary can supply energy. The growing contribution is absent exactly when the displayed integral defining $c_0$ vanishes.

For real positive weight $w$, the weighted inner product isThe functions must have finite weighted norm; the conjugate is omitted for real functions. Multiplying the differential equation by $x{e}^{-x}$ givesThus the Sturm-Liouville form has $p=x{e}^{-x}$, $q=0$ and $w=e^{-x}$. With the positive-sign spectral convention defineThe existing Laguerre differential operator on polynomials uses the opposite operator $L=-\mathcal{A}$, with eigenvalues $-n$; keeping this sign convention explicit avoids confusing the two descriptions.

For admissible $f,g$, integration by parts givesThe boundary form must vanish for every pair in the self-adjoint domain. In the usual regular-at-zero Laguerre realization one imposes ${\lim }_{x\downarrow 0}x{e}^{-x}{y}^{\prime }=0$, excluding the logarithmic branch; at infinity one uses the weighted square-integrable domain with vanishing boundary form. A condition $y(\infty )=0$ would be inappropriate for polynomial eigenfunctions. For any pair of polynomials, the factor $x$ makes the boundary form zero at zero, while the exponential beats all polynomial growth at infinity. These functions also have finite weighted norms.

Self-adjointness yields $(n-m)⟨y_m,y_n{⟩}_w=0$, giving orthogonality of Sturm-Liouville eigenfunctions:For the specified polynomials, direct differentiation gives $\mathcal{A}u=u$ and $\mathcal{A}v=3v$. Hence their eigenvalues in the displayed convention are 1 and 3, respectively, or −1 and −3 for $L=-\mathcal{A}$.

To find the quadratic, put $h=a{x}^2+bx+c$. The equation $x{h}^{\prime \prime }+(1-x)h^{\prime }+2h=0$ reduces to $(4a+b)x+(b+2c)=0$. Thus $b=-4a$, $c=2a$. Since $h(1)=-a=1$,These are scalar multiples of the Laguerre polynomials of degrees one, three and two.

For $n=0$, the first integral is $x{e}^{-x}{y}^{\prime }=C$, so all solutions with the prescribed interior normalization areTwo independent choices areThey both equal one at $x_0$, and their Wronskian is $e^x/x\ne 0$. For sufficiently large $x$, the last unit interval contributes at least $e^{x-1}/x$, proving that $y_2\to \infty$. These are the zero-eigenvalue solutions of the Laguerre differential equation, with the second expressible through the exponential integral. It also has a logarithmic singularity at zero and is not an admissible polynomial or weighted square-integrable eigenfunction. An interior value alone does not enforce the endpoint conditions of the spectral problem.