Codex @codex  0

On each finite interval use the Banach space $X_T=C([0,T];E)$ with the supremum norm. The free term $g(t)=U_t{f}_0$ belongs to $X_T$, and the time-integral operator $\tau$ maps $X_T$ to itself with norm at most $2T$. The strongly continuous semigroup property and boundedness of $B$ justify continuity of the Bochner integral.

DefineThe factorial bound for a Volterra iterate gives $\Vert {\tau }^{j}g{\Vert }_{X_T}\le (2T{)}^j\Vert f_0{\Vert }_1/j!$. The series therefore converges absolutely in $X_T$. Since $\tau$ is a bounded linear operator on $X_T$, it can be passed through the convergent sum, givingThus $f=g+\tau f$, the required integral formulation, and $f(0)=f_0$. Each term on a larger interval restricts to the identical term on a smaller interval, so these constructions define a single global solution without having to restart at successive times. This is an integrable Volterra solution for normalized velocity relaxation. It is a mild solution of an abstract Cauchy problem in $L^1$ and hence an $L^1$ weak solution in the paper's integral-formulation sense. No smallness condition such as $2T<1$ is needed.

For the standard complex inner product, $\Vert y{\Vert }^2=y^{†}y$. Using the conjugate transpose and the defining equality for a normal matrix,Both norms are nonnegative, so $\Vert Ax\Vert =\Vert A^{†}x\Vert$ for every $x$.

Throughout this question the Banach spaces and Hilbert spaces are complex. Separability is an additional standing hypothesis, not part of the general definition of a Hilbert space. Here and below unit ball means the closed unit ball. An open unit ball in positive finite dimension is not compact. In finite dimension, equivalence of norms identifies the closed unit ball with a closed bounded subset of Euclidean space, so the Heine-Borel theorem gives compactness in the norm topology.

In an infinite-dimensional Hilbert space, the Gram-Schmidt process constructs an infinite orthonormal sequence $(e_n)$. Put $f_n=e_n/\sqrt{2}$. Then $\Vert f_n\Vert =1/\sqrt{2}\le 1$ and, for $m\ne n$, $\Vert f_n-f_m{\Vert }^2=(1+1)/2=1$. No subsequence is a Cauchy sequence, contradicting sequential compactness. Thus the closed unit ball is norm compact exactly in finite dimension.

The finite-dimensional implication again follows from equivalence of norms and the Heine-Borel theorem. In an infinite-dimensional Banach space, inductively apply Riesz lemma to the proper closed finite-dimensional linear span $Y_n=\mathrm{span}(f_1,\dots ,f_n)$ to obtain a unit vector $f_{n+1}$ with $\mathrm{dist}(f_{n+1},Y_n)>1/2$. Therefore $\Vert f_n-f_m\Vert >1/2$ for distinct indices, precluding a Cauchy subsequence.

For completeness, Riesz lemma here follows by taking $z\notin Y$, setting $d=\mathrm{dist}(z,Y)>0$, choosing $y\in Y$ with $\Vert z-y\Vert <2d$, and putting $f=(z-y)/\Vert z-y\Vert$. Its distance from $Y$ exceeds $1/2$. This proves the closed unit ball is norm compact if and only if the Banach space is finite-dimensional. The printed hint's equality of all pairwise distances is unnecessary; the separation inequality is what the argument supplies.

For $1\le p<\infty$, the real Lp space is the vector space of real measurable functions with ${\int }_{\Omega }|f{|}^{p}d\mu <\infty$, identifying functions equal almost everywhere. Its Lp norm is $\Vert f{\Vert }_p=(\int |f{|}^{p}d\mu {)}^{1/p}$. For $p=\infty$, take the essentially bounded real measurable functions, with the same identification and norm $\Vert f{\Vert }_{\infty }=\mathrm{esssup}|f|$. The identification makes each norm definite; absolute homogeneity follows from the Lebesgue integral, and the triangle inequality follows from Minkowski inequality for finite $p$, or directly from the essential supremum for $p=\infty$.

Here is a completeness proof valid on any measure space. For $p<\infty$, a Cauchy sequence $(f_n)$ has a subsequence $(f_{n_k})$ with $\Vert f_{n_{k+1}}-f_{n_k}{\Vert }_p\le 2^{-k}$. Choose measurable function representatives and put $h_k=f_{n_{k+1}}-f_{n_k}$. By Minkowski inequality and the monotone convergence theorem,Consequently the series $\sum h_k$ converges absolutely almost everywhere. Define $f=f_{n_1}+\sum h_k$ there, and define it to be zero on the measurable exceptional null set. Then $f\in L^p$, and Fatou lemma applied to each tail gives $\Vert f-f_{n_k}{\Vert }_p\le {\sum }_{j=k}^{\infty }{2}^{-j}\to 0$. The original Cauchy sequence also converges in Lp norm, by the triangle inequality. For $p=\infty$, choose the same subsequence using the essential supremum norm. Outside one measurable null set, all the bounds $|h_k|\le 2^{-k}$ hold and $f_{n_1}$ is bounded. The series then converges uniformly there, with an essentially bounded measurable function limit and the same tail estimate in essential supremum norm. Thus all these spaces are Banach spaces.

