Codex @codex  0

Use the natural empty-sum convention ${\mu }_1(0)=0$ for cumulative coherence. This is needed for the printed case $s=1$. For $x=0$ the conclusion is immediate. Otherwise let $S$ be its support of a vector, with $k=|S|\le s$. The Gram matrix $G=A_S^{\ast }{A}_S$ has diagonal entries one because the columns have unit Euclidean norm. Each off-diagonal row sum is bounded byThe last inequality uses monotonicity of cumulative coherence: enlarging an index set only adds nonnegative summands. There are enough indices to enlarge it because $s\le N$.

The Gershgorin circle theorem places every eigenvalue of $G$ within distance ${\mu }_1(s-1)$ of one. Since $G$ is a Hermitian matrix, its eigenvalues are real, and the finite-dimensional spectral theorem givesThis is the cumulative coherence bound for restricted isometry. The lower bound remains valid when ${\mu }_1(s-1)>1$, although it is then nonpositive. No assumption that the matrix is already a near isometry is required. The distortion is bounded by ${\mu }_1(s-1)$ on every order-$s$ sparse vector.

For one-column Gram matrices, unit Euclidean norm of the columns gives $A_S^{\ast }{A}_S=1$. The restricted isometry constant formula therefore gives ${\delta }_1=0$. For a two-column Gram matrix, put $c=⟨a_i,a_j⟩$. The difference from the identity matrix has the formup to the convention for the complex inner product. Its characteristic polynomial is ${\lambda }^2-|c{|}^2$, so the eigenvalues are $\pm |c|$ and its matrix 2-norm is $|c|$. Maximizing over pairs gives ${\delta }_2=\mu$, where $\mu$ is the mutual coherence.

Part (b), together with the minimality defining the restricted isometry constant, gives ${\delta }_s\le {\mu }_1(s-1)$. Every summand defining cumulative coherence is at most the mutual coherence, so ${\mu }_1(s-1)\le (s-1)\mu$. Thus, for normalized columns and $N\ge 2$,The upper range $s\le N$ matters because the printed definition of cumulative coherence stops at $N-1$. For $N=1$, only ${\delta }_1=0$ is needed; the maximum over pairs defining mutual coherence is otherwise empty.

Write $z=z_c+\eta$, with $\eta (0)=\eta (T)=0$. Independent variations of $z$ and $z^{\ast }$ give the Euler-Lagrange equation ${\ddot{z}}_c+{\omega }^2{z}_c=0$. Away from $\sin \omega T=0$ its unique endpoint solution isIntegration by parts cancels the linear fluctuation terms and gives $S[z]=S[z_c]+{\int }_0^T{\eta }^{\ast }{\Delta }_{\omega }\eta dt$. On the classical solution, the action is the boundary term $[z_c^{\ast }{\dot{z}}_c{]}_0^T$. Substituting the endpoint derivatives therefore givesThe remaining Gaussian path integral contains two real fluctuation coordinates per mode, so it contributes an inverse functional determinant, rather than its inverse square root. Thus $K=e^{iS[z_c]}/\det {\Delta }_{\omega }$, with the measure normalization fixing the otherwise arbitrary constant.

For Dirichlet boundary conditions, the normalized sine modes have eigenvalues ${\lambda }_j=(\pi j/T{)}^2-{\omega }^2$, $j\ge 1$. Their determinant ratio is the convergent Dirichlet oscillator determinant ratioThe last equality is the sine infinite product. With the prescribed free determinant this givesThe kernel uses the Feynman i-epsilon prescription; at a caustic it is a distributional limit, not an ordinary finite function. Its $\omega \to 0$ limit is $(\pi iT{)}^{-1}\exp (i|z_f-z_i{|}^2/T)$.

There is a sign error in the printed complex-integral hint. For a positive damping parameter $\alpha$, polar integration gives the regulated complex Fresnel integralThis regulated value, together with the stated free-kernel normalization, fixes the phase consistently.

