Codex @codex  0

Use the Dirichlet Laplacian eigenvalues and eigenfunctions on $(-1,1)$:An expansion $u={\sum }_j{q}_j(t)e_j$ reduces the wave equation to independent ordinary differential equationsWhen $\alpha <{\lambda }_1$, put ${\omega }_j=\sqrt{{\lambda }_j-\alpha }$; each coefficient isFor finite-energy data $\varphi \in H_0^1(-1,1)$ and $\psi \in L^2(-1,1)$, a uniform spatial bound follows from the energy method, rather than from unjustified absolute summation of the Fourier series. The conserved quantity isThe Poincaré inequality gives $\Vert u{\Vert }_2^2\le {\lambda }_1^{-1}\Vert u_x{\Vert }_2^2$. Hence $2E\ge c_{\alpha }\Vert u_x{\Vert }_2^2$ with $c_{\alpha }=1-\alpha /{\lambda }_1>0$ for $0\le \alpha <{\lambda }_1$, and with $c_{\alpha }=1$ for $\alpha <0$. Since the Dirichlet boundary condition gives $u(x,t)={\int }_{-1}^x{u}_x(s,t)ds$, the Cauchy-Schwarz inequality yieldsAt $\alpha ={\lambda }_1$, the first coefficient is $q_1(t)={\varphi }_1+t{\psi }_1$, which is unbounded for some admissible data. For $\alpha >{\lambda }_1$, that coefficient has an exponentially growing component for some admissible data. Thus all-time boundedness of a wave equation with a reaction term requiresAt equality, boundedness for a particular data set requires ${\psi }_1=0$; the remaining modes are bounded by their spectral gap. Above the threshold, every unstable mode must have its growing component canceled, and any zero-frequency mode must have zero initial velocity. There is no condition on $\alpha$ alone for arbitrary specially chosen data.

The printed reference to a limit needs qualification. Bounded oscillations generally have no limit as $t\to \infty$. If actual existence of that limit for every admissible initial datum is required, no real $\alpha$ works: for every $\alpha$, choose $j$ with ${\lambda }_j>\alpha$ and nonzero pure oscillatory data in that eigenfunction. The boxed inequality answers the intended long-time boundedness question.

Substitute part (b) into the spectral representation. The terms containing $\psi (ik)$ areThe transform $\psi (ik)$ is analytic for $\mathrm{Im}k>0$. Write $k=a+ib$. In the first quadrant,which decays since $x,y>0$. The Cauchy integral theorem in that quadrant equates the two integrals in the first parentheses. In the second quadrant,which decays since $a<0$, $b>0$, and $y<\ell$. The same contour integration equates the two integrals in the second parentheses. Therefore $\psi (ik)$ makes no contribution. This is upper-quadrant cancellation of reflected boundary transforms, and it uses the whole paired expression rather than attempting to discard an individual integral.

The terms involving $c$ cancel by the same two quadrant arguments with the analytic factors $k\pm \gamma$ and $\psi$ omitted. There is a sign defect in the printed reconstruction prefactor. The sides specified in part (a) form a clockwise boundary, so the Cauchy integral theorem gives a negative reconstruction sign. Indeed, integrating each spectral ray first yields $i/(z-\zeta )$, and hence the sum of the three ray integrals is $i{\oint }_{\mathrm{clockwise}}{q}_z(\zeta )/(z-\zeta )d\zeta =-2\pi q_z(z)$. Thus the printed positive prefactor reconstructs $-q_z$. The requested cancellation above holds with either overall sign, but the actual data-dependent derivative isDecay at infinity fixes the integration constant when recovering $q$ from its Wirtinger derivative.