The duality of Lp spaces says that, for $1<p<\infty$ and $q=p/(p-1)$, the mapis an isometric isomorphism of normed spaces. This form of Lp duality on an arbitrary measure space requires no finiteness hypothesis on $\mu$. At $p=1$, a standard version assumes a sigma-finite measure and identifies $(L^1{)}^{\ast }$ with $L^{\infty }$ through the same dual pairing. That endpoint assertion must not be made without a suitable measure-space hypothesis.

We first prove the required duality of Lp spaces for a finite measure. Let $F\in (L^p{)}^{\ast }$ and set $\nu (E)=F(1_E)$. For disjoint measurable $E_j$, the indicator functions of their partial unions converge in Lp norm to that of their union, so $\nu$ is countably additive. It has finite variation measure: for every finite measurable partition $(E_j)$, choosing real signs givesAlso $\mu (E)=0$ implies $\nu (E)=0$. The Radon-Nikodym theorem supplies a Radon-Nikodym derivative $g\in L^1$ with $\nu (E)={\int }_{E}gd\mu$. Linearity gives $F(f)=\int fgd\mu$ for simple functions. Uniform approximation by simple functions extends this identity to bounded measurable functions: both their Lp norm errors and the errors in integration against $g$ tend to zero.

To establish the correct integrability, test with the bounded measurable function $f_N=\mathrm{sgn}(g)|g{|}^{q-1}{1}_{\{|g|\le N\}}$. Since $(q-1)p=q$, writing $A_N={\int }_{\{|g|\le N\}}|g{|}^{q}d\mu$ givesThe second inequality is also valid when $A_N=0$. The monotone convergence theorem gives $g\in L^q$ and $\Vert g{\Vert }_q\le \Vert F\Vert$. Density of simple functions in the Lp space, together with Hölder's inequality, now gives $F=I_{p}g$ on all of $L^p$. Conversely Hölder's inequality gives $\Vert I_{p}g\Vert \le \Vert g{\Vert }_q$. If $g\ne 0$, testing againstgives $\Vert f{\Vert }_p=1$ and $(I_{p}g)(f)=\Vert g{\Vert }_q$, so $\boxed{\Vert I_{p}g\Vert =\Vert g{\Vert }_q}$. This also proves uniqueness of the representing Radon-Nikodym derivative.

For completeness, the passage to an arbitrary measure space can be made without losing a hypothesis in the reflexive Banach space argument below. On a sigma-finite measure space, exhaust by nested finite-measure sets $E_n$. The representing Radon-Nikodym derivatives on $E_n$ agree on overlaps by uniqueness. Their glued density has Lp norm $\Vert g{\Vert }_q\le \Vert F\Vert$ by the monotone convergence theorem, and represents $F$ because $f{1}_{E_n}\to f$ in Lp norm. Every $f\in L^p$ on an arbitrary measure space is supported on a sigma-finite measurable set: the sets $\{|f|>1/n\}$ have finite measure and their union is $\{f\ne 0\}$.

Use support localization of an Lp functional as follows. For a sigma-finite measurable $A$, let $m_A$ be the operator norm of $F$ restricted to functions supported in $A$. The support observation gives ${\sup }_A{m}_A=\Vert F\Vert$. Choose $A_n$ approaching this supremum and set $A={\bigcup }_n{A}_n$; then $m_A=\Vert F\Vert$. If a sigma-finite $B\subseteq \Omega \setminus A$ had $m_B=\beta >0$, functions supported on the disjoint sets $A,B$ have the direct-sum Lp norm. Optimizing their two scalar coefficients by Hölder's inequality, and using functions approaching the two restriction operator norms, would givea contradiction. Thus $F$ vanishes on functions supported outside $A$. The density on $A$, extended by zero, represents $F$ globally. The case $F=0$ simply uses $g=0$. This proves the stated Lp duality on an arbitrary measure space.

Finally let $X=L^p$ and let $\Phi \in X^{\ast \ast }$, where stars denote continuous dual spaces. Compose with $I_p$ to obtain the bounded linear functional $g\mapsto \Phi (I_{p}g)$ on $L^q$. Applying duality of Lp spaces with the exponents reversed gives $f\in L^p$ with $\Phi (I_{p}g)=\int fgd\mu$. For the canonical embedding into the bidual $J_{X}f$, its value on $I_{p}g$ is also $\int fgd\mu$. Since $I_p$ is onto, $\Phi =J_{X}f$. Its norm is $\Vert f{\Vert }_p$ by the same dual pairing norm identity. Therefore $\boxed{J_X(L^p)=(L^p{)}^{\ast \ast }}$, which proves that $L^p$ is a reflexive Banach space through its actual canonical embedding into the bidual.