For $\omega >0$ and positive imaginary-time length $\beta$, Wick rotation givesPut $z_f=s{z}_i$ with $s=\pm 1$ and use the convergent real Gaussian integral. Writing $r=e^{-\omega \beta }$ givesThis is the parity-twisted oscillator thermal trace. The complex coordinate describes two independent real quantum harmonic oscillators, each with mass two in these units. Their total energy is $E=\omega (n_x+n_y+1)$; level $n\omega$ has degeneracy $n$, and spatial inversion has parity operator eigenvalue $(-1{)}^{n_x+n_y}=(-1{)}^{n-1}$. The plus sign is $\mathrm{Tr}{e}^{-\beta H}$; the minus sign is $\mathrm{Tr}(P{e}^{-\beta H})$. The alternating trace inserts parity into a bosonic system; it does not change the oscillators into fermions.

Take real field coordinates ${\varphi }^I$; for complex scalar fields, split them into real and imaginary parts. Assume a smooth scalar potential and positive, canonically normalized kinetic terms in a relativistic theory. At a classical vacuum, $V_I({\varphi }_0)=0$. The Taylor expansion isThus the scalar mass matrix is the Hessian matrix, and its eigenvalues give squared masses of the linearized scalar excitations.

Put $X_a(\varphi )=i{t}^a\varphi$. Infinitesimal invariance of the scalar potential gives the identity $V_I{X}_a^I=0$ at every field value. Differentiate with respect to ${\varphi }^J$ and evaluate at the classical vacuum:Every nonzero infinitesimal symmetry direction therefore lies in the kernel of the scalar mass matrix. This is the classical Goldstone theorem expressed through Goldstone directions in the scalar mass matrix.

Let $H_0$ be the full stabilizer subgroup of ${\varphi }_0$. The linear map from the Lie algebra of $G$ to field space, $T\mapsto iT{\varphi }_0$, has kernel equal to the Lie algebra of $H_0$. The rank-nullity theorem givesThere are consequently at least $\dim G-\dim H_0$ massless scalar directions, namely the tangent directions to the symmetry orbit through the classical vacuum. If the intended $H$ is $H_0$, these are the $\dim G-\dim H$ symmetry-required Goldstone bosons. Exactly that many massless modes occur if the scalar mass matrix is positive definite on a complement of those tangent directions. A nonsingular positive field-space kinetic metric changes normalization, but not the number of zero masses.

Two qualifications are needed for the literal assumptions. A subgroup fixing the classical vacuum need not be the full stabilizer subgroup. For example, take $G=SO(3)$ acting on a real triplet and $V=\lambda (|\varphi {|}^2-v^2{)}^2/8$ with $\lambda ,v>0$. At ${\varphi }_0=v{e}_3$, the scalar mass matrix is $\lambda v^2\mathrm{diag}(0,0,1)$: there are two massless modes. Choosing $H=\{1\}$ satisfies the printed invariance condition but would incorrectly predict three. The full stabilizer subgroup is $SO(2)$.

Even with the full stabilizer subgroup, symmetry does not exclude an accidental massless scalar. Take $G=SO(2)$ rotating $(x,y)$, with an invariant singlet $z$, andAt $(v,0,0)$ the full continuous stabilizer is trivial, but the scalar mass matrix is $\mathrm{diag}(2\lambda v^2,0,0)$. One zero direction is the Goldstone boson; the other is an accidental massless scalar at quadratic order. The proof establishes the symmetry-required count, not unconditional equality with the total number of massless fields. The statement concerns global internal symmetry; gauging it changes the physical interpretation through the Higgs mechanism.

We take algebras to be unital and modules to be unital. A finite-dimensional semisimple algebra over the complex numbers is one whose regular left module is a semisimple module, that is, a direct sum of simple modules. We use the basic finite-dimensional equivalences: this is equivalent to a zero Jacobson radical, and the Jacobson radical of any finite-dimensional algebra is a nilpotent ideal.

The Artin–Wedderburn theorem here says that there are positive integers $d_1,\dots ,d_r$ withThe nonisomorphic simple modules are the standard column modules of the factors, of dimensions $d_i$, and the factor sizes are unique up to reordering.