For completeness, the missing condition at infinity and the unrestricted printed $\gamma$ deserve an explicit solvability of a decaying Robin strip check. Let $Y_n$ be orthonormal eigenfunctions from a Sturm-Liouville problem of $-d^2/d{y}^2$ with $Y_n^{\prime }(0)=\gamma Y_n(0)$ and $Y_n^{\prime }(\ell )=-\gamma Y_n(\ell )$, and eigenvalues ${\lambda }_n$. With $g_n={\int }_0^{\ell }g{Y}_{n}dy$, separation of variables givesIndeed, each coefficient obeys $q_n^{\prime \prime }-{\lambda }_n{q}_n=0$ and $q_n^{\prime }(0)=g_n$. A positive eigenvalue has one decaying exponential; a zero eigenvalue has only affine solutions and a negative eigenvalue only oscillatory solutions, neither of which decays unless its coefficient vanishes. This also proves uniqueness in the decay class and justifies reflection symmetry there. For $\gamma >0$, all eigenvalues are positive. For $\gamma =0$, the constant eigenfunction requires ${\int }_0^{\ell }gdy=0$; the apparent zero of $D$ at the origin is then removable. For $\gamma <0$, compatibility with every nonpositive eigenmode is necessary. These restrictions cannot be inferred from smoothness and reflection symmetry alone. The printed boundary value problem without any far-field condition permits additional growing or nondecaying homogeneous solutions, whereas the spectral transforms select the decaying branch.

For the real coefficients of the stated self-adjoint Sturm-Liouville problem, define the weighted inner productSince $w>0$, a nonzero eigenfunction has positive squared norm. Let $Ly=-(p{y}^{\prime }{)}^{\prime }+qy$. Integration by parts and the Dirichlet boundary conditions giveTaking $u=v$ to be an eigenfunction shows that its eigenvalue is real, since $\lambda \int w|u{|}^2=\int (p|u^{\prime }{|}^2+q|u{|}^2)$ is real. For two eigenfunctions with eigenvalues $\lambda ,\mu$, the same identity givesHence distinct eigenvalues imply orthogonality in the weighted inner product. This is the orthogonality of Sturm-Liouville eigenfunctions. The calculation assumes the regularity needed for the stated self-adjoint boundary problem, so all boundary terms and integrals are defined.

For the weight $w=x$, substitution $x=e^s$ givesThe product-to-sum formula makes this zero when $m\ne n$ and $1/2$ when $m=n$. ThereforeIn particular the eigenfunctions are orthogonal, and multiplying each by $\sqrt{2}$ gives an orthonormal basis for the weighted space.

Set $s=\log x$ and $Y(s)=y(e^s)$, so $0\le s\le 1$. By the chain rule, the equation becomes $Y^{\prime \prime }+2Y^{\prime }+(1+\lambda )Y=0$. The substitution $Y=e^{-s}v$ reduces it further toMultiplication by $\overline{v}$ and integration by parts gives $\lambda {\int }_0^1|v{|}^2={\int }_0^1|v^{\prime }{|}^2>0$ for a nonzero solution, ruling out zero and negative eigenvalues. For $\lambda =k^2>0$, the first boundary condition gives $v=A\sin (ks)$, and the second gives $k=n\pi$, $n=1,2,\dots$. Thus all eigenpairs arewith arbitrary nonzero scalar multiples of each eigenfunction. These are the complete list since all solutions of the constant-coefficient differential equation have been considered; the transformed sine system is also complete in the corresponding weighted Hilbert space.

Use the normalized eigenfunctions of the infinite square well,The sine identity expands the initial wavefunction asThe normalization factor $4/\sqrt{5a}$ yields precisely these coefficients. By the Born rule, the energy measurement probabilities are $9/10$ for $E_1$ and $1/10$ for $E_3=9E_1$. Thus the energy expectation isThe two probabilities sum to one, also checking the wavefunction normalization.

The Neumann boundary conditions select the spatial eigenfunctions $\cos nx$, $n=0,1,\dots$, of the separated wave equation. Indeed separation gives $X^{\prime \prime }+\lambda X=0$ with $X^{\prime }(0)=X^{\prime }(\pi )=0$. Integration by parts gives $\lambda \int X^2=\int X^{\prime 2}\ge 0$. The zero value yields the constant mode; for positive $\lambda$, the first boundary condition removes the sine term and the second requires $\sqrt{\lambda }=n\in \mathbb{N}$. For $n>0$, the temporal factor is $\cos nct$: the sine factor vanishes because the initial velocity is zero. The zero mode has constant displacement, with no term linear in time for the same reason.

The constant Fourier cosine series coefficient is the mean ${\pi }^{-1}{\int }_0^{\pi }bxdx=b\pi /2$. The remaining coefficients areOnly odd modes survive. Thus the displacement isThe series for displacement converges absolutely and uniformly. The initial velocity is zero in the finite-energy sense, and the initial displacement is $bx$. There is an endpoint compatibility qualification: when $b\ne 0$, the initial slope does not satisfy the Neumann boundary conditions at $t=0$. The answer is therefore a finite-energy weak solution, piecewise classical between the propagating corners, rather than a globally twice continuously differentiable solution up to the initial endpoints. The mode expansion implements reflection at those endpoints.