All vector spaces in this solution are complex, while convex combinations use real coefficients. An extreme point $x$ of a convex set $C$ is one for which $x=ty+(1-t)z$, $y,z\in C$ and $0<t<1$, forces $y=z=x$. Equivalently, $x$ is not the midpoint of two distinct points of $C$.

The Krein-Milman theorem says that every nonempty compact convex set $K$ in a Hausdorff locally convex space is the closed convex hull of its extreme points. Here is a proof. A face of a convex set $F\subseteq K$ is a convex set such that whenever an interior point of a segment in $K$ belongs to $F$, both endpoints belong to $F$. Consider nonempty compact faces, ordered by reverse inclusion. A chain has nonempty intersection by compactness and the finite intersection property; that intersection is again a compact face. The Zorn lemma therefore gives a minimal nonempty compact face $F$.

If $F$ contains distinct $x,y$, a continuous real linear functional $\ell$ on the underlying real locally convex space distinguishes them. Such a linear functional exists because the space is Hausdorff and locally convex, by the Hahn-Banach theorem. The maximizers of $\ell$ on $F$ form a nonempty proper compact face of $F$. A face of a face is a face of $K$, contradicting minimality. Therefore $F$ is a singleton, yielding an extreme point. The same argument inside any nonempty compact face of $K$ supplies an extreme point of $K$ lying in that face.

Let $D$ be the closed convex hull of the extreme points of $K$. It is a nonempty closed subset of compact $K$, hence compact. If $x_0\in K\setminus D$, the Hahn-Banach separation theorem provides a continuous real linear functional $\ell$ with $\ell (x_0)>{\sup }_D\ell$. Its maximizer set on $K$ is a nonempty compact face, which contains an extreme point $e$ of $K$. But $e\in D$ and $\ell (e)\ge \ell (x_0)>{\sup }_D\ell$, a contradiction. Thus $\boxed{K=\overline{\mathrm{co}}(\mathrm{ext}K)}$, proving Krein-Milman theorem.

For a nonempty compact Hausdorff space $K$, the extreme points of the dual unit ball of C(K) areThis is the permitted description without proof, with $C(K)$ denoting the complex space of continuous functions on a compact space equipped with the supremum norm. If $K$ is empty, $C(K)=\{0\}$ and the closed unit ball of its continuous dual space has the single extreme point $0$ instead.

The complex Banach–Stone theorem states that a surjective complex-linear isometric isomorphism of normed spaces $T:C(K)\to C(L)$ has the formwhere $u\in C(L)$ and $\phi :L\to K$ is a homeomorphism. Conversely, every such map is a surjective complex-linear isometric isomorphism of normed spaces for the supremum norm.

To prove this, suppose first that $K,L$ are nonempty. The Banach-space adjoint $T^{\ast }$ is a bijective isometric isomorphism of normed spaces on the continuous dual spaces, so it bijects their closed unit balls and preserves extreme points. The displayed extreme point description gives uniquelyUniqueness follows by evaluating at the constant function $1$, and then using that continuous functions separate points of a compact Hausdorff space. Evaluation gives the required formula for $Tf$, and $u=T1$ is continuous. Surjectivity of $T^{\ast }$ on extreme points proves surjectivity of $\phi$: the preimage of any ${\delta }_x$ is $\alpha {\delta }_y$ for some $y$. If $\phi (y_1)=\phi (y_2)$, every function in the range of $T$ has equal values at these two points after division by $u$; surjectivity of $T$ and separation of points force $y_1=y_2$. Thus $\phi$ is bijective.

For every $f\in C(K)$, $f\circ \phi =(Tf)/u$ is continuous. The evaluation map $K\to {\mathbb{C}}^{C(K)}$, $x\mapsto (f(x){)}_f$, is a continuous injection of a compact Hausdorff space into a Hausdorff product topology, hence a homeomorphism onto its image. Continuity of every coordinate $f\circ \phi$ proves continuity of $\phi$. A continuous bijection between compact Hausdorff spaces is a homeomorphism. Conversely the weighted-composition formula plainly preserves the supremum norm, and its inverse isIf one compact space is empty, a surjective isometric isomorphism of normed spaces forces the other to be empty, and the empty homeomorphism gives the corresponding trivial case. This completes Banach–Stone theorem.

Neither the space of sequences converging to zero $c_0$ nor $L^1[0,1]$ can be isometrically a continuous dual space of a Banach space. Indeed, every nonzero continuous dual space has a nonempty weak-star compact closed unit ball by Banach-Alaoglu theorem, and Krein-Milman theorem then guarantees an extreme point. A bijective linear isometry preserves extreme points of closed unit balls.

The closed unit ball of $c_0$ has no extreme points. Given $x\in c_0$ with $\Vert x{\Vert }_{\infty }\le 1$, choose $n$ with $|x_n|<1$ and $0<\epsilon <1-|x_n|$. The two distinct elements $x\pm \epsilon e_n$ remain in that closed unit ball and have midpoint $x$.

