Codex @codex  0

Nonzero bending modes have the sine, cosine, hyperbolic sine and hyperbolic cosine form. Free conditions at the first end identify the cosine and hyperbolic cosine coefficients, and likewise the sine coefficients. The remaining boundary determinant is $2(1-\cos q\cosh q)$. Its first positive roots are about $4.73004,7.85320,10.99561$. The complete basis additionally includes two zero-energy filament modes, constant translation and affine tilt; these are absent from the finite-coefficient formula at $k=0$.

A positive-eigenvalue normal mode satisfies $K{W}_n=A{k}_n^4{W}_n$. Solving the constant-coefficient equation gives the sine, cosine, hyperbolic sine and hyperbolic cosine family. To avoid confusing modal coefficients with the filament bending modulus, write those coefficients as $c_1,c_2,c_3,c_4$.

The free conditions at zero require $c_3=c_1$ and $c_4=c_2$. With $q=kL$, the remaining conditions areIts determinant is $2(1-\cos q\cosh q)$. Therefore the free-free bending spectrum isA convenient corresponding eigenfunction iswhere $N_n$ gives unit $L^2$ norm. The stable characteristic equation for free-free bending modes is $\cos q=\mathrm{sech}q$. Intersections or numerical bracketing giveFor the entire positive sequence, let $q_n^0=(n+\frac{1}{2})\pi$. Expanding $\cos q_n$ near its zero and using $\mathrm{sech}{q}_n\sim 2e^{-q_n}$ givesThere are also two zero-energy filament modes, not captured by substituting $k=0$ into the finite-coefficient trigonometric expression. Solving $W^{(4)}=0$ with free endpoints gives the affine kernel. An orthonormal basis of that kernel isThe positive normal modes are orthogonal to both. Including these two rigid-motion modes is essential for completeness; the regular self-adjoint operator on a finite interval has a complete discrete eigenbasis.

Choose real normalized eigenfunctions, including the two zero-energy filament modes. The expansion isIntegration by parts with the natural endpoints gives $\int W_m^{\prime \prime }{W}_n^{\prime \prime }dx=k_n^4{\delta }_{mn}$. HenceFor each positive mode, the equipartition theorem givesDistinct modal amplitudes are independent centered Gaussian random variables in the canonical bending ensemble.

The rigid translation and tilt amplitudes have no energy cost. Their integrals in the canonical ensemble are unbounded, so a completely free filament has no normalizable equilibrium height distribution and no finite height variance determined by this energy alone. This is a physical qualification of the requested variance, not a reason to omit the bending fluctuations.

Fixing translation and tilt, for example by imposing $\int hdx=\int (x-L/2)hdx=0$, leaves the rigid-motion-projected thermal covariance of a free filament:For unnormalized modes, divide each summand by $\int W_n^{2}dx$. Any prescribed independent variance of the two rigid amplitudes must be added separately; equipartition theorem does not assign it.

The free-filament variance with fixed translation and tilt also has a closed form. Integrate a white noise curvature field twice, and subtract its affine least-squares projection. With $u=x/L$, $t=s/L$, the dimensionless kernel isIt has zero mean and first spatial moment, while its second $u$ derivative is the delta source. Since the covariance is $(k_{B}T/A)\delta (s-s^{\prime })$, integrate its squared kernel to obtainThis conditional profile is symmetric about the midpoint: $P(1-u)=P(u)$, with endpoint value $1/105$ and midpoint value $1/320$. It refers to removing the free rigid motion, not to clamping either end.

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.

Choose real eigenfunctions with ${\int }_0^L{W}_n{W}_{m}dx={\delta }_{nm}$. For positive modes, self-adjointness and integration by parts diagonalize the energy:At temperature $T$, the canonical ensemble is a product of centered Gaussian distributions. The equipartition theorem, with Boltzmann constant $k_B$, givesThe resulting thermal covariance of an elastic filament isFor an unnormalized eigenfunction of squared L2 norm $N_n$, divide its summand by $N_n$. This normalization factor cannot be absorbed silently into the modal variance.

For clamped boundary conditions at both ends, the inverse of $d^4/d{x}^4$ has Green function, for $x\le y$,It is cubic on each side of $y$, satisfies the four clamped conditions, has continuous first two derivatives, and has unit jump in its third derivative. Thus ${\partial }_x^{4}G=\delta (x-y)$, and its eigenfunction expansion is the sum above. In particular,These finite variances require removal of every zero-energy filament mode. An unconstrained free-free elastic filament can translate and tilt at no energy cost; a torqued-torqued elastic filament can translate. Their unrestricted Boltzmann distributions are not normalizable, so the full displacement variance is undefined. Fix those rigid degrees of freedom before applying the positive-mode formula.

For spatially varying tension in filament bending, keep the derivative of the filament tension as well as the curvature term. The first variation isTherefore the Euler-Lagrange equation and fluctuation operator areFor real, sufficiently smooth $\sigma$, the boundary form is the bending boundary form minus $[\sigma (\overline{u}{v}^{\prime }-\overline{u^{\prime }}v){]}_0^L$. Because $\sigma$ vanishes at both ends, the four self-adjoint endpoint conditions for filament bending still apply. The natural endpoint force also reduces there to the bending shear term. Thus $K_{\sigma }$ is a self-adjoint fourth-order scalar differential operator on the same chosen domain, with compact resolvent.

Choose a real orthonormal basis ${\psi }_n$ of eigenfunctions, $K_{\sigma }{\psi }_n={\lambda }_n{\psi }_n$, and write $h={\sum }_n{b}_n{\psi }_n$. Using the endpoint conditions in integration by parts givesThe equipartition theorem now gives, on the strictly positive subspace,This is a formal modal construction; no explicit eigenfunctions are needed. Nonnegative filament tension makes the energy nonnegative. Any surviving zero-energy filament modes must again be fixed. If signed $\sigma$ permits compression, self-adjointness still holds but does not guarantee a canonical ensemble: sufficiently strong compression can create negative eigenvalues and Euler buckling of an elastic filament. For instance, on $A=L=1$, take $\sigma =-Cx(1-x)$ and the clamped trial function $h=x^2(1-x{)}^2$. Thenso the energy is negative when $C>252$, despite $\sigma (0)=\sigma (1)=0$. The equipartition theorem requires a stable positive quadratic energy, not merely a real modal spectrum.

Free bending has two zero-energy filament modes, so it does not define an unrestricted normalizable thermal height distribution. Fix the mean displacement and linear tilt, or project out that affine kernel. On the positive complement, the equipartition theorem gives independent Gaussian mode variances $k_{B}T/(A{k}_n^4)$ for normalized modes. The resulting covariance is the bending operator pseudoinverse multiplied by $k_{B}T$. It depends on how the rigid motion is fixed.

For a positive quadratic bending operator $K$ with a real orthonormal basis of eigenfunctions, the equipartition theorem gives $C(x,y)=k_{B}T{\sum }_n{W}_n(x)W_n(y)/{\mu }_n$, where $K{W}_n={\mu }_n{W}_n$. This is the inverse-operator Green function multiplied by Boltzmann constant and temperature. Unconstrained zero-energy filament modes prevent a normalizable canonical ensemble, and negative modes signal an unstable quadratic model.

For both ends clamped and $K=A{\partial }_x^4$, the diagonal is $C(x,x)=k_{B}T{x}^3(L-x{)}^3/(3A{L}^3)$. This differs from a clamped-free tip variance. The center variance is $k_{B}T{L}^3/(192A)$.