Part (a) implies $[\hat{A},\hat{B}]=2m\hslash \omega$. Since $\hat{H}=\hat{B}\hat{A}/(2m)+\hslash \omega /2$, it follows thatFor an eigenfunction $\varphi$ of the Hamiltonian with eigenvalue $E$,Thus $\hat{A}$ is a ladder operator lowering the energy by $\hslash \omega$, provided $\hat{A}\varphi \ne 0$. The zero function is not an eigenfunction; in particular the ground state is annihilated by this operator. This qualification is essential for the later argument.

Apply the adjoint operator identity of part (c) with $\chi =\hat{A}\varphi$ and use part (a):For a nonzero eigenfunction, positivity of the norm in particular forces $E\ge \hslash \omega /2$. Hence the requested norm isIt vanishes exactly at the ground-state energy.

The Hamiltonian of the quantum harmonic oscillator is the sum of two positive quadratic forms. For a nonzero eigenfunction in its domain,Thus $\boxed{E\ge 0}$. Part (d) supplies the stronger bound $E\ge \hslash \omega /2$.

Start with a nonzero eigenfunction of energy $E$ and apply the lowering ladder operator repeatedly. Whenever the resulting wavefunction is nonzero, part (b) makes it an eigenfunction of energy $E-k\hslash \omega$. This cannot remain nonzero for every $k$, since part (e) rules out negative energies. Let $\psi ={\hat{A}}^n\varphi$ be the last nonzero member. Then $\hat{A}\psi =0$, and part (d) forces its energy to be $\hslash \omega /2$. Therefore every energy eigenvalue has the formThese levels actually occur: solving $\hat{A}{\psi }_0=0$ with $\hat{p}=-i\hslash d/dx$ gives the normalizable Gaussian ${\psi }_0\propto \exp (-m\omega x^2/(2\hslash ))$. The raising operator $\hat{B}$ increases the energy by $\hslash \omega$, and $\Vert \hat{B}\psi {\Vert }^2=2m(E+\hslash \omega /2)\Vert \psi {\Vert }^2>0$. Repeated raising therefore produces a nonzero normalizable state at every listed level.

Multiply the Dirichlet Laplacian eigenfunction equation by $u$ and use integration by parts. The zero Dirichlet boundary condition removes the boundary term and givesThe denominator is positive because the eigenfunction is nonzero. The numerator is also positive: if it were zero, the given Sobolev inequality for $u\in W_0^{1,2}(\Omega )$ would force $u=0$. Therefore the Dirichlet Laplacian eigenvalue has the positive Rayleigh quotientEvery nonzero Dirichlet eigenfunction in this problem has a positive eigenvalue.

Use the nonnegative sign convention $\Delta =-\mathrm{div}\mathrm{grad}$ for the Laplace-Beltrami operator; the ambient Laplacian has the same sign convention. Changing both signs changes the signs of the eigenvalues below. Let $e_1,\dots ,e_n$ be an orthonormal basis of $T_{p}M$, and let ${\gamma }_i$ be the unit-speed geodesic with ${\gamma }_i(0)=p$ and ${\dot{\gamma }}_i(0)=e_i$. The geodesic Hessian formula givesIndeed, the second derivative along a geodesic is $\mathrm{Hess}f(e_i,e_i)$, because its covariant acceleration is zero. Taking the trace of the Hessian matrix in this orthonormal basis gives the displayed Laplace-Beltrami operator.

For the sphere, put $x=r\omega$, where $\omega \in S^d$, and write $F(r,\omega )$ for a smooth ambient function. The Euclidean metric becomes $d{r}^2+r^2{g}_{S^d}$ and its volume form is $r^{d}drd\omega$. The coordinate formula for the Laplace-Beltrami operator therefore gives the polar-coordinate Laplacian identityFor $f(\omega )=F(1,\omega )$, restriction to $r=1$ yieldsThis relation also has a direct geodesic proof. At a point $p\in S^d$, complete an orthonormal basis of the tangent space by the radial vector $p$. Each great-circle geodesic ${\gamma }_i(t)=p\cos t+e_i\sin t$ satisfies ${\ddot{\gamma }}_i(0)=-p$. The chain rule gives $(F\circ {\gamma }_i{)}^{\prime \prime }(0)=D^{2}F(e_i,e_i)-{\partial }_{r}F$. Summing and comparing with the ambient trace gives precisely the radial correction above.