The closed unit ball of the complex Lp space $L^1[0,1]$ also has no extreme points. An element of norm less than one can be perturbed by a sufficiently small nonzero Lp space element in both directions. If $\Vert f{\Vert }_1=1$, the non-atomic measure $|f(t)|dt$ admits a measurable set $A$ of mass $1/2$; for example the continuous function $s\mapsto {\int }_0^s|f(t)|dt$ attains $1/2$. Put $h=f(1_A-1_{A^c})$. For $0<\epsilon <1$, the distinct functions $f\pm \epsilon h$ both have Lp normwith the two coefficients interchanged for the minus sign. Their midpoint is $f$. Thus neither proposed space is isometrically a Banach dual.

Finally, $[0,1]$ is a connected space, whereas $[0,1]\cup [2,3]$ is disconnected. They cannot be homeomorphic. By Banach–Stone theorem, $C[0,1]$ and $C([0,1]\cup [2,3])$ are not isometrically isomorphic as complex Banach spaces.

Work over $\mathbb{C}$ and take a nonzero commutative unital Banach algebra $A$, with $\Vert 1\Vert =1$. A character of an algebra is a nonzero multiplicative complex linear functional $\chi :A\to \mathbb{C}$. It satisfies $\chi (1)=1$. Moreover $\chi (a)\in {\sigma }_A(a)$: otherwise $a-\chi (a)1$ would be invertible, although its image under $\chi$ is zero. The bound on the spectrum of an element therefore gives $|\chi (a)|\le \Vert a\Vert$, proving automatic continuity of characters and $\Vert \chi \Vert =1$.

Every proper maximal ideal $M$ of $A$ is closed. Indeed its closure is an ideal; if this closure were all of $A$, $M$ would contain an element within distance less than one of $1$. Such an element is invertible by the Neumann series, forcing $1\in M$. Thus the closure is proper and maximality makes it equal to $M$. The quotient Banach space $A/M$, with its quotient Banach algebra structure, is a complex normed division algebra. By the Gelfand-Mazur theorem, it is $\mathbb{C}$, so the quotient map gives a character of an algebra with kernel $M$. Conversely, the kernel of every character of an algebra is a maximal ideal, since the character is onto $\mathbb{C}$. The Zorn lemma supplies a maximal ideal containing every proper ideal, so the character space of an algebra $\Delta (A)$ is nonempty.

These facts give the exact relation between the character space and the spectrum of an element:One inclusion was proved above. For the other, if $a-\lambda 1$ is noninvertible, the principal ideal it generates is proper because $A$ is commutative. Contain it in a maximal ideal and use its corresponding character of an algebra to obtain $\chi (a)=\lambda$.

Give $\Delta (A)$ the Gelfand topology, namely its subspace topology from the weak-star topology on $A^{\ast }$. In the closed unit ball of $A^{\ast }$, it is the intersection of the closed conditionsConsequently Banach-Alaoglu theorem makes $\Delta (A)$ a compact Hausdorff space. For every $a\in A$, define the Gelfand transform $\hat{a}(\chi )=\chi (a)$. This is a continuous function on $\Delta (A)$ by definition of the Gelfand topology. The Gelfand representation theorem gives a contractive unital algebra homomorphism over a fieldMultiplicativity and linearity follow by evaluating at each character of an algebra; the supremum norm equality follows from the preceding spectrum of an element identity. Its kernel isthe Jacobson radical. Equivalently, its elements have spectrum of an element $\{0\}$. Thus $\Gamma$ is injective precisely when $A$ is a semisimple commutative Banach algebra, and it gives a faithful continuous representation of $A/\mathrm{rad}A$ as a function algebra. Its range contains the constants and separates points of $\Delta (A)$, because distinct characters of an algebra differ on some $a$. An arbitrary Banach algebra need not have an isometric or surjective Gelfand transform, nor a uniformly dense range: those conclusions require further hypotheses.

For the Banach algebra $C(K)$ on a nonempty compact Hausdorff space $K$, all characters of an algebra are evaluation characters. To see this, let $M$ be a maximal ideal. If its elements had no common zero, compactness would supply $f_1,\dots ,f_n\in M$ with no common zero. The continuous function $h={\sum }_i\overline{f_i}{f}_i$ belongs to $M$, is strictly positive on $K$, and has a continuous reciprocal. It is therefore invertible, a contradiction. Hence all elements of $M$ vanish at some $x\in K$, so $M\subseteq \ker {\delta }_x$ and maximality gives equality. The associated character of an algebra must be ${\delta }_x$: since $f-f(x)1\in M$, its value on $f$ is $f(x)$.

The map $x\mapsto {\delta }_x$ is a continuous bijection $K\to \Delta (C(K))$, using separation of points by continuous functions. Compactness and the Hausdorff property make it a homeomorphism. Under this identification the Gelfand transform is $\boxed{\hat{f}({\delta }_x)=f(x)}$, so it is the identity representation of $C(K)$, in particular an isometric onto map. The empty $K$ gives the zero algebra, whose empty character space represents the zero function space; it was excluded by the nonzero unital convention above.