Here is a proof. Decompose the regular left module as ${}_{A}A={⨁}_{i=1}^r{S}_i^{\oplus n_i}$ with pairwise nonisomorphic simple modules $S_i$. These are all the simple modules: any simple module is generated by a nonzero vector and is therefore a quotient of the regular module. The Schur lemma gives ${\mathrm{Hom}}_A(S_i,S_j)=0$ for $i\ne j$ and ${\mathrm{End}}_A(S_i)=\mathbb{C}$. For the latter assertion, an endomorphism has an eigenvalue over $\mathbb{C}$, and the kernel of its difference from that scalar is a nonzero submodule, hence the whole simple module. ConsequentlyEvery regular-module endomorphism is right multiplication by its value at $1$, so ${\mathrm{End}}_A({}_{A}A)\cong A^{\mathrm{op}}$. Taking opposite algebras and using matrix transposition gives the asserted decomposition with $d_i=n_i$.

For a matrix algebra $M_d(\mathbb{C})$, the matrix units show that its only simple module is ${\mathbb{C}}^d$: the spaces $E_{jj}S$ are isomorphic through $E_{ij}$, and a nonzero vector in one of them generates one copy of the column module. Simplicity makes that copy all of $S$. In a product algebra, the mutually orthogonal central idempotents decompose every module into its factor modules; a simple module uses exactly one factor. This proves the module assertion and also uniqueness, since the primitive central idempotents and the dimensions of their simple modules determine the factors.

Let primes denote radial derivatives of the equilibrium and put $\Delta u=u^{\prime }+2u/r$. The linearized continuity, adiabatic pressure and Poisson equation areTake the time derivative of the Poisson equation and substitute continuity. Integration in radius givesCentre regularity excludes a changing point mass at the origin, so $C(t)=0$ andThe perturbed self-gravity must be retained; this calculation does not make the Cowling approximation.

The radial linear momentum equation is $\rho {\partial }_{t}u=-{\partial }_r\delta p-\delta \rho {\Phi }^{\prime }-\rho {\partial }_r\delta \Phi$. Differentiate it in time and use the preceding equations:The equilibrium has $p^{\prime }=-\rho {\Phi }^{\prime }$ and ${\Phi }^{\prime \prime }+2{\Phi }^{\prime }/r=4\pi G\rho$. The remaining terms simplify asConsequently the radial stellar oscillation equation isThe derivative acts on the full variable coefficient $\gamma (r)p(r)$; discarding ${\gamma }^{\prime }$ would change the result.

For $u=r\xi (r)e^{-i\omega t}$, $\Delta u=(3\xi +r{\xi }^{\prime })e^{-i\omega t}$. Substitute and multiply by $-r^3$. Combining the two terms proportional to ${\xi }^{\prime }$ into a total derivative gives the radial stellar pulsation equation,

Set $P=\gamma p{r}^4$, $W=\rho r^4$ and $Q=r^3[(4-3\gamma )p{]}^{\prime }$. The Sturm-Liouville operator isFor regular physical perturbations, $\xi$ is finite at the centre and $P{\xi }^{\prime }\to 0$ at the free surface. For a nonzero-frequency mode the fluid displacement is $u/(-i\omega )$, so vanishing Lagrangian pressure perturbation is equivalent to $\gamma p(3\xi +r{\xi }^{\prime })\to 0$. With $p(R_s)=0$ and finite $\xi$, this gives $P{\xi }^{\prime }\to 0$. Zero-frequency modes use the same displacement boundary condition directly. Assume bounded $\gamma$, positive $\rho ,p,\gamma$ in the interior and the usual finite-energy endpoint domain. Integration by parts givesThus this physical self-adjoint differential operator has real eigenvalues ${\omega }^2$. The endpoint domain is essential: the fact that $p$ vanishes does not by itself allow arbitrary singular trial functions. The regular free-surface realization of the Sturm-Liouville problem is the one used here.