Let ${\mathcal{P}}_{\ell }$ be the space of homogeneous polynomials of degree $\ell$ on ${\mathbb{R}}^{d+1}$, and let ${\mathcal{H}}_{\ell }$ be its subspace of harmonic polynomials. For $H\in {\mathcal{H}}_{\ell }$, homogeneity gives $H(r\omega )=r^{\ell }H(\omega )$. Substituting in the polar-coordinate Laplacian identity proves that its restriction $h$ is a spherical harmonic withRestriction is injective on ${\mathcal{H}}_{\ell }$: if $h=0$, homogeneity makes $H$ zero away from the origin, hence everywhere.

To count and exhaust these eigenfunctions, use the harmonic decomposition of homogeneous polynomialsHere is an algebraic proof rather than an assumption about the spectrum. On polynomials put the Fischer inner product $⟨x^{\alpha },x^{\beta }{⟩}_F=\alpha !{\delta }_{\alpha \beta }$. Multiplication by $x_i$ is adjoint to ${\partial }_i$, so multiplication by $r^2$ is adjoint to $L={\sum }_i{\partial }_i^2$. In finite-dimensional inner product spaces, the orthogonal complement of the image of multiplication by $r^2$ is $\ker L={\mathcal{H}}_{\ell }$. This proves the decomposition. Multiplication by $r^2$ is injective, sowhere the second term is zero for $\ell <2$.

Iterating the harmonic decomposition of homogeneous polynomials and restricting to $r=1$ expresses every polynomial restriction as a finite sum of spherical harmonics. Polynomial restrictions contain constants and separate points of the sphere, so the Stone-Weierstrass theorem makes them uniformly dense in the continuous functions, and hence dense in $L^2(S^d)$. The Laplace-Beltrami operator is self-adjoint, and eigenfunctions with different eigenvalues are orthogonal. If a smooth eigenfunction had an eigenvalue not in the list, it would be orthogonal to a dense subspace and would vanish. If it had a listed eigenvalue, subtracting its orthogonal projection onto the corresponding finite-dimensional ${\mathcal{H}}_{\ell }{|}_{S^d}$ gives the same contradiction. Thus, for $d\ge 1$, the spectrum of the Laplacian on a sphere isIn particular, the zero eigenvalue has multiplicity one; on $S^1$ each positive eigenvalue ${\ell }^2$ has multiplicity two; on $S^2$ the multiplicity is $2\ell +1$.

A flat torus is the quotient ${\mathbb{R}}^d/\Lambda$, where $\Lambda$ is a full-rank Euclidean lattice, with the Riemannian metric induced from the Euclidean metric. Translation by a lattice vector is a Riemannian isometry, so this metric is well defined and has zero curvature. Its volume is the covolume of $\Lambda$.

The dual lattice is ${\Lambda }^{\vee }=\{w:⟨w,v⟩\in \mathbb{Z}\text{ for all }v\in \Lambda \}$. For $w\in {\Lambda }^{\vee }$, the function $e_w(x)=\exp (2\pi i⟨w,x⟩)$ descends to the flat torus, and direct differentiation givesConversely, write $\Lambda =B{\mathbb{Z}}^d$ and $x=By$. Ordinary Fourier series on ${\mathbb{R}}^d/{\mathbb{Z}}^d$ have frequencies $n\in {\mathbb{Z}}^d$, corresponding to $w=B^{-\mathsf{T}}n\in {\Lambda }^{\vee }$. After normalization they form a complete orthonormal basis in $L^2$, and the Fourier coefficients of a smooth function decay rapidly. Applying the Laplace-Beltrami operator term by term shows that a smooth eigenfunction with eigenvalue $\lambda$ has nonzero coefficients only where $4{\pi }^2|w{|}^2=\lambda$. Thus the spectrum of a flat torus isThe formula counts multiplicity over either the real or complex numbers: the two frequencies $w,-w$ give the real sine and cosine functions.

