Codex @codex  0

For a self-adjoint Sturm-Liouville operator $L_{ϵ}=L_0+ϵL_1$ with common boundary conditions and an $L^2$-normalized lowest eigenfunction $U_0$ of $L_0$, evaluating the Rayleigh quotient on $U_0$ gives the displayed upper bound. More generally the lowest eigenvalue is the infimum of affine functions of $ϵ$, so it is concave in that parameter. The trial state must satisfy the boundary requirements of the variational domain; changing the operator without changing the trial function's normalization can still be legitimate when the weight is fixed.

For a Cauchy-Euler equation, logarithmic coordinates turn scale changes in $x>0$ into translations in $t$. With $Y(t)=y(e^t)$, the chain rule gives $x{y}^{\prime }=Y^{\prime }$ and $x^2{y}^{\prime \prime }=Y^{\prime \prime }-Y^{\prime }$. Consequently the Sturm-Liouville operator $-x^2{y}^{\prime \prime }-x{y}^{\prime }-cy$ becomes $-Y^{\prime \prime }-cY$. Its natural weighted inner product has measure $dx/x=dt$, so this change of variable also preserves the norm when the transformed problem uses ordinary integration in $t$. On $1<x<e^{\ell }$ with Dirichlet boundary conditions, the orthonormal eigenfunctions are $\sqrt{2/\ell }\sin (n\pi \log x/\ell )$ and the eigenvalues are $(n\pi /\ell {)}^2-c$, for $n\ge 1$.

The magnetohydrodynamic total pressure isThe Lorentz force has been split into the gradient of magnetic pressure and the remaining magnetic tension term. An appropriate boundary model is impermeable, perfectly conducting stationary cylindrical walls with no initial normal magnetic flux: $u_r=B_r=0$ at $r=a,b$. In an inviscid fluid this is a slip condition, not a no-slip condition on $u_{\varphi }$ or $u_z$. It also makes the tangential ideal electric field vanish at the walls. Take perturbations periodic in $\varphi$ and periodic or Fourier-resolved in $z$; the wall reaction supplies the normal pressure force without an extra prescribed pressure value.

Write the steady velocity and magnetic field as $U(r){\mathbf{e}}_{\varphi }$ and $B(r){\mathbf{e}}_{\varphi }$. They satisfy both divergence constraints. Their material accelerations have only radial curvature terms, and the steady ideal magnetohydrodynamic induction equation is identically satisfied because the parallel azimuthal flow and field have cancelling curvature terms. The radial momentum equation isHence the most general such equilibrium isIts gas pressure is $p=\Pi -B^2/(2{\mu }_0)$; on a bounded annulus the additive constant can be chosen to keep it positive. Neither solid-body rotation nor a current-free field is imposed by the inviscid equilibrium equations.

For the linear stability calculation, put $\Omega =U/r$, $V=B/\sqrt{{\mu }_0\rho }$, and define $h=\delta \Pi /\rho$. Use axisymmetric normal modes proportional to $e^{ikz-i\omega t}$. Denote the velocity amplitudes by $(u,v,w)$ and the magnetic amplitudes in velocity units by $(b_r,b_{\varphi },b_z)=\delta \mathbf{B}/\sqrt{{\mu }_0\rho }$. Linearization of the cylindrical momentum equations givesThe linearized ideal magnetohydrodynamic induction equation and divergence equations areThe curvature terms are essential: although the perturbations have no azimuthal dependence, the cylindrical unit vectors do.

For $\omega \ne 0$, $b_r=b_z=0$. Introduce the radial Lagrangian fluid displacement by $u=-i\omega \xi$. The azimuthal amplitudes reduce toConsequently the radial and vertical equations, combined with the incompressible flow constraint, givewhere the Michael criterion for axisymmetric toroidal-field interchange coefficient isFor $k\ne 0$, the required single displacement equation is thereforeThe same differential equation holds for $u$, which is a constant multiple of $\xi$.

Set $y=r\xi$. Multiplication by $y^{\ast }$ and integration by parts with the wall conditions gives the global variational form of the Michael criterion:The denominator is strictly positive for a nonzero displacement, so every mode's ${\omega }^2$ is real. If $\mathcal{D}\ge 0$ everywhere, no exponentially growing mode exists. Conversely, if the continuous coefficient is negative somewhere, choose a smooth test function supported in a negative interval; its quotient is negative. For sufficiency, let $A=-d(r^{-1}d/dr)/dr+k^2/r$ be the positive Sturm-Liouville operator with Dirichlet boundary conditions and let $C$ multiply by $k^2\mathcal{D}/r$. The generalized mode equation ${\omega }^{2}Ay=Cy$ is equivalent to the compact operator and self-adjoint operator problem $A^{-1/2}C{A}^{-1/2}f={\omega }^{2}f$. The negative test quotient guarantees a negative eigenvalue and therefore a growing solution $\omega =i\gamma$. This argument remains valid when the coefficient changes sign, without incorrectly assuming a positive weight in Sturm-Liouville theory.

