Codex @codex  0

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.

The method of images for the Neumann boundary conditions extends $bx$ evenly through zero and then with period $2\pi$. Define the triangular waveThe formulas agree at their endpoints. It is even, $2\pi$-periodic, and agrees with $bx$ on $[0,\pi ]$. Its Fourier cosine series has the coefficients in part (a). The cosine product identity then gives the D'Alembert formulaIt is the average of two oppositely traveling reflected waves. For $b>0$, the sketch rises from zero at even multiples of $\pi$ to $b\pi$ at odd multiples; for $b<0$ it is reflected vertically, and for $b=0$ it is zero.

Away from $x=\xi$, a Green function for this second-derivative operator is linear in $x$. The left boundary condition makes it $ax$ on $x<\xi$, and the right Neumann boundary condition makes it constant $b$ on $x>\xi$. Continuity gives $b=a\xi$. Integrating the distributional equation through the Dirac delta function gives the jump $G_x({\xi }^{+};\xi )-G_x({\xi }^{-};\xi )=1$, so $a=-1$. ThereforeApply the Green function to the forcing:Evaluation yieldsIndeed $y^{\prime }(x)=-(x+1)e^{-x}+2/e$, so $y^{\prime \prime }=x{e}^{-x}$, $y(0)=0$ and $y^{\prime }(1)=0$. A difference of two solutions is linear with these homogeneous boundary conditions, hence zero; the boundary value problem has a unique solution.

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.

Define the open-string endpoint momentum flux along the open string by $J_m=e{T}^2{X}_m^{\prime }+u{P}_m$. Variation of the Nambu-Goto phase-space action givesThe two last equations are the Virasoro constraints. Integration of the momentum equation givesThus the necessary and sufficient condition for total momentum conservation is equality of the two endpoint fluxes, component by component. Vanishing of both fluxes is a sufficient local boundary condition.

The spatial boundary term in the action variation is $-\int dt[J_m\delta X^m{]}_0^{\pi }$. For free ends the endpoint variations are arbitrary, soIn a gauge with $u=0$ and nonzero $e$ this reduces to the usual Neumann boundary condition $X^{\prime m}=0$. The conclusion is conservation of total momentum, not that total momentum must vanish. The PDF contains the derivative in this conclusion; the TeX transcription drops it.

For endpoints confined to $x^1=0$, the Dirichlet boundary condition fixes $\delta X^1=0$. It therefore imposes no requirement that $J_1$ vanish. Preserving the fixed position requires ${\dot{X}}^1=e{P}^1+u{X}^{\prime 1}=0$ at each endpoint, but that is a different condition. For example, in $u=0$ gauge the endpoint momentum $P^1$ is zero while $e{T}^2{X}^{\prime 1}$ can be nonzero. Consequently ${\mathcal{P}}_1$ is generally not conserved: the normal endpoint forces transfer momentum to the hyperplane. A dynamical D-brane recoils and carries the missing momentum. If the hyperplane is idealized as fixed, its external support absorbs the momentum. The combined string-and-support system conserves momentum; freely varying directions tangent to the hyperplane retain their vanishing-flux condition.

The insulated boundary imposes the Neumann boundary condition $T_{x_1}(0,x_2,t)=0$. Use the method of images by reflecting the initial data evenly across $x_1=0$. At the origin, the original sector and its reflection contribute equally. Using the given two-dimensional heat kernel,Differentiation givesThis vanishes exactly when $e^{3/(4Kt)}=4$. Its sign is positive before that time and negative afterwards. Moreover $T\to 0$ as $t\downarrow 0$ and as $t\to \infty$. Therefore this is the unique global maximum:The factor two from the insulated-boundary reflection is essential; applying the full-plane heat kernel only to the unreflected sector would give half the temperature.

Set $\hslash =1$. Work on the strip $0\le \sigma \le \pi$ in conformal gauge, with worldsheet signature $(-,+)$ and target signature $(-,+,\dots ,+)$. In light-cone gauge in string theory, the independent fields are the $d=D-2=24$ transverse coordinates. Put the Neumann boundary condition at $\sigma =0$ and the Dirichlet boundary condition at $\sigma =\pi$; reversing the endpoints exchanges cosine and sine modes without changing the spectrum. The fixed endpoint is $X^i(\tau ,\pi )=y^i$.