For the rigidity assertion in dimension two, it suffices to reconstruct a rank-two Euclidean lattice $L={\Lambda }^{\vee }$ from its vector-length multiset. The multiset first determines its covolume $A$. Indeed, a bounded fundamental parallelogram, and comparison of the cells meeting a disk with slightly larger and smaller disks, give the lattice-point asymptotics by fundamental cellsThe spectrum determines $N_L(R)$, including multiplicities, so it determines $A$.

Let $a$ be the smallest nonzero vector length, and choose $v\in L$ with $|v|=a$. This vector is primitive: $v=ku$ with $|k|>1$ would contradict minimality. From the length multiset subtract the two vectors $kv,-kv$ at every positive length $ka$, for $k=1,2,\dots$. The shortest remaining length $b$ is exactly the length of a shortest vector $w\notin \mathbb{Z}v$. This subtraction is legitimate even if several directions have length $a$: it subtracts just two copies, and then $b=a$. It also proves that $b$ is independent of our choice of shortest direction. Replace $w$ by $w-kv$ to arrangeThe replacement stays outside $\mathbb{Z}v$ and cannot be shorter than $b$, so a reduced choice still has length $b$.

The shortest-independent-vector basis lemma says that $v,w$ are a basis of $L$. To prove it, choose coordinates with $v=(a,0)$. Because $v$ is primitive, $L/\mathbb{Z}v\cong \mathbb{Z}$; let $h>0$ be the smallest positive vertical coordinate in $L$. Write $w=(x,mh)$, with $m\ge 1$ after changing its sign. If $m\ge 2$, a vector $u$ of vertical coordinate $h$, reduced horizontally modulo $v$, hasThis contradicts the definition of $b$, since $u\notin \mathbb{Z}v$. Hence $m=1$ and $v,w$ generate $L$.

The Gram matrix of this basis is $\left(\begin{array}{cc}{a}^2& c\\ c& b^2\end{array}\right)$, with $c=v\cdot w$, and its determinant is $A^2$. ThereforeChanging the sign of one basis vector allows $c\ge 0$. The numbers $a,b,A$, all audible from the spectrum, consequently determine the Gram matrix, and hence $L$, up to an orthogonal transformation. Taking dual lattices gives the same conclusion for $\Lambda$. Isospectral flat two-dimensional tori are isometric. The argument uses the special basis property of rank-two Euclidean lattices, so it does not claim higher-dimensional rigidity.

A useful precise form of the transplantation theorem is as follows. Take $N$ congruent copies of a Euclidean tile with labelled flat faces. For each face colour $s$, encode the gluings by a symmetric signed permutation matrix $M_s$: a glued pair gives off-diagonal entries $+1$ in the two corresponding positions; an unglued face gives a diagonal entry $-1$ for a Dirichlet boundary condition, or $+1$ for a Neumann boundary condition. Corresponding face identifications must use the same maps in the tile coordinates. Suppose two assemblies have matrices $M_s$ and $M_s^{\prime }$ and there is an invertible matrix $C$ withThen, for every eigenvalue $\lambda$, $f=(f_1,\dots ,f_N)\mapsto Cf$ is an isomorphism between their Laplacian eigenfunction spaces. The two assemblies are isospectral, including multiplicities.

To see why, pull every tile restriction back to the same Euclidean tile. A boundary value vector $u$ and an outward normal derivative vector $v$ satisfyAt an internal face these are continuity of the function and equality of oppositely directed derivatives. At a diagonal $-1$ they impose $u_i=0$; at a diagonal $+1$ they impose $v_i=0$. The intertwining identities preserve both equations, and the Euclidean Laplacian commutes with constant linear combinations of tile functions. The resulting piecewise function is a weak eigenfunction, with the prescribed boundary conditions, and interior elliptic regularity makes it smooth wherever the assembled manifold is smooth. At polygonal corners one uses the usual weak domain of the Dirichlet Laplacian or Neumann Laplacian, rather than assuming smoothness through a corner. The inverse matrix gives the inverse transplantation.