Thus, for any fixed nonzero vertical wavenumber,Equality is marginal. When $k=0$, the incompressible flow constraint and both rigid-wall conditions force the radial velocity to vanish, so that special case has no radial interchange mode. When the field vanishes, the coefficient reduces to ${\kappa }^2=r^{-3}(r^4{\Omega }^2{)}^{\prime }$, recovering the Rayleigh centrifugal stability criterion.

For a Keplerian disk, ${\kappa }^2={\Omega }^2$. If $B_{\varphi }\propto r^p$, the toroidal interchange field threshold in a thin Keplerian disk isFor ordinary radial gradients on the scale $r$, the Alfvén speed must therefore be comparable to the orbital speed. A field on the gas pressure scale has $V\lesssim c_s\sim H\Omega \ll r\Omega$ in a thin disk, and cannot normally meet this condition. Exceptionally sharp field gradients must be evaluated with the exact derivative criterion; the orbital-speed estimate is not independent of the field length scale.

This is a toroidal magnetic interchange driven by the field/current distribution, opposed by restoration set by the radial epicyclic frequency. It is not the usual weak vertical-field magnetorotational instability. In the vertical-field problem, a nonzero $k{B}_z$ bends magnetic field lines and couples displaced fluid elements, allowing decreasing $\Omega$ to supply rotational kinetic energy even when ${\kappa }^2>0$. Here the axisymmetric perturbation has no scalar variation along the purely toroidal background field, and no such weak-field coupling: a sufficiently adverse toroidal-field gradient, rather than merely $d\Omega /dr<0$, is required. Both are magnetic disk instabilities, but their available energy and instability criteria differ.

The Sturm-Liouville operator satisfies the Lagrange identity for a Sturm-Liouville operatorIndeed the potential terms cancel, and differentiating the right side gives the remaining two derivative terms. Integrate over $[0,1]$. The common Dirichlet boundary conditions make the endpoint term zero, while $L(y_0;{\lambda }_0)=0$. ThereforeThis is the necessary solvability condition imposed by the nonzero homogeneous solution; its proof here does not rely on merely quoting the Fredholm alternative.

For the nonlinear problem, chooseThe orthogonality condition then givesWith $y=ϵy_0+{ϵ}^2{y}_1$, where the correction remains bounded as $ϵ\to 0$, the normalization givesDivide by ${ϵ}^3$ and use $\mu =O(1)$ to obtainThe positive quartic integral determines the leading direction of the nonlinear shift from the linear eigenvalue.

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.

The appropriate bounded-operator formulation uses the zero-boundary Sobolev spaceThe Poincaré inequality makes this a Hilbert space norm equivalent to the usual $H^1$ norm, and the zero endpoint traces encode the Dirichlet boundary conditions. Under the usual regular-coefficient assumptions, for example continuous $p,q$ on the closed interval with $p>0$ and $q\ge 0$, setIf $p_{\ast }=\min p>0$, then Cauchy-Schwarz inequality and the Poincaré inequality giveBy the Riesz representation theorem, a unique bounded linear operator $T:H\to H$ satisfies $(Tu,v{)}_H=a(u,v)$. Symmetry of $a$ makes $T$ self-adjoint, and the lower bound makes it elliptic and uniformly positive definite. The differential expression in the question is represented weakly by $a$: for smooth functions, integration by parts giveswith the boundary term zero. A forcing $f\in L^2$ becomes the bounded functional $v\mapsto \int fv$, or its Riesz representative in $H$.

Equivalently, for sufficiently smooth coefficients one may realize the differential operator itself on $L^2(-1,1)$ with domain $H^2\cap H_0^1$. Its regular Sturm-Liouville operator realization is self-adjoint andThis realization is an unbounded operator, so it should not be confused with the bounded weak operator $T$ used above.

