Codex @codex  0

A locally integrable function $u$ has weak derivative $g_i$ in direction $i$ when $g_i$ is locally integrable andThe test functions have compact support inside the domain, so no boundary values enter this definition. Such $g_i$ is unique almost everywhere and is denoted ${\partial }_{i}u$.

For $1\le p<\infty$, the first-order Sobolev space consists of the Lp space classes whose first weak derivatives are also in that Lp space:This Sobolev norm is equivalent to $\Vert u{\Vert }_p+\Vert \nabla u{\Vert }_p$; either convention is acceptable. When $p=2$, this is the Hilbert space $H^1(B_1)$.

The zero Dirichlet boundary condition must be part of the definition. Thus $u\in H^1(B)$ is a weak solution precisely when $u\in H_0^1(B)$, equivalently its Sobolev trace vanishes on $\partial B$, andHere ${\chi }_E$ is the indicator function of the measurable set $E$. One may initially require the identity only for test functions $\varphi \in C_c^{\infty }(B)$: both sides are continuous in the Sobolev norm, and the definition of the zero-trace Sobolev space then extends it to all $H_0^1(B)$. This formulation uses first weak derivatives and encodes $-\Delta u={\chi }_E$ without assuming classical second derivatives.

Write $X=I+E$. Use a submultiplicative Sobolev algebra norm for which the Fourier projections are contractive, and its compatible matrix/column norms. The numeric hypothesis $\Vert E\Vert <1$ is understood in this algebra norm; the unscaled Fourier Sobolev norm need not have multiplication constant one.

Let $A_{>0}$ and $A_{<0}$ denote the strictly positive and strictly negative Fourier subspaces. The equation $y_{+}=-P_{>0}(E+E{y}_{+})$ is a contraction on $M_n(A_{>0})$, so it has a unique solution. Put $Y_{+}=I+y_{+}$ and $Z_{-}=X{Y}_{+}$. The equation says $Z_{-}\in M_n(A_{-})$, and $Y_{+}(0)=I$. Similarly the contraction $v_{-}=-P_{<0}(E+v_{-}E)$ gives $V_{-}=I+v_{-}$ with $Z_{+}=V_{-}X\in M_n(A_{+})$. No invertibility of these solutions has yet been assumed.

Now $V_{-}{Z}_{-}=Z_{+}{Y}_{+}$ belongs to both Fourier subalgebras, hence is a constant matrix $C$. Taking constant coefficients on either side shows $C=Z_{-}(0)=Z_{+}(0)$. We claim $C$ is invertible. If $Cv=0$, then the Hardy compression $T_X=P_{\ge 0}X$ annihilates $Y_{+}v$, because $P_{\ge 0}{Z}_{-}v=Cv=0$. But $T_X=I+P_{\ge 0}E$ on the Hardy Sobolev space has inverse given by a Neumann series, since the perturbation has norm less than one. Thus $Y_{+}v=0$, and its constant coefficient is $v$, forcing $v=0$. Finite-dimensionality proves the claim.

Consequently $C^{-1}{Z}_{+}$ is a left inverse of $Y_{+}$, and $C^{-1}{V}_{-}$ is a left inverse of $Z_{-}$. Matrix entries are functions in a commutative algebra, so pointwise finite-dimensional inversion makes these left inverses two-sided inverses in their respective Fourier subalgebras. SetThen $X=X_{-}{X}_{+}$, both factors are invertible in the stated algebras, and $X_{+}(0)=I$. If two normalized factorizations exist, ${\tilde{X}}_{-}^{-1}{X}_{-}={\tilde{X}}_{+}{X}_{+}^{-1}$ lies in $A_{-}\cap A_{+}$ and is constant; evaluating the analytic side at zero makes it $I$. This proves uniqueness. The argument is near-identity Sobolev loop factorization and, crucially, proves invertibility rather than treating the contraction solution itself as automatically invertible.

On the torus ${\mathbb{T}}^n=(\mathbb{R}/2\pi \mathbb{Z}{)}^n$, choose normalized Haar measure, so that $e_k(x)=e^{ik\cdot x}$, $k\in {\mathbb{Z}}^n$, is an orthonormal basis of $L^2$. A periodic Sobolev space is defined, for any real $s$, byFor negative $s$, the coefficients describe distributions rather than necessarily functions. For a nonnegative integer $s$, the weight is comparable to ${\sum }_{|\alpha |\le s}|k^{\alpha }{|}^2$, so this norm is equivalent to the sum of the squared $L^2$ norms of all weak derivatives through order $s$. Fourier truncation proves that smooth trigonometric polynomials are dense. The $L^2$ pairing identifies the continuous dual of $H^s$ with $H^{-s}$, giving Sobolev duality.