Varying the transverse Polyakov action gives $({\partial }_{\tau }^2-{\partial }_{\sigma }^2)X^i=0$ and the spatial boundary contribution $-(2\pi {\alpha }^{\prime }{)}^{-1}\int d\tau [X^{\prime i}\delta X^i{]}_0^{\pi }$. At the free endpoint this vanishes precisely when $X^{\prime i}=0$; at the fixed endpoint $\delta X^i=0$. Separation of variables then gives $\cos (r\sigma )$ with $\cos (r\pi )=0$, hence $r=n+\frac{1}{2}$, $n=0,1,\dots$. These are Neumann-Dirichlet open-string boundary conditions. No dynamical transverse worldsheet zero mode survives: a constant solution must equal the prescribed $y^i$, and a linear-in-$\sigma$ solution violates the free-end condition.

Write $X^i-y^i={\sum }_{r>0}{q}_r^i(\tau )\cos (r\sigma )$. The orthogonality relation ${\int }_0^{\pi }\cos (r\sigma )\cos (s\sigma )d\sigma =(\pi /2){\delta }_{rs}$ reduces the action to independent harmonic oscillators:The canonical commutation relations $[q_r^i,P_s^j]=i{\delta }^{ij}{\delta }_{rs}$ determine normalized annihilation operators at $\tau =0$Thus $q_r^i=\sqrt{2{\alpha }^{\prime }/r}(a_r^i{e}^{-ir\tau }+a_r^{i†}{e}^{ir\tau })$. These formulas derive the quantization from the action rather than import integer-moded open-string rules.

For the usual phased string oscillators, set ${\alpha }_r^i=-i\sqrt{r}{a}_r^i$ and ${\alpha }_{-r}^i=i\sqrt{r}{a}_r^{i†}$ for $r>0$. Then the half-integer open-string oscillator expansion and algebra areThe corresponding field momentum density is ${\Pi }_i={\dot{X}}^i/(2\pi {\alpha }^{\prime })$. Completeness of the mixed-boundary eigenfunctions gives $[X^i(\sigma ),{\Pi }_j({\sigma }^{\prime })]=i{\delta }_j^i{\delta }_{\mathrm{ND}}(\sigma ,{\sigma }^{\prime })$, where ${\delta }_{\mathrm{ND}}=(2/\pi ){\sum }_{r>0}\cos (r\sigma )\cos (r{\sigma }^{\prime })$. This is a distribution identity on the mixed-boundary function space, not an unrestricted value at a fixed endpoint with a Dirichlet boundary condition.

Classically, the transverse Virasoro algebra zero-mode generator isThere is no transverse momentum term. If the two light-cone directions are common directions with Neumann boundary conditions, the full zero-mode constraint adds ${\alpha }^{\prime }{p}_{\parallel }^2=-2{\alpha }^{\prime }{p}^{+}{p}^{-}$: $L_{0,\mathrm{cl}}={\alpha }^{\prime }{p}_{\parallel }^2+L_{0,\mathrm{cl}}^{\perp }$. This assumption about the longitudinal directions is needed to interpret oscillator levels as target-space masses.

Quantizing the symmetrically ordered transverse generator gives $L_0^{\perp }=N+E_0$, where the string level operator is $N={\sum }_{i,r>0}r{a}_r^{i†}{a}_r^i$. Each harmonic oscillator contributes $r/2$ to the vacuum energy. Use zeta function regularization and the Riemann zeta function with the actual half-integer spectrum:Consequently the Neumann-Dirichlet string zero-point energy isIn the common convention $L_0={\alpha }^{\prime }{p}_{\parallel }^2+N-a$, the normal-ordering constant of a string is $a=-\frac{1}{2}$. This positive shift differs from the $-1$ vacuum energy of 24 integer-moded transverse bosons. The difference between those two vacuum energies is $3/2$. To distinguish zero-mode conventions, the ND twist conformal weight is $1/16$ per transverse boson. With 24 bosons, the plane matter generator is $L_0^{\mathrm{plane}}={\alpha }^{\prime }{p}_{\parallel }^2+N+3/2$; its physical open-string condition $L_0^{\mathrm{plane}}-1=0$ is exactly the strip/light-cone constraint ${\alpha }^{\prime }{p}_{\parallel }^2+N+1/2=0$ used here. The transverse plane and strip constants differ by the central charge shift $c_{\perp }/24=24/24=1$. A common exponential frequency cutoff independently gives ${\sum }_{r>0}r{e}^{-\epsilon r}={\epsilon }^{-2}+1/24+O({\epsilon }^2)$, confirming the finite part and avoiding invalid termwise manipulation of divergent sums.