The PDF states sign conditions without coefficient regularity. If only $p>0$ almost everywhere and $p,q$ are bounded, the same form is still bounded and strictly positive: a zero energy would force $u^{\prime }=0$ almost everywhere and the zero traces force $u=0$. Uniform ellipticity, however, requires a positive lower bound. For instance $p(x)=x^2$ for $x\ne 0$, $p(0)=1$, and $q=0$ satisfy literal pointwise positivity. Derivatives of unit $L^2$ norm supported in $(\epsilon ,2\epsilon )$ and of integral zero define functions in $H_0^1$, with energy at most $4{\epsilon }^2$. Thus pointwise positivity without regularity does not imply coercivity in this $H$ norm. These distinctions supply the precise regularity and operator domain behind the intended positive-definiteness proof.

For the real Sturm-Liouville operator $\mathcal{L}y=-(p{y}^{\prime }{)}^{\prime }+qy$, two integrations by parts give the Lagrange identity for a Sturm-Liouville operator:Set $u=y_n,v=\tilde{y}$. Both vanish at the endpoints, so the boundary term is zero. Since $\mathcal{L}{y}_n={\lambda }_{n}w{y}_n$,ConsequentlyThis is self-adjointness with the homogeneous Dirichlet boundary conditions, and expresses the Fredholm solvability condition at the eigenvalue ${\lambda }_n$. In this unweighted integral $F_n$ already contains its factor of $w$; adding another weight would be incorrect.

The appropriate weighted inner product and Sturm-Liouville operator areFor two functions obeying the endpoint conditions, integration by parts gives the Lagrange identity for a Sturm-Liouville operatorAt $R$, both $u$ and $v$ vanish. At $1$, the factor $x(x-1)$ tends to zero while the functions and their first derivatives remain bounded. Thus the boundary form vanishes at both ends.

This gives the regular self-adjoint differential operator realization: its singular-end condition is ${\lim }_{x\downarrow 1}x(x-1)v^{\prime }=0$, together with $v(R)=0$. The bounded classical eigenfunctions satisfy exactly this no-log choice. More precisely the operator domain is the maximal differential-operator domain in the weighted space with these separated endpoint conditions; the preceding calculation verifies its self-adjoint boundary form, rather than merely the differential expression's formal symmetry.

For the weighted inner product, the factor $w$ cancels the $w^{-1}$ in the Sturm-Liouville operator. Assuming the coefficients are real, two applications of integration by parts giveThis is the boundary form in the Lagrange identity for a Sturm-Liouville operator. The potential terms cancel; if $q$ were complex, an additional integral of $(q-q^{\ast })f^{\ast }g$ would remain, so suitable real coefficients are part of formal self-adjointness. The endpoint restrictions in part (a) eliminate the displayed form for every pair of functions in the domain.

Self-adjoint form. Grouping the two derivative terms gives $[(1+ϵx)u^{\prime }{]}^{\prime }-ϵu=-\lambda u$. Thus a regular Sturm-Liouville operator isThe boundary conditions are mixed Neumann boundary conditions at zero and Dirichlet boundary conditions at $\pi$. For $ϵ\ge 0$, $p>0$, $q\ge 0$ and $w>0$, as required. The boundary term $[p(f{g}^{\prime }-f^{\prime }g){]}_0^{\pi }$ vanishes for two functions satisfying these conditions, confirming self-adjointness on that domain.

Constrained variation. Introduce a Lagrange multiplier $\Lambda$ for the normalization and varyIts first variation isThe endpoint term vanishes because $v^{\prime }(0)=0$ and $\delta v(\pi )=0$. The Euler-Lagrange equation is thereforewhich is exactly the original eigenvalue problem for the boxed choices of $p,q,w$. Multiply by $v$ and integrate by parts. The same endpoint conditions giveThus the stationary functional values are the eigenvalues, with $\Lambda =\lambda$. More generally, before normalization the relevant Rayleigh quotient is

Unperturbed ground mode. At $ϵ=0$, the equation is $-u^{\prime \prime }=\lambda u$. The quadratic form is nonnegative, and zero eigenvalue would force a constant function incompatible with the endpoint at $\pi$. For positive $\lambda$, $u=A\cos (\sqrt{\lambda }x)+B\sin (\sqrt{\lambda }x)$. The derivative condition forces $B=0$, while the other endpoint gives $\cos (\pi \sqrt{\lambda })=0$. HenceThe lowest eigenvalue and normalized eigenfunction are