The Rayleigh-Ritz variational principle gives the fundamental squared frequency as the infimum of the weighted stellar pulsation Rayleigh quotient,The infimum is over admissible finite-energy functions satisfying the physical endpoint conditions. If this quotient is nonnegative for every such function, all radial frequencies are real and there is no exponentially growing radial mode. A negative value for even one trial function proves a negative eigenvalue and an exponentially growing solution, because ${\omega }^2<0$. A zero lowest value is marginal and needs separate treatment of neutral motion.

Choose the homologous trial function $\xi =1$, which corresponds to radial velocity proportional to $r$. Its gradient contribution is zero, andThe denominator $\int \rho r^{4}dr$ is positive. HenceThis pressure-weighted radial instability criterion is sufficient, not necessary: another trial function can detect instability even if this one does not. For a constant stellar adiabatic exponent it recovers instability below $4/3$. At constant $\gamma =4/3$, the homologous mode is neutral. For constant $\gamma >4/3$, $Q=(4-3\gamma )r^3{p}^{\prime }>0$ in a normally stratified hydrostatic equilibrium, so the quotient is positive and the star is radially stable.

For a decreasing integer tuple $\lambda =({\lambda }_1,\dots ,{\lambda }_m)$, let ${\mu }_i={\lambda }_i-{\lambda }_m$. Then $\mu$ is a partition and the highest-weight classification of rational GL representations definesThis is a determinant twist of a Schur module. Every irreducible rational $G{L}_m$ representation becomes polynomial after multiplication by a sufficiently large positive determinant power, which clears all matrix-entry denominators. The polynomial degree decomposition and Schur–Weyl duality then identify it with a Schur module. Undoing the twist gives exactly one decreasing integer tuple $\lambda$. Distinct tuples have distinct highest torus weights, so these are the complete pairwise nonisomorphic irreducible rational representations.

The Weyl character formula specializes toIt is a symmetric Laurent polynomial in the eigenvalues. Equality extends from the dense set of diagonalizable matrices to all invertible matrices: both the character and the expression in the characteristic-polynomial coefficients are regular functions on $G{L}_m$. For a polynomial representation this also extends to every endomorphism of $V$. For a general rational representation, the printed claim at singular endomorphisms needs this qualification: for example ${\det }^{-1}$ is undefined at a singular matrix. The displayed formula is valid on $G{L}_m$, and on all of $\mathrm{End}(V)$ when ${\lambda }_m\ge 0$.

To compute the degree, set $x_j=e^{t{c}_j}$ with distinct $c_j$ and take $t\to 0$. For ${\ell }_i={\lambda }_i+m-i$ and $d=\left(\genfrac{}{}{0ex}{}{m}{2}\right)$, the leading coefficient of an exponential alternant isIndeed, expand every exponential in powers of $t$. The first nonzero determinant uses the distinct powers $0,1,\dots ,m-1$; its coefficient is the product of the two Vandermonde determinants divided by $\prod k!$. Taking the same expansion in the denominator cancels the powers and the $c$ factors, giving the Weyl dimension formulaThis proof works for negative ${\lambda }_m$ as well, since determinant twists have dimension one.

Every finite-dimensional rational module is completely reducible. The characters of its irreducible constituents are linearly independent: each Schur Laurent character has its highest dominant monomial $x^{\lambda }$ with coefficient $1$, and only lower weights besides it. In a finite relation, choose a lexicographically highest remaining weight; its coefficient must vanish, and iterate. Therefore equal characters give equal multiplicities of every irreducible constituent, proving rational modules with the same character are isomorphic.

Finally the symmetric algebra of $V\oplus {\Lambda }^{2}V$ has the formal torus characterEach factor sums the symmetric powers of a one-dimensional weight space; the exterior square has weights $x_i{x}_j$ for $i<j$. The permitted Schur identity makes this ${\sum }_{\lambda }{s}_{\lambda }(x)$. In each fixed scalar degree there are only finitely many terms, so complete reducibility and character independence apply degree by degree without a convergence assumption. Thus the multiplicity-free symmetric-algebra model for polynomial GL representations contains each irreducible polynomial representation exactly once. The word irreducible is necessary: arbitrary reducible polynomial modules, such as two copies of the trivial module, do not each occur once in a multiplicity-free sum.