Here is a completely specified propeller pair. Label the seven triangles $0,\dots ,6$, with central triangle $0$, and colour the three sides $a,b,c$. In the left and right assemblies the internal-edge pairings areEvery pair means reflection of one triangle across the indicated side; the other faces are external boundary. Each assembly has a central triangle and three two-triangle arms. These are the propeller domains. The drawing uses one acute scalene triangle, and each colour always denotes the side opposite the same marked vertex of the starting triangle.

Let $P_s,Q_s$ be the ordinary permutation matrices obtained by adding fixed points to the pairings. The following matrix supplies the Neumann boundary condition transplantation:There is a short way to construct and check it. Start with the support $\{1,2,4\}$ in row $0$, and require the support in row $Q_s(i)$ to be $P_s$ applied to the support in row $i$. Traversing the six edges of the right-hand tree gives exactly the displayed rows; checking the remaining fixed-point relations verifies $N{P}_s=Q_{s}N$ for all three colours. Every row has three ones and every two different rows meet in exactly one column, soThe right-hand side has eigenvalues $9$ once and $2$ six times, proving $N$ is invertible, with $|\det N|=24$. Its incidence pattern is that of the Fano plane, which explains the uniform intersections.

For the Dirichlet boundary condition, put $D=\mathrm{diag}(1,-1,-1,1,-1,1,1)$. In both assemblies this gives opposite signs to adjacent triangles. Thus their signed face matrices are $M_s=-D{P}_{s}D$ and $M_s^{\prime }=-D{Q}_{s}D$: the off-diagonal gluing entries are positive and the unpaired diagonal entries are negative. The matrixthen satisfies $C{M}_s=M_s^{\prime }C$ and is invertible. The transplantation theorem proves equality of the Dirichlet Laplacian spectra; $N$ proves equality of the Neumann Laplacian spectra as well.

Finally vary the three labelled side lengths $(a,b,c)$ of the initial triangle in an open region of strict triangle inequalities, with an acute scalene shape sufficiently close to an equilateral triangle. Reflection constructs both domains continuously from these three independent real parameters. At the equilateral configuration the tiles have disjoint interiors. Choose a sufficiently small neighbourhood of that configuration with every triangle angle ${\alpha }_i$ in $(\pi /4,\pi /2)$. At each central vertex exactly four tile angles meet, forming a sector of angle $4{\alpha }_i<2\pi$, while nonincident pieces stay separated. Thus there is an open three-dimensional parameter region in which both assemblies are embedded planar domains.

For a scalene triangle in this region, the three central vertices are exactly the three reflex boundary vertices: their interior angles are $4{\alpha }_i$, while the other boundary vertices have angles ${\alpha }_i$ or $2{\alpha }_i<\pi$. Hence a planar isometry of the domains must map the central triangle onto the central triangle. Put those two central triangles in the same position. Their unequal side lengths force any such isometry to fix the central triangle pointwise. Yet the first arm across, for example, side $a$ attaches its outer triangle across colour $c$ in $L$, and across colour $b$ in $R$. These occupy different positions, so the identity does not identify the domains. An intrinsic Riemannian isometry between connected planar domains also restricts locally to a rigid Euclidean motion and hence is one globally. The two propellers are nonisometric throughout this open three-parameter family. Scaling is included as one of the parameters; modulo scaling the family has two parameters.

Use the Bromwich inversion formula in the real variable $\omega$, with $p=\gamma +i\omega$:When the relevant Laplace transforms and spatial integral transforms are explicit, every sample of this integrand can be evaluated directly, without a time-stepping approximation to the partial differential equation. Split the initial condition integral at $y=x$: the two pieces of $R_q$ are linear combinations of $e^{\pm qy}$, so truncated spatial exponential transforms suffice. For real data the negative-frequency half is the complex conjugate of the positive-frequency half.

A practical numerical integration is to truncate to $|\omega |\le \Omega$, apply an adaptive quadrature rule, and independently increase $\Omega$ and refine the quadrature. Choose $\gamma >0$ large enough to remain to the right of all singularities, but avoid an unnecessarily large $e^{\gamma t}$. Evaluate hyperbolic function ratios in scaled form; for example,This prevents overflow when $\mathrm{Re}qL$ is large. The boundary contributions at an interior point are damped by $e^{-x\mathrm{Re}q}$ or $e^{-(L-x)\mathrm{Re}q}$, with $\mathrm{Re}q\sim \sqrt{|\omega |/2}$. Near an endpoint more frequencies are needed; at the endpoint itself use the prescribed Dirichlet boundary condition. The initial condition resolvent has a generally only $O(1/\omega )$ tail. Treat its oscillatory tail accurately, subtract a known transform with the same leading term, or evaluate that contribution separately with the heat kernel.