Trial-function bound. The Rayleigh-Ritz method says that the minimum eigenvalue is no larger than the quotient of any admissible trial function. The same $U_0$ still satisfies both boundary conditions and retains unit norm. Its derivative is $U_0^{\prime }=-\sqrt{2/\pi }\sin (x/2)/2$, givingUsing ${\sin }^2(x/2)=(1-\cos x)/2$,Therefore the affine-parameter Sturm-Liouville eigenvalue bound isThis is an upper bound for every $ϵ\ge 0$, not merely a small-parameter approximation. As a lower-bound check, $p\ge 1$ and the unperturbed variational minimum imply ${\mathcal{R}}_{ϵ}[v]\ge 1/4+ϵ$ for every admissible normalized $v$. Thus

Write $Y={\sum }_{n\ge 0}{b}_n{y}_n$ in the complete Sturm-Liouville eigenfunction expansion. To determine $b_m$, multiply the differential equation by $y_m$ and integrate. The Sturm-Liouville operator has the symmetry identitybecause both functions obey the homogeneous Dirichlet boundary conditions. Consequentlyso that $({\lambda }_m-\alpha )b_m=a_m$. The stipulated exclusion of every eigenvalue makes all denominators nonzero, and givesThis is the resolvent operator applied to $f$ in the weighted space. The same projection shows uniqueness: a homogeneous solution would have $({\lambda }_m-\alpha )b_m=0$ for every $m$, hence would vanish by completeness. For the regular problem the formula gives the solution in the operator domain when $f$ is square-integrable.

Take real admissible functions with finite $F,G$, positive $G$, and a self-adjoint domain for the Sturm-Liouville operator. Compactly supported variations $\eta$ avoid any boundary contribution at infinity. Formal self-adjointness givesFor $\Lambda =F/G$, the quotient rule therefore givesVanishing for arbitrary variations requiresConversely an eigenfunction in the self-adjoint domain makes this variation zero. Multiplying its equation by $y$ and integrating also shows directly that its stationary quotient is $\boxed{\Lambda [y_n]={\lambda }_n}$.

Thus the stationary values are the admissible eigenvalues. Under the usual positive-weight, semibounded hypotheses the lowest is the global minimum, while higher eigenfunctions generally give saddle points. The assumption $y\to 0$ alone is not a complete domain specification for arbitrary $\rho ,\sigma$: the required integrals and boundary form must exist. The rapidly decaying oscillator trials in part (b) meet these conditions.

A weighted Sturm-Liouville operator has the formwhere $p,q,w$ are real, $p>0$ and $w>0$ in the interval's interior, and the coefficients have the regularity required by the derivatives. The corresponding eigenvalue equation is $-(p{y}^{\prime }{)}^{\prime }+qy=\lambda wy$. Its natural weighted inner product is $⟨f,g{⟩}_w=\int wf\bar{g}dx$. Suitable endpoint conditions must also annihilate $[p(f{\bar{g}}^{\prime }-f^{\prime }\bar{g})]$ for a self-adjoint differential operator; the displayed differential expression alone does not specify those conditions.

Two integrations by parts give the boundary form of the Sturm-Liouville operator:The interior terms cancel because $p,q$ are real. Thus the required endpoint condition isFor a regular interval, a standard self-adjoint domain is specified by one real homogeneous Robin boundary condition at each endpoint:where neither pair of real constants is identically zero. At each end, these conditions force the two boundary data vectors of $f,g$ to be proportional, making the displayed form zero. Dirichlet boundary conditions and Neumann boundary conditions are the cases $B=0$ and $A=0$. Coupled conditions also work, for example $y(a)=y(b)$ and $p(a)y^{\prime }(a)=p(b)y^{\prime }(b)$, whose two endpoint contributions cancel. Self-adjointness uses the maximal common domain with these conditions, not an arbitrary smaller symmetric restriction.

The weighted inner product must be nondegenerate. The usual interpretation is $w>0$ almost everywhere, allowing isolated endpoint zeros as in part (b), with the operator defined on its appropriate weighted domain. Merely allowing a weight that is zero on an entire interval would not by itself define the stated operator or a positive-definite inner product on ordinary functions.

For the two sine modes, orthogonality givesThe diagonal potential integrals vanish, while the mixed one isThus ${\int }_0^1\cos (\pi x){\tilde{\psi }}^{2}dx=\alpha /2$, giving the two-sine variational bound for a Dirichlet Mathieu equation:Its derivative vanishes when $1+6{\pi }^2\alpha -{\alpha }^2=0$. Hence the two stationary parameters are $3{\pi }^2\pm \sqrt{9{\pi }^4+1}$. To identify the minimum, use the normalized sine basis: the restricted self-adjoint operator has real symmetric matrixIts smaller eigenvalue is the minimum of the Rayleigh quotient, attained at the negative stationary parameter. ThereforeThe estimate improves on the single-sine value ${\pi }^2$. The Rayleigh-Ritz variational principle remains valid here even though $q(x)=\cos (\pi x)$ changes sign: on this finite interval the Sturm-Liouville operator is self-adjoint and bounded below, which is sufficient for the upper bound.