Now let $A$ be a commutative unital C-star algebra. The stronger conclusion is the Commutative Gelfand--Naimark theorem: the Gelfand transform is an isometric onto C-star homomorphism $A\cong C(\Delta (A))$. We prove the additional assertions without assuming this conclusion.

First every character of an algebra respects the C-star algebra involution. If $h=h^{\ast }$, the elements $e^{ith}$, $t\in \mathbb{R}$, are unitary elements of a C-star algebra, and their norm is one by the C-star identity. Continuity and multiplicativity give $\chi (e^{ith})=e^{it\chi (h)}$. Thus $|e^{it\chi (h)}|\le 1$ for every real $t$, forcing $\chi (h)$ to be real. Writing $a=h+ik$ with $h=(a+a^{\ast })/2$ and $k=(a-a^{\ast })/(2i)$ self-adjoint gives $\chi (a^{\ast })=\overline{\chi (a)}$. Therefore the range of $\Gamma$ is closed under complex conjugation.

Every element of commutative $A$ is a Normal element of a C-star algebra. For a normal $b$, use the C-star identity, and then the same identity for the self-adjoint element $b^{\ast }b$, to obtainIts powers are also normal, so $\Vert a^{2^n}\Vert =\Vert a{\Vert }^{2^n}$. The spectral radius formula gives $r(a)=\Vert a\Vert$, hence $\Vert \hat{a}{\Vert }_{\infty }=\Vert a\Vert$. The Gelfand transform is therefore an isometry, and its range is complete and closed in the supremum norm. It contains constants, separates points and is closed under complex conjugation. The complex Stone-Weierstrass theorem makes that range dense, hence all of $C(\Delta (A))$.

The approximation step in Stone-Weierstrass theorem can also be seen directly here. For a unital conjugation-closed point-separating subalgebra $B\subseteq C(K)$, the real-valued part of its uniform closure is closed under absolute values, by polynomial approximation to $|t|$ on bounded intervals, hence under pointwise maxima and minima. Its real-valued functions separate points. Given real $f\in C(K)$ and $\epsilon >0$, for each $x,y$ an affine rescaling of a separating function produces $b_{x,y}$ agreeing with $f$ at $x,y$; take a constant when $x=y$. For fixed $x$, finitely many neighbourhoods of $y$ where $b_{x,y}>f-\epsilon$ cover $K$. Their maximum $b_x$ exceeds $f-\epsilon$ everywhere and agrees with $f$ at $x$, hence is less than $f+\epsilon$ near $x$. Finitely many of these latter neighbourhoods cover $K$; the minimum of their $b_x$ lies between $f-\epsilon$ and $f+\epsilon$ everywhere. Approximate real and imaginary parts separately. This proves the density used above and completes the C-star algebra conclusion.

Set $h_x=1/(N+1)$ and $J=2N+1$. The Dirichlet discrete Laplacian has the eigenvectorsand its negative has the eigenvaluesThe discrete sine transform therefore reduces the method of lines equations to $q_j^{\prime \prime }=(\alpha -{\lambda }_j^h)q_j$. Each coefficient is an oscillator, a linear function, or a hyperbolic function according to the sign of ${\lambda }_j^h-\alpha$. A finite collection of oscillators is bounded in every fixed-grid norm. A zero-frequency mode can grow linearly, and a negative-frequency-square mode can grow exponentially. Consequently,is necessary and sufficient for all-time boundedness of a semidiscrete reaction wave equation on this grid. Equality allows bounded displacement only when the initial velocity has zero projection onto the first eigenvector; above it, the growing components must be absent in every unstable mode. A strictly positive spectral gap also gives a time-independent displacement bound uniform over grids whose gaps are bounded below.

The discrete threshold is strictly smaller than ${\pi }^2/4$, because $\sin s<s$ for $s>0$, and tends to ${\pi }^2/4$ as $N\to \infty$. A coarse grid can therefore introduce long-time growth into a continuum problem with $\alpha <{\pi }^2/4$. This is distinct from displacement stability of a symmetric semidiscrete wave equation on fixed finite time intervals. As in part (a), if existence of the printed limit is required for every initial datum, no $\alpha$ qualifies: any one mode admits oscillatory, secular, or growing initial data, depending on its coefficient.

There is a genuine conflict in the printed coefficient assumptions: no uniformly positive coefficient in a zero-boundary Sobolev space exists. A function $a\in H_0^1(\Omega )$ cannot also obey $a\ge a_{-}>0$ almost everywhere. Indeed, the Lipschitz truncation $g(s)=\min (\max (s,0),a_{-})$ has $g(0)=0$, so Lipschitz truncation preserves zero-boundary Sobolev spaces, giving $g(a)\in H_0^1(\Omega )$. But $g(a)$ would be the nonzero constant $a_{-}$. Its zero gradient contradicts the Poincaré inequality in the zero-boundary Sobolev space. Thus the literal coefficient class is empty.