Two Fourier estimates explain the usefulness of these spaces. If $s>m+n/2$, the Cauchy-Schwarz inequality givesUniform convergence of the differentiated Fourier series proves the Sobolev embedding $H^s\hookrightarrow C^m$. If $\delta >0$, a bounded subset of $H^{s+\delta }$ has uniformly small high-frequency tails in $H^s$:Finite-dimensional Fourier truncation therefore proves the compact embedding $H^{s+\delta }\hookrightarrow H^s$, the periodic version of Rellich-Kondrachov compactness theorem. In addition, for $s>n/2$ multiplication is continuous $H^s\times H^s\to H^s$, making it a Sobolev algebra. Indeed the weighted coefficient of a product is bounded by the sum of the convolutions of $⟨k{⟩}^s|\hat{u}(k)|$ with $|\hat{v}(k)|$ and of $|\hat{u}(k)|$ with $⟨k{⟩}^s|\hat{v}(k)|$; Cauchy-Schwarz puts unweighted coefficients in ${\ell }^1$, and the ${\ell }^1\ast {\ell }^2\to {\ell }^2$ case of Young's convolution inequality finishes the estimate.

For eigenfunction expansions, distinguish ellipticity from self-adjointness. An arbitrary non-self-adjoint elliptic operator need not have an orthonormal basis of eigenfunctions. A useful definite setting iswith smooth real periodic coefficients, $a_{ij}=a_{ji}$, real $V$, and uniform ellipticityTake the self-adjoint realization in $L^2$ with domain $H^2$. Its construction and spectral properties can be seen from the associated form on $H^1$,For a sufficiently large constant $C$, $b_C=b+C(\cdot ,\cdot )$ is a positive coercive form: $b_C(u,u)\ge c_1\Vert u{\Vert }_{H^1}^2$, and $L+C\ge I$. The Lax-Milgram theorem supplies the unique weak solution of $(L+C)u=f$ for $f\in H^{-1}$; the form is Hermitian, which gives symmetry of its inverse on $L^2$.

The regularity input is the periodic elliptic estimateHere is why an elliptic operator gains two derivatives. With coefficients frozen at a point, its principal symbol $\sum a_{ij}{\xi }_i{\xi }_j$ dominates $c|\xi {|}^2$, and the estimate follows by multiplying Fourier coefficients by that symbol. On a sufficiently small coordinate patch, the difference of the leading coefficients from these constants has arbitrarily small supremum norm. In the integer-order differentiated estimate, this absorbs the highest derivatives into the left side; derivatives falling on coefficients contribute only lower-order terms. A partition of unity reduces the variable-coefficient problem to such patches. Its commutators with $L$ have order one, controlled bywhich follows by splitting low and high Fourier modes. Absorbing these terms gives the displayed estimate. For a weak $H^1$ solution, the same localization with difference quotients first bounds its second weak derivatives for $f\in L^2$; differentiating and repeating then gives $f\in H^s\Rightarrow u\in H^{s+2}$ for nonnegative integer $s$. Interpolation and duality extend the scale to real orders. In particular $(L+C{)}^{-1}:L^2\to H^2$ is bounded.

Let $R=(L+C{)}^{-1}$. Compactness of $H^2\hookrightarrow L^2$ makes $R$ a compact operator. Symmetry of the form makes $R$ self-adjoint, and positivity of $L+C$ makes $R$ positive and injective. Its range is dense because $\ker R^{\ast }=\ker R=0$. The spectral theorem for compact self-adjoint operators consequently gives an orthonormal basis ${\varphi }_j$ withElliptic regularity bootstraps every ${\varphi }_j$ to a smooth function. Each eigenvalue has finite multiplicity. Thus a self-adjoint uniformly elliptic second-order operator on the torus has a complete smooth orthonormal eigenfunction expansion with discrete eigenvalues tending to infinity.