The weight cancels the denominator in the Sturm-Liouville operator. Two integrations by parts, with real $p,q$, give the Lagrange identity for a Sturm-Liouville operatorThus the boundary form must vanish for every pair in the operator domain:Separated real Robin boundary conditions $A_{s}y(s)+B_{s}p(s)y^{\prime }(s)=0$ at $s=a,b$, with $(A_s,B_s)\ne (0,0)$, make each endpoint contribution vanish. This includes Dirichlet boundary conditions and Neumann boundary conditions. Coupled conditions can also work: periodic $y(a)=y(b)$ and $p(a)y^{\prime }(a)=p(b)y^{\prime }(b)$ cancel the two endpoint contributions. For the regular problem, imposing a complete pair of such self-adjoint boundary conditions gives a self-adjoint differential operator; the displayed identity is the boundary-form criterion.

Factor out $x$ from the differential expression:Comparing with $\mathcal{L}y=w^{-1}[-(p{y}^{\prime }{)}^{\prime }+qy]$ givesMultiplying all three by the same positive constant would give the same Sturm-Liouville operator.

Multiplication by $x^2$ gives the Sturm-Liouville formso the weight is $w(x)=x^2$. Use the weighted inner product $⟨f,g{⟩}_w={\int }_0^{\pi /2}{x}^2\overline{f}gdx$ and the Sturm-Liouville operator $L=-x^{-2}(x^{2}d/dx{)}^{\prime }$. Integration by parts gives its boundary form:The upper endpoint term vanishes by the Dirichlet boundary conditions, and the lower one vanishes for the regular solutions. More precisely, the radial Dirichlet spectrum in three dimensions identifies $L$ under the unitary transformation $v=xu$ with $-v^{\prime \prime }$ on $(0,\pi /2)$ with zero endpoint values; this specifies the self-adjoint realization at the singular endpoint as well as proving the formal identity.

For eigenfunctions $u_a,u_b$ with different eigenvalues, the same identity gives $({\lambda }_b-{\lambda }_a)⟨u_a,u_b{⟩}_w=0$. Taking real eigenfunctions proves the required orthogonality of Sturm-Liouville eigenfunctions.

Write the equation as $(L+\beta )y=-\cos x/x$ with the positive Sturm-Liouville operator from part (a). A Sturm-Liouville eigenfunction expansion therefore gives $a_n=-f_n/(4n^2+\beta )$, where weighted orthogonal projection givesHere the product-to-sum formulas reduce the integral to half the sum of the integrals of $\sin ((2n+1)x)$ and $\sin ((2n-1)x)$. ConsequentlyAt $x=0$, use the continuous limits of the terms. The coefficients are $O(n^{-3})$, so the resulting endpoint series converges. Completeness is the ordinary Fourier sine series completeness after $v=xy$. For an independent check, solving $v^{\prime \prime }-\beta v=\cos x$ with zero endpoint values givesIts numerator vanishes at zero, and it has the same orthogonal projection coefficients. Positivity of $L+\beta$ proves uniqueness.

These Robin boundary conditions give the spatial Sturm-Liouville operator $-d^2/d{x}^2$ one negative eigenvalue, $-1$, with eigenfunction $e^x$. The other eigenvalues are $n^2{\pi }^2$, with eigenfunctions $V_n=n\pi \cos (n\pi x)+\sin (n\pi x)$, $n\ge 1$; there is no zero eigenvalue. Hence the diffusion equation has a growing contribution $c_0{e}^{t+x}$ in addition to decaying modes.

There is no contradiction with ordinary diffusive smoothing. The outward normal derivative is $-u_x$ at the left endpoint and $u_x$ at the right endpoint, so imposing the same coordinate derivative at both ends does not impose dissipative outward-normal data at both ends. The energy identity is $\frac{1}{2}{\partial }_t\int u^2=[u{u}_x{]}_0^1-\int u_x^2$, whose boundary term is $u(1{)}^2-u(0{)}^2$ here and can supply energy. Orthogonal projection of initial data $f$ gives $c_0=2{\int }_0^{1}f{e}^{x}dx/(e^2-1)$; it vanishes precisely when the growing mode is absent.