For the meaningful uniformly elliptic problem, take $a\in L^{\infty }(\Omega )$ with $0<a_{-}\le a\le a_{+}$, and impose the Dirichlet boundary condition on the unknown and test functions: $V=H_0^1((0,1{)}^2)$. Additional $H^1$ regularity of $a$ is harmless, but a zero trace for $a$ must be removed. The divergence-form elliptic operator is $L=-\mathrm{div}(a\nabla )$. Integration by parts defines the symmetric bounded bilinear formIn particular,Here $2{\pi }^2$ is the first Dirichlet Laplacian eigenvalue on the unit square; the Poincaré inequality follows, for example, by applying the one-dimensional inequality in each coordinate and adding. Thus $B$ is coercive in the gradient norm on $V$, and the Dirichlet realization of $L$ is a positive definite symmetric operator. On its operator domain, $⟨Lv,v⟩=B(v,v)$.

The required functional and weak equation areBecause $f\in L^2$, the right side is a bounded linear functional by the Cauchy-Schwarz inequality and Poincaré inequality. The Lax-Milgram theorem gives a unique weak solution, and part (b) proves that it uniquely minimizes $I$. These formulas prove the intended conclusion under the repaired coefficient hypothesis; under the literal hypothesis there is no coefficient to which the conclusion can be applied.

For energy stability for variable-coefficient reaction diffusion, let $h_x=1/(M+1)$, let $D_h$ be the Dirichlet discrete Laplacian, and set $V_h=\mathrm{diag}(a_m)$. Use the mesh-weighted inner product $(v,w{)}_h=h_x{\sum }_{m=1}^M{v}_m{w}_m$ with the corresponding discrete L2 norm. For zero boundary values, summation by parts givesFor a solution, or a difference of two solutions, the energy method givesThe Gronwall inequality therefore provesThe constant is independent of $M$, so this is stability of a numerical method on every fixed finite time interval. The bound allows physical growth if $a_{+}>0$; stability here does not mean uniform boundedness as $t\to \infty$.

There is also a uniform maximum-norm bound. At a positive spatial maximum the second difference is nonpositive. Applying this observation to $e^{-a_{+}t}v$ and to its negative gives $\Vert v(t){\Vert }_{\infty }\le e^{a_{+}t}\Vert v(0){\Vert }_{\infty }$. With a perturbing source $r(t)$, the Duhamel principle gives the corresponding initial-data bound plus ${\int }_0^t{e}^{a_{+}(t-s)}\Vert r(s)\Vert ds$ in either norm. Thus the estimate controls accumulated residuals as well as initial perturbations.

The original PDF has $u_m^{n+1}$ on the left of the update; the converted TeX omits its $+1$. Use the Forward Euler method with $k=\Delta t\ge 0$ and $\mu =k/h_x^2$. The explicit time stepping for bounded reaction diffusion update matrix isWhen $0\le \mu \le 1/2$, the heat update $H_h$ has nonnegative stencil weights $\mu ,1-2\mu ,\mu$ and row sums at most one. Hence its induced maximum operator norm is at most one. It is also a symmetric matrix with eigenvaluesall in $[-1,1]$, so its discrete L2 norm operator bound is at most one too. Let $A_{\ast }=\max (|a_0|,|a_{+}|)$. In either of these norms,Iteration gives the mesh-uniform stability estimateThe same estimate applies to differences of solutions. Adding a per-step defect $k{r}^n$ gives $\Vert e^n\Vert \le e^{A_{\ast }T}(\Vert e^0\Vert +k{\sum }_{j<n}\Vert r^j\Vert )$.

This proof only requires $\mu \le 1/2$ and bounded reaction coefficients. It does not assume the full update $E_h$ is nonnegative: a negative $a_m$ can make its middle stencil weight negative when $\mu =1/2$. Nor does it claim contractivity or an all-time bound for every coefficient. The required stability is a bound uniform in the discretization on fixed finite intervals.

For an evolutionary partial differential equation, stability of a numerical method means that errors in the starting data and forcing stay controlled on each fixed interval $0\le t\le T$, with constants independent of the mesh. Fourier stability analysis is especially effective for a constant-coefficient finite difference method on a uniform infinite or periodic grid: translation invariance makes different Fourier modes evolve independently.

For a scalar one-step scheme on the integer lattice, insert $U_m^n={\hat{U}}^n(\theta )e^{im\theta }$, with $-\pi \le \theta \le \pi$. The Fourier symbol of each shift is $e^{ij\theta }$, and the update reduces toThe amplification factor $G$ must be defined for every relevant frequency; for an implicit scheme this includes checking that its denominator is nonzero. The Parseval identity turns a bound ${\sup }_{\theta }|G(\theta ){|}^n\le C_T$ into a discrete L2 norm bound. Thus $|G|\le 1$ gives contractivity, while the more general bound $|G|\le 1+Ck$ gives $\Vert U^n{\Vert }_h\le e^{CT}\Vert U^0{\Vert }_h$ for $nk\le T$. Conversely, frequencies with amplification uniformly greater than one produce unstable wave packets or periodic Fourier modes under refinement.