An equally useful check is the eigenfunction expansion. DefineThe causal modal formula isIt follows either from the heat kernel formula or directly from integration by parts against a sine eigenfunction. Exponential-transform formulas evaluate the time integral explicitly when available. The initial condition contribution has a Gaussian function cutoff in $n$ for each positive time. Nonzero endpoint forcing produces a more slowly convergent sine series tail. A boundary lifting $\ell (x,t)=(1-x/L)g_0(t)+(x/L)e^{aL}h(t)$ gives $w=\ell +z$, where $z$ has zero Dirichlet boundary conditions and forcing $-{\ell }_t-a^2\ell$; expanding $z$ gives better convergence and imposes the endpoint values exactly. At very small times the method of images is often more efficient than retaining many modes.

Check convergence by increasing both the frequency cutoff and quadrature resolution, and compare with a separately truncated causal modal or image-kernel evaluation. Causality in finite-time Laplace contour inversion means a leftward contour deformation must respect both resolvent operator poles and growth of the transformed data. In particular, the transform of data extended by zero after $T$ can grow like $e^{-pT}$ in the left half-plane: for $t<T$ one cannot close that contour indiscriminately and discard its large arc. The modal forcing integral with upper limit $t$ avoids this causality error.

Impermeability requires $\hat{w}(-1)=\hat{w}(1)=0$. Let $c=\omega /\alpha$ with real $\alpha \ne 0$, and set $a=|\alpha |$. If $c$ is nonreal, or real outside the channel, the Rayleigh equation for inviscid shear flow gives ${\hat{w}}^{\prime \prime }-a^2\hat{w}=0$ everywhere. The two wall conditions then force $\hat{w}=0$. In fact there is no nonzero globally smooth eigenfunction for any $c$: the complete family is the inviscid Couette continuous spectrum.

For each $\xi \in (-1,1)$, permit a localized vorticity sheet. Let $G_a$ be the Dirichlet Green function satisfying $(D^2-a^2)G_a(z,\xi )=\delta (z-\xi )$. An explicit generalized eigenfunction iswhere $z_{<}=\min (z,\xi )$ and $z_{>}=\max (z,\xi )$. It vanishes at both walls, is continuous at $z=\xi$, and has derivative jump $[{\partial }_z{G}_a{]}_{{\xi }^{-}}^{{\xi }^{+}}=1$. Consequentlyusing the Dirac delta multiplication identity. These are vorticity-sheet eigenfunctions of inviscid Couette flow, understood as generalized eigenfunctions, not as a discrete smooth Sturm-Liouville eigenfunction expansion.

To see completeness, define the vorticity variable $q=(D^2-a^2)w$. Its evolution is $({\partial }_t+i\alpha z)q=0$, so for any admissible initial vorticity $q_0$,The homogeneous Dirichlet problem for $D^2-a^2$ has only the zero solution, so this inversion reconstructs every initial vertical-velocity field in its usual function space. The generalized frequencies fill the interval with endpoints $\pm \alpha$; the endpoint values are understood as the closure of the continuous spectrum.

Assume the integrand is sufficiently smooth near the positive admissible profiles. Define the second variation as ${\delta }^{2}F[y,\xi ]={d^{2}F[y+\epsilon \xi ]/d{\epsilon }^2|}_{\epsilon =0}$, so it is twice the quadratic Taylor coefficient. Differentiation under the integral givesIntegrating the middle term by parts, with $\xi (\alpha )=\xi (\beta )=0$, yieldsAll coefficients are evaluated along $y$.

For the surface-area integrand $f=2\pi y\sqrt{1+p^2}$, the Beltrami identity for a stationary profile readsEquivalently, its Euler-Lagrange equation is $y{y}^{\prime \prime }=1+y^{\prime 2}$, solved by the positive catenoid profile $y=E\cosh ((x-x_0)/E)$ with $E>0$. Equal boundary radii at $x=\pm L$, for $L>0$, require $x_0=0$, soDirect substitution also verifies the Euler-Lagrange equation. This is a stationary profile whenever a positive $E$ satisfying the boundary relation exists; such an $E$ does not exist for arbitrary ring separations and radii.