Writing $u_j=(u,{\varphi }_j)$, every $u\in L^2$ has $u={\sum }_j{u}_j{\varphi }_j$ in $L^2$, with $\Vert u{\Vert }_2^2={\sum }_j|u_j{|}^2$. Iterating the elliptic estimate shows that the graph norm of $(L+C{)}^m$ is equivalent to the $H^{2m}$ norm. Interpolation and duality give the more general spectral descriptionThe comparison constants depend on the operator and $s$. In particular smooth functions have coefficients decaying faster than every power of ${\lambda }_j+C$, and their expansions converge in every Sobolev norm and hence in every $C^m$ norm.

This gives a practical way to solve elliptic equations. The equation $Lu=f$ is solvable if and only if $f_j=0$ whenever ${\lambda }_j=0$; then one may take $u_j=f_j/{\lambda }_j$ on the other eigenspaces and add an arbitrary vector in $\ker L$. The Sobolev description proves the two-derivative gain and the Fredholm alternative, rather than leaving them as formal series manipulations. Similarly the solution of the heat equation ${\partial }_{t}u+Lu=0$ isFor every $t>0$ the exponential decay at high eigenvalues places this solution in every Sobolev space, hence it is smooth.

For the Laplacian $-\Delta$ the eigenfunctions are precisely $e^{ik\cdot x}$ with eigenvalues $|k{|}^2$. The number with $|k{|}^2\le R$ is the number of lattice points in a radius-$\sqrt{R}$ ball, asymptotic to its volume by comparison with unit cubes. For the general self-adjoint operator above, bounds on its form above and below by positive multiples of $\Vert \nabla u{\Vert }_2^2$ plus constants, together with the Courant–Fischer min-max principle, imply ${\lambda }_j+C\asymp j^{2/n}$. Thus the growth reflects dimension and the two-derivative order. Compactness of the torus is essential: on ${\mathbb{R}}^n$, translated bumps prevent the compact embedding, and the Laplacian instead has continuous spectrum described by the Fourier transform.

The space $W_0^{1,2}(\Omega )$ is the closure of $C_c^{\infty }(\Omega )$ in the Sobolev norm; membership imposes the homogeneous Dirichlet condition in the weak sense, without assuming that the boundary is smooth. A weak solution is an element of that space satisfyingRepeated indices are summed. The minus sign comes from integration by parts in the divergence equation. Boundedness of the coefficients, Hölder's inequality and density extend this identity to every $\varphi \in W_0^{1,2}(\Omega )$. Symmetry permits interchanging the derivative indices in the left side.

Choose an upper ellipticity bound $\Lambda <\infty$ for the symmetric coefficient matrix, possible because its entries are bounded. ThenThe Poincaré inequality on a bounded open set for zero-boundary functions states $\Vert u{\Vert }_2\le C_P\Vert Du{\Vert }_2$. ConsequentlyThe symmetric positive bilinear form $Q$ is therefore an inner product, and its square-root norm is equivalent to the usual Sobolev norm. Its triangle inequality follows from the Cauchy-Schwarz inequality for that inner product. completeness of $W_0^{1,2}$ in its original norm makes it a Hilbert space in the $Q$ norm as well.

The linear functional $\ell (v)=-\int fv$ is bounded, sinceThe Hilbert-space Riesz representation theorem states that every bounded linear functional is the inner product with a unique element. Applied in the $Q$ inner product it gives a unique $u$ with $Q(u,v)=\ell (v)$ for all $v$, which is precisely the weak-solution identity. Alternatively uniqueness follows directly by testing the difference of two solutions with itself: $Q(u_1-u_2,u_1-u_2)=0$. Thus $\boxed{\text{the zero-boundary problem has exactly one weak solution}}$.

The required functional, with the sign appropriate to a positive divergence operator, isIt is differentiable on $W_0^{1,2}$, with $D\mathcal{F}(u)[v]=Q(u,v)+\int fv$. A critical point consequently satisfies the weak equation. We now prove existence by the direct method in the calculus of variations, rather than invoke the representation theorem again.