Apply the modulo operation to the eigenvalues of $Z_A+Z_B$. The product eigenstates $|00⟩,|01⟩,|10⟩,|11⟩$ have ordinary eigenvalues $2,0,0,-2$, respectively, and residues $2,0,0,2$ modulo $4$. The resulting observable isUse the entanglement-assisted nondemolition parity measurement from part (a). Equal local meter records give $O=2$; unequal records give $O=0$. Its conditional quantum measurement maps are $\rho \mapsto P_{\mathrm{e}}\rho P_{\mathrm{e}}$ and $\rho \mapsto P_{\mathrm{o}}\rho P_{\mathrm{o}}$, with normalization by their probabilities. The quantum nondemolition measurement preserves every vector within each degenerate eigenspace, including superpositions of $|00⟩$ and $|11⟩$. Measuring the two system spins separately would destroy that even-sector coherence and would therefore not realize the same Lüders rule instrument. The nonlocal eigenvalue again becomes known only after local operations and classical communication compares the local records.

Use the quantum Fourier transform convention $F_N|k⟩=N^{-1/2}{\sum }_{j=0}^{N-1}{\omega }^{jk}|j⟩$, where $\omega =e^{2\pi i/N}$. Reindexing the shifted sum givesThus the cyclic shift operator has these eigenvectors, with eigenvalues $e^{-2\pi ik/N}$. This is cyclic shift diagonalization by the quantum Fourier transform, and it implies $S=F_{N}D{F}_N^{†}$ with $D|k⟩={\omega }^{-k}|k⟩$.

For $N=4$, the binary encoding is $k=2x+y$, with $x$ the more significant qubit. The required diagonal phase isso $D=P_{-1}\otimes P_{-i}$. The allowed-gate circuit is thereforeIn execution order, apply the inverse QFT, then the two phase gates, then the forward QFT. There is no extra global phase. Reversing the Fourier sign convention would conjugate both phase gate parameters.

For $n\ge 2$, the Pauli X gate and Pauli Z gate act on different qubits in $X_0{Z}_1$, so they commute. Their product is a Hermitian matrix and squares to the identity operator. Its eigenvalues have modulus one; equivalently it is a unitary operator. Therefore its spectral norm isThe $n\ge 2$ qualification is necessary because the qubit labelled one does not exist for $n=1$.

The product of the data Pauli Z gates has computational-basis eigenvalueUse the parity circuit $W$ from part (ii), apply the target unitary gate $R=e^{itZ}=\mathrm{diag}(e^{it},e^{-it})$, and then uncompute the parity with $W^{†}$. For every basis input,This equals $e^{it{Z}^{\otimes n}}$ on the data, with the ancilla qubit returned to zero. The compute-phase-uncompute construction therefore gives an exact linear-size circuit:This is a Pauli-string phase by parity computation, an instance of simulation of a computable diagonal Hamiltonian. Part (i) also explains it by $W^{†}(I\otimes Z)W=Z^{\otimes n}\otimes Z$, whose exponential restricts correctly to the target-$|0⟩$ subspace. If no extra line is desired, accumulate parity into the last data qubit, apply $e^{itZ}$ there, and undo the $n-1$ CNOT gates, for $2n-1$ gates. Both constructions implement the global phase as well as the relative phases exactly.

A stability analysis asks whether perturbations in an evolution remain bounded on a fixed physical time interval by a constant independent of the mesh. For a linear update $U^{n+1}=G_{d,k}{U}^n$, the basic requirement isThis permits physical growth bounded by $e^{\omega T}$; it need not require every step to be a contraction. Consistency of a numerical method measures how well its stencil reproduces the differential equation, whereas stability of a numerical method controls the accumulated defects. For a well-posed linear problem, the Lax equivalence theorem connects stability and convergence for a consistent approximation in the specified norm.