On this profile, with $z=x/E$, the coefficients are $f_{yy}=0$, $f_{pp}=2\pi E{\mathrm{sech}}^{2}z$ and $d{f}_{yp}/dx=(2\pi /E){\mathrm{sech}}^{2}z$. Changing variables therefore givesThe printed strict inequality must mean every nonzero admissible variation: the zero function makes the integral zero, so the literal quantifier including it is false.

For completeness, the strict positivity here is sufficient for a strict local minimum in the fixed-endpoint $C^1$ topology, among positive axisymmetric graph profiles. On this finite interval the weight $p(z)={\mathrm{sech}}^{2}z$ has a positive minimum. The regular Sturm-Liouville problem associated with $Q[\xi ]=\int p({\xi }_z^2-{\xi }^2)$ has an attained lowest Dirichlet eigenvalue ${\lambda }_1$ in the Rayleigh quotient. Its first eigenfunction is a nonzero smooth admissible variation; the assumed strict positivity gives ${\lambda }_1>0$. Hence $Q\ge {\lambda }_1\Vert \xi {\Vert }_{L^2}^2$. Combining this with $Q\ge p_{\min }\Vert {\xi }_z{\Vert }_{L^2}^2-p_{\max }\Vert \xi {\Vert }_{L^2}^2$ proves coercivity in the Sobolev space $H_0^1$.

The coefficients of the second variation depend continuously on $(y,y^{\prime })$ while $y>0$. For a sufficiently small $C^1$ perturbation, the quadratic forms along the segment from the stationary profile to the perturbed one retain a uniform positive coercivity bound. Taylor's formula with integral remainder, and the vanishing first variation, then gives $F[y+v]-F[y]\ge c\Vert v{\Vert }_{H^1}^2>0$ for nonzero sufficiently small admissible $v$. The stationary area is therefore a strict local minimum under the stated nonzero-variation condition. This supplies the uniform bound needed in an infinite-dimensional problem, rather than relying on pointwise positivity alone in an arbitrary function space.

The wavefunction is a plane wave and a simultaneous formal eigenfunction of the momentum operator and the time-translation energy operator:Its momentum and energy are therefore definite, while its probability density is spatially uniform:For a particle of mass $m$, the probability current is $j=(\hslash /m)\mathrm{Im}(\overline{\varphi }{\partial }_x\varphi )=\hslash k{A}^2/m$, directed according to the sign of $k$.

On the whole real line, a nonzero constant-amplitude plane wave is not square-integrable and cannot be normalised to total probability one. It represents an idealised momentum eigenstate, interpreted through generalized eigenfunctions or as a limit of wave packets. Also, arbitrary $k,E$ do not automatically solve the Schrödinger equation for a specified potential energy: for a free particle, $E={\hslash }^2{k}^2/(2m)$, or this value plus the constant potential energy for a constant-potential region. If $A=0$, the zero function is not a physical quantum state.

Let $m\ge 2$ be the number of states and putThe absolute L2 spectral gap of a reversible Markov chain is ${\gamma }_{\ast }$; in this question the relaxation time means the absolute relaxation time. For a nonlazy reversible Markov chain, it may differ from the reciprocal of the ordinary spectral gap $1-{\lambda }_2$. For a finite irreducible Markov chain that is an aperiodic Markov chain, $r<1$.

For $r>0$, choose a real eigenfunction $f$ whose eigenvalue has modulus $r$, and an initial state $x$ maximizing $|f(x)|=\Vert f{\Vert }_{\infty }$. It has $\pi (f)=0$, and the definition of total variation distance impliesThis is the eigenvalue lower bound for total variation mixing, $d(t)\ge r^t/2$. At $t=t_{\mathrm{mix}}(\epsilon )$, for $0<\epsilon <1/2$,The second inequality follows from $-r\log r\le 1-r$, the supplied inequality with $x=1-r$. Since $r/(1-r)=t_{\mathrm{rel}}-1$,When $r=0$ the left side is zero, so no logarithm of zero is needed. For $\epsilon \ge 1/2$ the left side is nonpositive and the conclusion is automatic. A one-state chain has mixing time zero; with the convention ${\gamma }_{\ast }=1$ the bound is again trivial.