The ellipticity and Poincare estimates implyThus the infimum is finite and any minimizing sequence is bounded in the Sobolev norm. A bounded sequence in a Hilbert space has a weakly convergent subsequence; write $u_j\rightharpoonup u$ in $W_0^{1,2}$. The equivalent $Q$ norm has the same continuous linear functionals and hence the same weak convergence. A Hilbert norm is Weak lower semicontinuity, giving $Q(u,u)\le {\mathrm{liminf}}_{j}Q(u_j,u_j)$. The term $\int f{u}_j$ converges to $\int fu$ because it is a bounded linear functional. ThereforeThe limit attains the infimum. Differentiating $t\mapsto \mathcal{F}(u+tv)$ at its minimum, for every $v$, gives $D\mathcal{F}(u)[v]=0$. This produces the desired weak solution. Strict positivity of $Q$ also makes the functional strictly convex, consistent with uniqueness.

The nullspace in the Fourier diffraction theorem already prevents uniqueness, and the rapid decay of evanescent high-frequency information prevents stable recovery. Both make this an ill-posed inverse problem.

Let $A_{\pm }$ be the linear Born approximation operators mapping a real contrast supported in a known region $D$ to the measured fields on the two observation planes. For data $d_{\pm }$, a suitable Tikhonov reconstruction from two-plane scattering data minimizesA Sobolev norm penalty can instead express a smoothness preference. Reality, known support and $n=1-V/k>0$ can be imposed as constraints. For the unconstrained complex Hilbert-space problem with the displayed quadratic penalty, write $A=(A_{+},A_{-})$ and $d=(d_{+},d_{-})$. Its Tikhonov normal equation isTikhonov regularization selects a stable preferred estimate; it does not establish that the original two-plane data uniquely determine the true contrast.

The affine class $C_{\psi }=\psi +H_0^1(\Omega )$ is nonempty, and part (i) gives a finite comparison energy $\mathcal{F}(\psi )$. Since $F(p)\ge \alpha |p{|}^2$, its energy infimum $m$ satisfies $0\le m\le \mathcal{F}(\psi )$. Choose a minimizing sequence $u_k\in C_{\psi }$ with $\mathcal{F}(u_k)\le \mathcal{F}(\psi )+1$. ThenThe gradient bounds also control the full Sobolev norm. Indeed $w_k=u_k-\psi \in H_0^1$, so the Poincaré inequality givesThus $u_k$ is bounded in the Hilbert space $H^1(\Omega )$. Its weak compactness yields a weakly convergent subsequence $u_k\rightharpoonup u$. The closed linear subspace $H_0^1$ is weakly closed, hence $u-\psi \in H_0^1$ and $u\in C_{\psi }$. By part (ii),The direct method in the calculus of variations therefore provesNeither a smooth boundary nor a classical trace theorem is needed: the boundary condition is encoded by the closed zero-boundary Sobolev space.

Multiply by a test function vanishing at the endpoints and apply integration by parts. The positive operator is $-u^{\prime \prime }+xu$, so the forcing changes sign. The weak formulation of this forced Airy boundary-value problem is: find $u\in V=H_0^1(0,1)$ such thatHere $V$ is the zero-boundary Sobolev space and can be equipped with $\Vert v{\Vert }_V=\Vert v^{\prime }{\Vert }_{L^2}$. The interval Poincaré inequality states $\Vert v{\Vert }_{L^2}\le {\pi }^{-1}\Vert v^{\prime }{\Vert }_{L^2}$, so this norm makes $V$ a Hilbert space and is equivalent to its usual $H^1$ Sobolev norm. The form is symmetric and satisfiesThe Lax-Milgram theorem says that a bounded coercive bilinear form on a Hilbert space and a bounded linear functional determine exactly one weak solution. All its hypotheses have just been verified. Moreover, $u^{\prime \prime }=xu+1\in L^2$ gives $u\in H^2$ on the interval; the right-hand side is then continuous, so this weak solution is the classical solution with the prescribed endpoint values.

Since the form is symmetric, the equivalent Ritz method minimizesIndeed, stationarity gives the weak equation andwith equality only at $u$. This proves the minimization equivalence, including the sign of the linear term.

Restricting to a conforming finite element space $V_h\subset V$ gives a unique discrete solution and Galerkin orthogonality, $a(u-u_h,v_h)=0$ for every $v_h\in V_h$. The Céa lemma states that if the continuity and coercivity constants are $M$ and $\alpha$, then $\Vert u-u_h{\Vert }_V\le (M/\alpha ){\inf }_{v_h\in V_h}\Vert u-v_h{\Vert }_V$. For this symmetric problem the exact best-approximation constant in the energy norm is one. The piecewise-linear finite element interpolation estimate $\Vert u-I_{h}u{\Vert }_{H^1}\le Ch\Vert u{\Vert }_{H^2}$ on regular finite element meshes therefore gives convergence of a numerical method of the Ritz method, with first-order error in the energy norm.