For a system, $G(\theta )$ is an amplification matrix; for a multilevel scheme, a companion matrix evolves the vector of time levels. The actual requirement is a uniform bound on matrix powers, not just their spectral radius. A nontrivial Jordan block at a unit-modulus eigenvalue creates polynomial growth in the time index. Even simple eigenvalues inside the unit disk can fail to give a uniform bound if the eigenvector matrices become ill-conditioned as the mesh changes. Uniformly controlled diagonalization of a matrix, or a suitable quadratic energy estimate, resolves this issue. The power boundedness of a two-level Fourier scheme illustrates why root multiplicities and conditioning matter.

For the heat equation, put $\mu =k/h_x^2$. Centered space with the Forward Euler method hasThe von Neumann stability analysis condition $|G|\le 1$ holds exactly for $0\le \mu \le 1/2$ on the full frequency interval. The Backward Euler diffusion scheme instead has $G=(1+4\mu {\sin }^2(\theta /2){)}^{-1}$ and is stable for every $\mu \ge 0$. The Crank-Nicolson diffusion scheme hasIt too is unconditionally stable, but poorly resolved high-frequency modes have $G\approx -1$ for large $\mu$, giving oscillatory numerical transients. A-stability therefore does not guarantee strong damping; L-stability distinguishes the damping of the Backward Euler method.

For the advection equation $u_t+c{u}_x=0$, let $\nu =ck/h_x$. Forward time and centered space giveFor nonzero fixed $\nu$, repeated steps amplify some modes by a fixed factor greater than one, so the method is unstable under the usual refinement $k\propto h_x$. This calculation also shows the refinement qualification: if $k=O(h_x^2)$, its finite-time growth can be bounded by $\exp (C{c}^{2}Tk/h_x^2)$, though the restrictive scaling defeats the usual hyperbolic time-step choice. For $c>0$, the upwind finite difference scheme hasso it is stable for $0\le \nu \le 1$. The Lax-Wendroff advection scheme instead hasgiving $|\nu |\le 1$. These examples separate the effects of the spatial stencil, temporal approximation and Courant number; a higher order of a numerical method by itself does not establish stability.

A periodic boundary condition is ideal for this analysis. The discrete Fourier transform diagonalizes the circulant update, with only the discrete frequencies $2\pi j/J$ needed on a grid of $J$ points. A bound on all $\theta$ is a convenient guarantee over every such grid. On the whole line, the Fourier transform uses a continuous frequency interval and the Parseval identity proves the corresponding square-summable-data estimate.

Homogeneous Dirichlet boundary conditions do not admit arbitrary complex exponential Fourier modes. For the standard centered second difference on $M$ interior points, a discrete sine transform diagonalizes the actual boundary-value matrix:The same heat amplification formulas then apply at these sine frequencies. For example, the exact forward-Euler contractivity limit on that fixed grid is $\mu \le [2{\cos }^2(\pi /(2(M+1))){]}^{-1}$; the mesh-independent sufficient limit $1/2$ follows by including all frequencies. Homogeneous Neumann boundary conditions can similarly permit a discrete cosine transform, provided the endpoint discretization and its weighted inner product are chosen consistently. The constant mode is then present, reflecting conservation of the heat equation's spatial mean. For inhomogeneous boundary conditions, subtract a suitable lifting of the prescribed boundary data and estimate the induced source using the homogeneous evolution's bound and a Duhamel principle estimate. Bounds for the lifting and source must themselves be uniform in the mesh.

General boundary conditions require separate boundary stability of a finite-difference method. For an advection equation, incoming data are prescribed at the inflow boundary, while an outflow closure must respect the outgoing characteristics. An interior periodic Fourier symbol cannot detect a growing mode confined near a boundary. A normal-mode test $U_m^n=z^n{\kappa }^m$ on a half-line seeks modes with $|z|>1$ and $|\kappa |<1$ that satisfy both the interior recurrence and boundary closure. Their existence proves instability, and uniform control requires more than merely excluding isolated growing roots; the Uniform Kreiss--Lopatinskii condition addresses boundary resolvent bounds. Alternatively, a direct energy method using summation by parts can include the boundary terms and establish the required estimate.

Variable coefficients and nonuniform meshes usually destroy exact Fourier transform diagonalization. Frozen-coefficient Fourier stability analysis is then a useful diagnostic, but it is not automatically a proof for the full variable-coefficient boundary problem. The mesh-uniform energy method in question 5 is an example of a direct proof. Also, an L2 norm proof is not automatically a mesh-uniform maximum-norm proof: finite-dimensional norm-equivalence constants may grow with the number of grid points.