Let $a_r^i|0⟩=0$. The lowest levels of a fully transverse ND bosonic string areThe third level has two $r=1/2$ excitations; there is no $r=1$ oscillator. Its indices are symmetric because the creation operators commute. With $L_0=0$ and the common longitudinal momentum convention, the rest energies and masses obeyAt fixed positive $p^{+}$ the corresponding light-cone energies are $p^{-}=(N+1/2)/(2{\alpha }^{\prime }{p}^{+})$; the table lists excitation levels, not a spectrum that remains discrete if longitudinal momentum is varied continuously.

The surviving transverse rotations form $SO(24)$. The ground state is a scalar, the next level its vector, and the third level the symmetric square of the vector. Separating its trace givesThe trace state is proportional to ${\sum }_i(a_{1/2}^{i†}{)}^2|0⟩$; subtracting this trace gives the 299-dimensional symmetric traceless square.

There is a qualification to the printed “little group”. The little group with mixed string boundary conditions must preserve the endpoints as well as momentum. These boundary conditions break the full 26-dimensional Lorentz group, so the massive states above cannot be classified as representations of the unbroken 26-dimensional massive little group $SO(25)$. In a D1–D25 realization the common worldvolume has Lorentz group $SO(1,1)$ and trivial connected massive little group; $SO(24)$ acts on the ND coordinates as an internal rotation group. The scalar, vector and symmetric trace-free decomposition is under the surviving transverse SO(24), with this boundary-background qualification.

Use a time Laplace transform and a spatial Green function. The first-order transport term must be retained: the Airy resolvent kernel for ${\partial }_x^3+p$ alone would solve a different equation. The governing linear partial differential equation is the linear dispersive Stokes equation, distinct from the viscous-fluid Stokes equation.

For the time interval of interest, extend the given boundary signal by zero after $T$ and putThe possible discontinuity at $T$ does not affect a solution at $t<T$. Equivalently, use any sufficiently regular continuation beyond $T$: the resulting causal Green function produces the same earlier solution. We seek the usual solution with spatial decay and a time transform, rather than an unrestricted spatially growing solution on the half-line.

Writing $Q(x,p)$ for the time Laplace transform, the Laplace transform of a derivative givesThe corresponding characteristic roots are the three roots of a polynomialFor $\mathrm{Re}p>0$, label the unique root of negative real part by $r_{-}(p)$, and the two roots of positive real part by $r_1(p),r_2(p)$. This spatial root splitting for the dispersive Stokes resolvent is easily checked: if $r=ik$, then $p=-i(k-k^3)$ is imaginary, so no root crosses the imaginary axis as $p$ varies in the right half-plane. For real $p>0$ the cubic is strictly increasing on the real line and has one negative real root, while its conjugate pair has positive real parts. The root count remains the same throughout that connected half-plane. A repeated root would satisfy $3r^2+1=0$ and again put $p$ on the imaginary axis. All roots are therefore simple there, and $r_{-}(p)$ is a nonzero holomorphic function.

Set $d_j=3r_j^2+1$, including $j=-$. Define the Stokes resolvent kernel with advectionThe sign on the negative half-line is essential. To verify the Green function, use the partial fraction decompositionIts expansion at infinity yields the reciprocal-polynomial root identitiesConsequently $R_p$ and $R_p^{\prime }$ are continuous at zero, and $R_p^{\prime \prime }$ has jump one. Each piece solves the homogeneous ordinary differential equation, so the jump conditions for a third-order Green function imply, in the sense of distributional derivatives,Both pieces decay in their respective spatial directions.

Now defineIt satisfies the transformed inhomogeneous equation. The difference $Q-V$ solves the homogeneous linear ordinary differential equation; spatial decay leaves only $e^{r_{-}x}$. Imposing the Neumann boundary condition fixes its coefficient, giving the half-line Stokes solution with a boundary derivativeNo division by a vanishing root is involved, since $r_{-}=0$ would require $p=0$. Moreover the remaining boundary trace is explicitly known from $q_0$:Thus introduce the completely explicit initial-data kernelFor any upward Bromwich contour $\mathrm{Re}p=\sigma >0$ to the right of any transform singularities, the requested integral representation isThe contour integral is understood in the standard Bromwich inversion formula sense, with truncation limits if it is not absolutely convergent. Every term is expressed in the supplied data and algebraic characteristic roots; there are no unknown Dirichlet or second-derivative boundary traces.