Let $I=(-1,1)$ and use the clamped second-order Sobolev spaceHere weak derivatives up to order two are square integrable, with Sobolev norm $\Vert v{\Vert }_{H^2}^2=\Vert v{\Vert }_2^2+\Vert v^{\prime }{\Vert }_2^2+\Vert v^{\prime \prime }{\Vert }_2^2$. The boundary values are the continuous traces given by the Sobolev trace theorem, not arbitrary pointwise representatives. This space is the $H^2$ closure of $C_c^{\infty }(I)$, so it is a closed Hilbert space. The derivative boundary conditions are essential, and $H_0^1(I)$ alone is insufficient.

Define the symmetric bilinear form and bounded linear functionalThen the variational problem isFor complex data replace $a$ by its sesquilinear version, take the real part of the forcing term and use absolute squares in $J$; the same argument applies.

Its first variation is $J^{\prime }(u)[v]=a(u,v)-\ell (v)$. Thus stationary points satisfy $a(u,v)=\ell (v)$ for every $v\in H$, the weak formulation of a clamped fourth-order equation. For smooth coefficients and functions, twice integrating the highest term and once the middle term gives the differential equation, with all boundary terms zero because $v=v^{\prime }=0$ at both ends. With only the stated continuous coefficients, this manipulation is instead understood distributionally: compactly supported test functions give $(p{u}^{\prime \prime }{)}^{\prime \prime }-(q{u}^{\prime }{)}^{\prime }+ru=f$ as a distribution. This avoids assuming nonexistent classical derivatives of $p$ or $q$.

Take the positive integer $n\ge 1$ and put $m=n+1\ge 2$. From the preceding estimate and the absent zero Fourier mode, for $t\ne 0$,Here $\zeta$ is the Riemann zeta function, and its displayed series is finite because $m>1$. One may replace the last derivative norm by the given full mixed Sobolev norm to obtain the requested constant depending only on $n$ and $f_0$.

An absolutely summable sequence of Fourier coefficients gives a uniformly and absolutely convergent Fourier series, with supremum bounded by the sum of their absolute values. Its sum agrees almost everywhere with $r_t$ by uniqueness of Fourier coefficients for integrable periodic functions. HenceThis is the uniform phase-mixing bound from velocity derivatives: uniform convergence for the continuous representative of the density, with rate $O(|t{|}^{-n-1})$. No uniform decay of the full phase-space distribution is asserted. The free-transport phase mixing acts by shifting nonzero spatial modes to large velocity frequency. If a convention allows $0\in \mathbb{N}$, that endpoint needs separate assumptions or an argument: the harmonic series in this proof would diverge.

We establish a quantitative almost global existence for wave equations estimate. Take $0<\epsilon \le 1$, fix an integer $m\ge 6$, and use the commutation vector fields for the wave equation from the preceding part. Put $\partial =({\partial }_t,\nabla )$ and define the commuted wave energyAll these L2 norms are finite on any smooth existence interval by finite propagation speed. At $t=0$, the polynomial coefficients of the vector fields are bounded on the fixed compact support of the Cauchy data. Whenever a higher time derivative occurs, use ${\varphi }_{tt}=\Delta \varphi -{\varphi }_t^2$ and its differentiated versions to express it in terms of initial spatial derivatives. Every term contains at least one factor of $\epsilon$; consequentlyfor a constant depending only on finitely many derivatives and the support radius of ${\varphi }_0,{\varphi }_1$.