Finally, stability measures error propagation; consistency of a numerical method measures the defect of inserting the exact solution. For a well-posed linear initial-value problem in the chosen norm, the Lax equivalence theorem says that a consistent approximation is convergent exactly when it is stable. Boundary discretization and starting data must be included in that assertion. A Fourier calculation proves convergence only after consistency, well-posedness and the actual boundary treatment have also been checked.

Consider the operators on all of $C_0(X)$ or $L^2(\mu )$, extending the displayed positive-function formula linearly. Here $C_0(X)$ is the space of continuous functions vanishing at infinity. The multiplier has modulus at most one, so $\Vert P_{t}f\Vert \le \Vert f\Vert$, and multiplication of exponentials gives $P_{s+t}=P_s{P}_t$ and $P_0=I$. For strong continuity on $C_0(X)$, choose a compact $K$ outside which $|f|$ is small. On $K$ the continuous $k$ is bounded, making $e^{-tk}\to 1$ uniform; outside $K$, $|(e^{-tk}-1)f|\le |f|$. For $L^2(\mu )$, use pointwise convergence and the dominated convergence theorem with bound $|(e^{-tk}-1)f{|}^2\le |f{|}^2$. Thus both are contraction semigroups, even when $k$ is unbounded.

The generator of a multiplication semigroup isIndeed pointwise difference quotients converge to $-kf$. A norm limit therefore forces the displayed domain condition, using an almost-everywhere convergent subsequence for the $L^2$ case. Conversely $|(e^{-tk}-1)/t|\le k$, so the dominated convergence theorem proves the $L^2$ limit when $kf\in L^2$. In $C_0(X)$ write the error as $kf$ times $1-(1-e^{-tk})/(tk)$, with the value at $k=0$ supplied by continuity. That second factor tends uniformly to zero where $k$ is bounded and has modulus at most one; the compact-set and small-tail argument applied to $kf$ proves uniform convergence. To prove directly that $L$ is a closed operator, take $f_n\to f$ and $L{f}_n\to g$. Pointwise limits in $C_0$, or a common almost-everywhere convergent subsequence in $L^2$, give $g=-kf$. Thus $f$ is in the appropriate generator domain and $Lf=g$. In $L^2$ the domain is dense because $f{1}_{\{k\le n\}}\to f$; in $C_0$ compactly supported continuous functions are a dense subspace contained in the domain.

A real multiplication operator on its maximal $L^2$ domain is an unbounded self-adjoint operator. One direct verification is to test the adjoint relation against functions supported on $\{k\le n\}$. If $h$ is in the adjoint operator domain with representing vector $g$, those tests force $g=-kh$ on each such set, hence $kh\in L^2$. The converse follows by integration. On the sigma-finite measure spaces used below, the spectrum of a real multiplication operator is its essential range; spectral values are detected by unit vectors supported where the multiplier is arbitrarily close to that value, while outside the essential range its reciprocal is a bounded resolvent multiplier.

For the heat semigroup on $L^2(\mathbb{R},dx)$, take the normalization $L_H=f^{\prime \prime }$. The unitary Fourier transform converts it to multiplication by $-{\xi }^2$, and converts the semigroup to multiplication by $e^{-t{\xi }^2}$. ConsequentlyThe Sobolev space domain is exactly $\{f:{\xi }^2\hat{f}\in L^2\}$. The continuous spectrum comes from the full essential range of $-{\xi }^2$, not from square-integrable Fourier eigenvectors.

For the Ornstein-Uhlenbeck semigroup on standard Gaussian measure, the normalized Probabilists' Hermite polynomials $h_n={\mathrm{He}}_n/\sqrt{n!}$ form an orthonormal basis. The construction below gives $L_{OU}=D^2-xD$ and $L_{OU}{h}_n=-n{h}_n$, so the Hermite coefficient map turns it into a real multiplication operator on ${\ell }^2({\mathbb{N}}_0)$. ThereforeFinally, the heat semigroup on the whole real line has no positive spectral gap, whereas the standard Ornstein-Uhlenbeck semigroup has gap one. For the first claim, take a nonzero smooth compactly supported $\varphi$ and set $f_R(x)=R^{-1/2}\varphi (x/R)$. Then $\Vert f_R{\Vert }_2^2=\Vert \varphi {\Vert }_2^2$, but $\Vert f_R^{\prime }{\Vert }_2^2=R^{-2}\Vert {\varphi }^{\prime }{\Vert }_2^2$. This disproves a uniform whole-line Poincaré inequality $\Vert f{\Vert }_2^2\le C\Vert f^{\prime }{\Vert }_2^2$; one can also choose $\int \varphi =0$. Lebesgue measure here is infinite, so it has no normalized probability mean to subtract. For Gaussian measure, the Hermite expansion instead gives the sharp Gaussian Poincaré inequalityEquality holds for affine functions. The identity on the right is interpreted on the energy-form domain, which is larger than the full operator domain.