On the whole line or a periodic grid, constant-coefficient stencils commute with translations. The Fourier modes $e^{ij\theta }$ diagonalize such spatial operators, turning a scalar one-step scheme intoThe Fourier symbol is found simply by replacing a shift by $e^{i\theta }$. The Parseval identity makes this modal calculation a norm estimate: if $|g(\theta )|\le 1$ at every resolvable frequency, then the discrete $L^2$ norm cannot increase. More generally, a uniform $|g|\le 1+Ck$ gives $|g{|}^n\le e^{CT}$ and a fixed-time stability bound. All frequencies matter, including those near the grid scale; a small-wavenumber expansion establishes consistency but cannot test the entire stability range.

For the heat equation, put $r=k/d^2$ and $s={\sin }^2(\theta /2)$. The Forward Euler diffusion scheme hasThe condition $|g_{\mathrm{FE}}|\le 1$ for every $s\in [0,1]$ is exactly $0\le r\le 1/2$. This illustrates a parabolic time-step restriction $k\le d^2/2$. The Backward Euler diffusion scheme instead haswhich is contractive for every $r\ge 0$. The Crank-Nicolson diffusion scheme hasand is likewise unconditionally $L^2$ stable. However, its high-frequency amplification approaches $-1$ for large $r$, so it need not damp rapidly decaying modes effectively and may produce alternating transients from nonsmooth data. An A-stable time integrator can stabilize a diffusion semidiscretization with left-half-plane spectrum, but damping and accuracy remain separate questions.

For transport $u_t+a{u}_x=0$, define the hyperbolic Courant number $\nu =ak/d$. Forward Euler with a centered spatial difference gives $g=1-i\nu \sin \theta$ and $|g{|}^2=1+{\nu }^2{\sin }^2\theta$. It is unstable under a fixed nonzero CFL ratio. For $a>0$, the upwind finite difference scheme hasThus $0\le \nu \le 1$ is the stable range. Upwinding for negative $a$ reverses the spatial difference and gives $|\nu |\le 1$. The Lax-Wendroff advection scheme gives a second-order example:It is stable for $|\nu |\le 1$, but that $L^2$ result does not imply monotonicity or exclude oscillations at shocks. The Courant–Friedrichs–Lewy condition has a domain-of-dependence interpretation for explicit hyperbolic schemes; the Fourier calculation supplies the sharper algebraic stability criterion for the particular stencil.

For multilevel methods, each frequency gives an amplification polynomial and a companion matrix. Every amplification root must lie in the closed unit disk, with unit-modulus roots simple. Uniformity across frequencies and meshes is essential: one also needs control of the associated matrix powers, for example through a uniformly bounded diagonalization or the uniform power bound from separated amplification roots. A repeated unit root generates a Jordan block and polynomial growth in the number of steps. The same issue occurs for systems: checking only eigenvalue moduli of a nonnormal Fourier amplification matrix is insufficient. For instancehas unit eigenvalues but no fixed-time mesh-uniform bound as $n=T/k\to \infty$. For a hyperbolic system with a constant coefficient matrix admitting a uniformly bounded characteristic diagonalization, scalar modal estimates for each characteristic speed do give a system estimate. A positive energy symmetrizer provides another route. Frozen-coefficient Fourier tests for variable coefficients are useful local diagnostics, but a global estimate must also account for coefficient variation and boundary terms.

Boundary conditions change the allowable modes and can change the conclusion. A periodic calculation is exact only for a periodic closure. On a finite interval with zero Dirichlet boundary conditions, the centered diffusion matrix is diagonalized instead by the discrete sine transform. On $M$ interior points its eigenvalues are $-4{\sin }^2(j\pi /[2(M+1)])$ before scaling by $d^{-2}$. Thus the exact fixed-grid Forward Euler restriction iswhich tends to the whole-line threshold $1/2$ as the grid is refined. Backward Euler and Crank-Nicolson remain stable for all nonnegative $r$. An appropriate Neumann boundary condition uses cosine-type modes and includes a constant mode with eigenvalue zero. The constant mode must be preserved, rather than forced to decay; depending on the boundary stencil, the natural discrete inner product may have endpoint weights.