The transformed equation verifies the linear partial differential equation and the initial term, while taking the spatial derivative gives $Q_x(0,p)=G(p)$. Bromwich inversion formula therefore supplies the prescribed derivative for $0<t<T$ and the initial value for $x>0$. For zero initial and boundary data, the transformed solution is a decaying homogeneous mode with zero boundary derivative, hence vanishes; this proves uniqueness in the stated transform/decay class. A change of the boundary continuation after $T$ has transform $e^{-pT}$ times another causal transform and contributes only after $T$, by the Laplace transform time-shift rule. The printed dot on $q_0$ is interpreted as its derivative in its sole variable $x$, so the stated corner compatibility is $q_0^{\prime }(0)=g_1(0)$. Additional differentiability at the corner requires the corresponding higher-order compatibility, but is not needed for the representation at positive $x,t$.

The drift makes the natural measure weighted. LetUnder homogeneous Neumann boundary conditions, the differential operator is self-adjoint in the weighted inner product $⟨u,v{⟩}_w={\int }_0^{L}w\overline{u}vdx$. Indeed integration by parts giveswhen the derivative boundary terms vanish. In particular its eigenvalues are nonpositive and it has the constant stationary eigenfunction. We first construct the weighted Neumann heat kernel with constant drift, then incorporate the two prescribed coordinate derivatives.

Put $b=\beta /2$. The Robin gauge transform for constant drift, $\varphi =e^{-bx}\psi$, turns $\mathcal{A}\varphi =-\lambda \varphi$ intoThese are coordinate-derivative Robin boundary conditions, with the same algebraic sign at both endpoints. The left condition for an oscillatory solution of wavenumber $k$ gives $\psi =\cos (kx)+(b/k)\sin (kx)$. The right condition then reduces to $(k+b^2/k)\sin (kL)=0$, giving $k_n=n\pi /L$. There is also a separate mode $\lambda =0$, $\psi =e^{bx}$, which must not be discarded by restricting to oscillatory $k$. For $0\le \lambda <b^2$, write $k=i\eta$: the endpoint determinant reduces to $\lambda \sinh (\eta L)=0$, leaving only that zero mode. At $\lambda =b^2$, the linear solution $\psi =A+Bx$ satisfies both endpoint conditions only if it vanishes, because $b>0$. The complete Sturm-Liouville eigenfunction expansion therefore usesTheir squared norms in the weighted measure areThe cancellation of the weight against $e^{-2bx}$ makes the latter formula an ordinary trigonometric integral. DefineFor $t>0$ this kernel is smooth in the interior. It integrates initial values against $w(y)dy$, not against unweighted Lebesgue measure. It is symmetric in $x,y$ as a weighted kernel, even though the drift operator does not look symmetric in an unweighted inner product.

To determine the boundary forcing, set $a_n(t)=⟨{\varphi }_n,q(\cdot ,t){⟩}_w/N_n$. Two applications of integration by parts leaveThe minus sign at the left endpoint is the sign of the boundary evaluation $[\cdot {]}_0^L$; the supplied $g_1$ is $q_x(0,t)$, not an outward normal derivative. Solving this linear differential equation by an integrating factor and recombining the modes gives the requested space-time integral representation:This is the Neumann boundary-forcing formula. The zero mode also gives the useful consistency checkOrdinary unweighted mass instead obeys $d{\int }_0^{L}qdx/dt=h_1-g_1+\beta [q(L,t)-q(0,t)]$. Replacing the weighted measure by $dy$ or deleting the constant mode would therefore give an incorrect solution.

There is a subtle point when checking nonzero derivative data: every homogeneous eigenfunction has zero endpoint derivative, but differentiating the boundary-forcing sum and its time integral term by term at an endpoint is not valid. The kernel is singular as $t-s\downarrow 0$, and the prescribed derivatives are interior limits of the full solution. An explicit Laplace transform form verifies those limits without this interchange and gives an alternative purely contour-integral version.