The commutators $[Z,{\partial }_{\mu }]$ are constant linear combinations of translations. Together with $[\square ,S]=2\square$ and the Leibniz rule, this shows that each commuted source is a finite linear combination of productsThis statement includes the extra copies of the original source produced by the scaling vector field. In each product put the factor with fewer commutations in the Lp norm $L^{\infty }$ and the other in the L2 norm. The lower order is at most $\lfloor m/2\rfloor$. Applying the Klainerman-Sobolev inequality to ${\partial }_{\mu }{Z}^J\varphi$ costs at most two additional commutations; commuting those past ${\partial }_{\mu }$ introduces only lower-order translations. Since $\lfloor m/2\rfloor +2\le m$,The inhomogeneous wave energy estimate now givesLet $T_{\epsilon }={\epsilon }^{-N}$. Use a bootstrap argument with $A_m(t)\le 2C_0\epsilon$ up to the smaller of $T_{\epsilon }$ and the maximal existence time. The energy estimate improves this toFor each fixed $N$,Choose ${\epsilon }_N\le 1$ so that $4C_m{C}_0\epsilon \log (1+{\epsilon }^{-N})\le 1/2$ for every $0<\epsilon \le {\epsilon }_N$. Then $A_m(t)\le (3/2)C_0\epsilon$, a strict improvement. A continuity argument closes the bootstrap argument.

Finally, the translation terms in $A_m$ control ordinary spatial Sobolev norms of $\partial \varphi$. The missing L2 norm of $\varphi$ satisfiesThus the full local-existence Sobolev norms remain bounded on this finite interval. The smooth continuation criterion for semilinear wave equations extends the solution past any finite endpoint before $T_{\epsilon }$. To see smooth persistence explicitly, the tame Sobolev product estimate gives $\Vert {\varphi }_t^2{\Vert }_{H^k}\le C_k\Vert {\varphi }_t{\Vert }_{\infty }\Vert {\varphi }_t{\Vert }_{H^k}$. Ordinary differentiated wave energy estimates therefore bound each higher derivative energy by its initial value times $\exp (C_k{\int }_0^t\Vert {\varphi }_t(s){\Vert }_{\infty }ds)$. This is finite on the interval already controlled by the base commuted wave energy; no separate ${\epsilon }_N$ is needed for each derivative order. ThereforeFor $\epsilon =0$ the zero solution is global. The same energy estimate in fact permits an exponential lower bound for the lifespan, which is stronger than any fixed inverse power.

Small data global regularity for wave maps holds in four spatial dimensions. Smallness is measured in sufficiently high weighted Sobolev norms relative to a constant map, with localized compatible Cauchy data. The vector field method for wave equations gives derivative decay $(1+t{)}^{-3/2}$. This is time-integrable, so commuted wave energy estimates close a small-data bootstrap argument for the derivative-quadratic semilinear wave equation. Higher regularity persists. No smallness of energy alone is asserted.

Small data global regularity for wave maps also holds in three spatial dimensions. For localized compatible data small in high weighted Sobolev norms, the key is the classical null condition for wave equations. Each derivative contraction is a null form for wave equations, vanishing for parallel null derivatives. The vector field method for wave equations exploits derivatives tangent to the light cone and weighted energy estimates to close the bootstrap argument; ordinary $(1+t{)}^{-1}$ decay alone is insufficient.

For $I=(-a,a)$, the first-order Sobolev space isThe weak derivative $v=u^{\prime }$ is unique as an equivalence class modulo equality almost everywhere. We use the standard Sobolev norm $\Vert u{\Vert }_{H^1(I)}^2=\Vert u{\Vert }_{L^2(I)}^2+\Vert u^{\prime }{\Vert }_{L^2(I)}^2$. Pointwise values of an arbitrary representative have no intrinsic meaning; the one-dimensional Sobolev representative supplies the continuous representative used later.

Reflect $u$ across the two endpoints of $I$ into a slightly larger interval, giving an $H^1$ extension, and multiply by a smooth cutoff function equal to one near $\overline{I}$. Convolve the resulting function on $\mathbb{R}$ with a mollifier. The convolutions $u_{\epsilon }$ are smooth; differentiation commutes with convolution of the weak derivative. Continuity of translations in $L^2$ proves both $u_{\epsilon }\to u$ and $u_{\epsilon }^{\prime }\to u^{\prime }$ in $L^2(I)$, hence convergence in the Sobolev norm.

To obtain convergence almost everywhere, select a sequence with $\Vert u_n-u{\Vert }_2\le 2^{-n}$. Then $\int {\sum }_n|u_n-u{|}^2<\infty$, so the summands tend to zero almost everywhere. Thus a smooth approximating sequence converges both in $H^1$ and almost everywhere. This is density of smooth functions in a Sobolev space; smooth functions need not vanish at the endpoints.