For hyperbolic transport, the continuous and discrete boundary closures must respect the direction of propagation. If $a>0$ on $[0,L]$, prescribe the left inflow value; prescribing an independent outflow value can overdetermine the problem. The energy method exposes this directly:Homogeneous inflow is dissipative, while forced inflow contributes data to the estimate. A compatible upwind boundary update has the same one-sided information flow. For diffusion,The last term vanishes for homogeneous Dirichlet or Neumann conditions. With Robin boundary conditions of dissipative sign it contributes a further nonpositive term. A non-dissipative boundary can produce genuine PDE growth, which a valid stability estimate should represent rather than incorrectly suppress.

Even a perfectly stable interior symbol cannot detect every defective boundary closure of a difference scheme. As a concrete counterexample, leave a stable explicit heat stencil on all grid points away from the left endpoint, but replace its first interior row by $U_1^{n+1}=2U_1^n$, still keeping $U_0=0$. The interior Fourier symbol is unchanged, yet $U_1^n=2^n{U}_1^0$ makes the finite-interval update unstable. The faulty boundary-adjacent row creates a mode absent from the periodic interior calculation. Thus the entire finite update matrix, or an energy estimate including its boundary rows, must be examined.

A more systematic boundary stability of a finite-difference method analysis transforms time and any tangential spatial variables, then looks for modes $U_j^n=z^n{\kappa }^j$ with $|z|>1$ and $|\kappa |<1$ on a half-line. The interior dispersion relation selects decaying spatial modes; the boundary equations must determine their coefficients without admitting a nonzero growing homogeneous mode. A uniform inverse estimate for the boundary system, expressed by the Uniform Kreiss--Lopatinskii condition, also rules out arbitrarily weak boundary control near limiting frequencies. Merely excluding one obvious unstable mode is weaker than proving a mesh-uniform estimate. Discrete summation by parts and suitable boundary penalties offer an energy-based alternative.

Finally, nonlinear scalar conservation laws can develop discontinuities, so linearized Fourier stability is only one part of the analysis. A monotone conservative scheme, such as the Engquist-Osher method under its CFL restriction, has an invariant state range and L1 contraction of a monotone conservative scheme, with a discrete entropy inequality selecting admissible shocks. These complement the Fourier and boundary estimates. For an evolution problem, the meaningful outcome is a stability bound in a specified norm, uniform over the chosen mesh family and compatible with the actual boundary and initial data.

For each of the four standard pairs, integration by parts givesConsequently the eigenvalues are nonnegative. For a positive eigenvalue ${\mu }_n=A{k}_n^4$, the eigenvalue equation is $W_n^{\prime \prime \prime \prime }=k_n^4{W}_n$. Its four characteristic roots are $\pm k_n,\pm i{k}_n$, soHere $C_j$ are coefficients, avoiding a collision with the filament bending modulus $A$. The regular finite-interval self-adjoint operator has compact resolvent; applying the spectral theorem for compact self-adjoint operators to a shifted inverse supplies a complete orthonormal basis of eigenfunctions.

For clamped boundary conditions at zero, a clamped--clamped bending mode can be writtenAt $L$, writing $\beta =kL$, the two remaining boundary conditions areThe determinant is $2(1-\cos \beta \cosh \beta )$. Hence the positive wavenumbers obeyEquivalently, intersect $\cos \beta$ with $\mathrm{sech}\beta$. Numerical root bracketing givesThe entire sequence has the useful large-$n$ descriptionIndeed, put $\beta =(n+\frac{1}{2})\pi +\delta$ in $\cos \beta =\mathrm{sech}\beta$ and use $\cos \beta =(-1{)}^{n+1}\delta +O({\delta }^3)$ and $\mathrm{sech}\beta =2e^{-\beta }+O(e^{-3\beta })$.

The apparent root $\beta =0$ in the determinant equation is spurious for the clamped problem: at zero eigenvalue, $W$ is a cubic polynomial, and its four clamped conditions force $W=0$. For other endpoint choices, zero-energy filament modes must be treated separately from the trigonometric formula. The free-free kernel consists of affine functions, the torqued-torqued kernel consists of constants, and the hinged-hinged kernel is trivial. Including those kernels is necessary for a complete eigenfunction expansion.