For $\mathrm{Re}p>0$, set $\kappa =\sqrt{p+b^2}$ using the principal square root, and defineBoth solve $(p-\mathcal{A})f=0$, with $u_p(0)=v_p(L)=1$, $u_p^{\prime }(0)=v_p^{\prime }(L)=0$. Their weighted Wronskian is the nonzero constantHence the resolvent kernel for Neumann advection-diffusion on an interval isThe first derivative in $x$ has jump $-e^{-\beta y}$, giving $(p-{\mathcal{A}}_x){\mathcal{R}}_p=e^{-\beta y}{\delta }_y$ in the sense of distributional derivatives. If $G(p)={\int }_0^T{e}^{-ps}{g}_1(s)ds$ and $H(p)={\int }_0^T{e}^{-ps}{h}_1(s)ds$, an equivalent transformed solution isIn fact $v_p^{\prime }(0)=-D_p$ and $u_p^{\prime }(L)=e^{-\beta L}{D}_p$, so $Q_x(0,p)=G(p)$ and $Q_x(L,p)=H(p)$ exactly. The inverse Bromwich contour integral of $e^{pt}Q(x,p)/(2\pi i)$ is an alternative final answer. Its poles at $p=0$ and $p=-{\lambda }_n$ reproduce the constant and decaying modes of $K_{\beta }$, respectively; apparent square-root singularities in the normalized hyperbolic factors are removable. This also independently checks the kernel normalization and boundary signs.

As in question 1, the printed dot on $q_0$ denotes its spatial derivative: the two corner conditions are $q_0^{\prime }(0)=g_1(0)$ and $q_0^{\prime }(L)=h_1(0)$. The initial eigenfunction expansion converges to $q_0$ in the weighted space and, for the stated smooth compatible data, has the appropriate classical initial and boundary limits. Differences of solutions with zero data satisfy $\frac{1}{2}d\Vert q{\Vert }_w^2/dt=-{\int }_0^{L}w|q_x{|}^{2}dx$, proving uniqueness in the usual regular solution class. In the limit $\beta \to 0$, $N_0\to L$, ${\varphi }_n\to \cos (n\pi x/L)$, and the ordinary Neumann heat kernel on an interval and its boundary signs are recovered.

For a scalar $\varphi$, the displayed basis identities giveTaking the divergence and differentiating the basis vectors yields the Laplacian in polar coordinatesConsequently,

The regular harmonic mode with angular dependence $\cos 2\theta$ is $C{r}^2\cos 2\theta$. Its radial derivative at $r=a$ is $2Ca\cos 2\theta$, so $C=1/(2a)$. A constant is invisible to the Neumann boundary condition, givingIf two solutions existed, their difference $u$ would satisfy ${\nabla }^{2}u=0$ and ${\partial }_{n}u=0$. Green's first identity gives ${\int }_D|\nabla u{|}^2={\int }_{\partial D}u{\partial }_{n}u-{\int }_{D}u{\nabla }^{2}u=0$, so $\nabla u=0$ and $u$ is constant. Thus these are all the solutions.

The positive uniform ellipticity condition from part (i) is needed here. It is not explicitly imposed in the printed preamble. Without it the conclusion is false: on $(0,\pi )$, take $L=u^{\prime \prime }$, $\beta =0$, and $f=0$. The Neumann boundary condition is satisfied by $u(x)=\cos (mx)$, and $(L+m^2)u=0$ for arbitrarily large shifts.

There is also an actual error in the PDF's supplied trace inequality: a nonzero constant has nonzero trace and zero gradient. Use instead the valid multiplicative trace inequality on a bounded smooth domainFor completeness, choose a smooth vector field $X$ agreeing with the outward normal vector on the boundary. The divergence theorem gives, for smooth $u$,Approximation and the Sobolev trace theorem extend it to $H^1(U)$. Thus the intended existence result remains true after correcting the hint.

Set $q=\Vert u{\Vert }_2$, $d=\Vert Du{\Vert }_2$, ${\beta }_{-}=\max \{-\beta ,0\}$, $c_{-}=\max \{-c,0\}$, andUniform ellipticity and the Young inequality giveFor $\lambda \ge C_0+1$, the bilinear form is coercive on the real Hilbert space $H^1(U)$ with coercivity constant $\min \{\theta /2,1\}$.

All coefficients are bounded, and the Sobolev trace theorem gives $\Vert Tu{\Vert }_2\le C_T\Vert u{\Vert }_{H^1}$. The Cauchy-Schwarz inequality applied to each volume and boundary term therefore proves $|B_{\lambda }[u,v]|\le C_{\lambda }\Vert u{\Vert }_{H^1}\Vert v{\Vert }_{H^1}$. Also $\ell (v)={\int }_{U}fv$ is a continuous linear functional, with $|\ell (v)|\le \Vert f{\Vert }_2\Vert v{\Vert }_{H^1}$. Every hypothesis of the Lax-Milgram theorem is now verified, and it gives