For smooth periodic initial data, Fourier expansion reduces the equation to $\dot{\hat{u}}(m,t)=-(1+m^2)\hat{u}(m,t)$. Therefore define the Fourier multiplier solution operators byTheir Sobolev norm satisfies $\Vert S(t)u_0{\Vert }_{H^s}^2\le e^{-2t}\Vert u_0{\Vert }_{H^s}^2\le \Vert u_0{\Vert }_{H^s}^2$, so they extend to contractive linear operators on $H_{\mathrm{per}}^s$. Multiplication of the Fourier factors gives $S(t+s)=S(t)S(s)$ and $S(0)=I$. Finallyby dominated convergence theorem, dominated by four times the summable norm terms. This proves a strongly continuous contraction semigroup for every requested $s$; the semigroup property extends continuity from zero to other times. For smooth data, termwise derivatives converge and recover the original Cauchy problem for a partial differential equation uniquely.

For an integer $s\ge 0$, the Sobolev space iswhere $D^{\beta }$ is a weak derivative. For $p<\infty$ one may use the Sobolev norm $({\sum }_{|\beta |\le s}\Vert D^{\beta }u{\Vert }_p^p{)}^{1/p}$; for $p=\infty$ use the maximum of the finitely many essential-supremum norms.

For $1\le p<d$, the Sobolev inequality is $\Vert u{\Vert }_{p^{\ast }}\le C_{d,p}\Vert \nabla u{\Vert }_p$, with Sobolev conjugate exponent $p^{\ast }=dp/(d-p)$. For $d<p<\infty$, Morrey's inequality supplies a continuous representative satisfyingAt $p=\infty$ the representative is Lipschitz continuous. At the critical exponent $p=d>1$, first-order Sobolev regularity gives every finite $L^q$ embedding for $q\ge d$, with an inhomogeneous norm, but generally no $L^{\infty }$ embedding. The one-dimensional endpoint $W^{1,1}(\mathbb{R})\hookrightarrow L^{\infty }(\mathbb{R})$ is an exception.

For the proof of Morrey's inequality, start with a smooth $u$ and write $u_{B(x,r)}$ for its average on a ball. Averaging the fundamental theorem of calculus along a line segment and changing radial variables givesThe last step is the Holder inequality; integrability of the kernel to power $p^{\prime }$ is exactly $(d-1)p^{\prime }<d$. For $r=|x-y|$, translate the averaging ball along the segment from $x$ to $y$. The fundamental theorem of calculus along a line segment and the Holder inequality giveCombining the two point-to-average bounds and this average-to-average bound proves the required Hölder estimate. The point-to-average bound with $r=1$, together with $|u_{B(x,1)}|\le |B_1{|}^{-1/p}\Vert u{\Vert }_p$, gives the supremum estimate. Density of smooth functions in a Sobolev space then gives a uniformly convergent sequence of smooth representatives, preserving both bounds. For $p=\infty$, mollification gives the Lipschitz version.

For the decay conclusion assume $d<p<\infty$. The representative is uniformly continuous. If $|u(x_j)|\ge 2\theta >0$ along points escaping to infinity, the Hölder bound gives a radius $\rho >0$, independent of $j$, on which $|u|\ge \theta$. A subsequence has disjoint radius-$\rho$ balls, each contributing at least ${\theta }^p|B_{\rho }|$ to $\Vert u{\Vert }_p^p$, a contradiction. This is uniformly continuous integrable functions vanish at infinity.

The finite-$p$ restriction is necessary. If the printed range includes $p=\infty$, its decay assertion is false: $u\equiv 1$ belongs to $W^{1,\infty }$ but does not tend to zero.

For real smooth functions of compact support in a planar domain $D$, use the Dirichlet inner productChanging the positive normalization factor does not change orthogonality. For complex functions, insert complex conjugation in the second factor to obtain the corresponding Hermitian form.

Let $f:D\to \tilde{D}$ be a conformal bijection, and let $\tilde{u},\tilde{v}$ be smooth functions of compact support on $\tilde{D}$. Its real Jacobian matrix is $Df=|f^{\prime }|R$, where $R$ is a rotation. By the chain rule,The change of variables formula has Jacobian determinant $|f^{\prime }{|}^2$, so this factor cancels:This proves conformal invariance of the planar Dirichlet inner product. By completion it is also an isometry between the corresponding Dirichlet energy spaces. It asserts invariance of the energy form, not of the inhomogeneous Sobolev norm, whose $L^2$ term has a different transformation rule.