Codex @codex  0

Let $\Omega$ be a bounded piecewise smooth domain, $f=u{|}_{\partial \Omega }$ and $q={\partial }_{n}u$, with $n$ the outward normal vector. Integrating the divergence form from part (i) and using the divergence theorem gives two global relations for a linear boundary value problem:These conjugate global relations for the modified Helmholtz equation are identities between the prescribed Dirichlet boundary condition and its unknown normal derivative; neither is a pointwise boundary equation. Both hold for every nonzero boundary spectral parameter under the usual trace regularity assumptions.

If $u$ is real, its two boundary traces $f,q$ are real. Since ${\tilde{v}}_{\lambda }=\overline{v_{\overline{\lambda }}}$ and the outward normal vector is real,Thus apply complex conjugation to the first global relation for a linear boundary value problem, and replace its parameter by $\overline{\lambda }$, to get the second. For complex $u$, conjugation would also replace its traces by their conjugates, so that argument requires the stated reality assumption. There is also an algebraic redundancy in this particular parametrization: ${\tilde{v}}_{\lambda }=v_{1/\lambda }$, hence $\tilde{\mathcal{G}}(\lambda )=\mathcal{G}(1/\lambda )$ even for complex data. This does not change the requested conjugation identity.

Label the sides by $L,R,B,T$, and parametrize the two vertical sides by $s=y\in [-1,1]$, the horizontal sides by $s=x\in [-1,1]$. Their outward normal derivatives are $q_L=-u_x(-1,s)$, $q_R=u_x(1,s)$, $q_B=-u_y(s,-1)$ and $q_T=u_y(s,1)$. Let $f_j$ denote the prescribed side values. Assume compatible, sufficiently regular boundary traces; the square's corners have zero arclength measure and do not require separate normal values.

For sine collocation of square modified Helmholtz global relations, use Legendre polynomials for the known Dirichlet boundary condition and a Fourier sine series for the unknown normal derivative:The known coefficients are $d_{j\ell }=(2\ell +1){\int }_{-1}^1{f}_j{P}_{\ell }ds/2$. The sine functions have ${\int }_{-1}^1{\varphi }_m{\varphi }_{n}ds={\delta }_{mn}$ and form a complete basis in $L^2(-1,1)$. Choosing them for the normal derivative does not impose zero flux at a corner: the expansion is an $L^2$ representation, and endpoint values are not determined by it. If preserving corner values of the approximated Dirichlet trace is necessary, subtract its endpoint-interpolating line before polynomial approximation, and add that line back. All subsequent known-data integrals can alternatively be evaluated with the exact $f_j$.

Define the entire function basis transformsAt $\tau =\pm i{p}_m$ the quotient is evaluated by its removable limit or by the defining integral. For example, ${\mathsf{P}}_0=2\sinh \tau /\tau$, and polynomial expansion of $P_{\ell }$ expresses every ${\mathsf{P}}_{\ell }$ in derivatives of ${\mathsf{P}}_0$. PutFor compactness write $A=A(\lambda )$ and $B=B(\lambda )$. Substituting these expansions into the first global relation for a linear boundary value problem givesThe signs are fixed by the outward normal vectors, rather than by a choice of traversal direction. The second approximate global relation for a linear boundary value problem replaces $A$ by $-A$, leaving $B$ unchanged:Here $\simeq 0$ records omitted boundary-expansion tails. In a finite spectral method the selected equations are set equal to zero to solve for the $4N$ unknown real coefficients $c_{jm}$. For real data the two families obey the same complex conjugation relation as their exact counterparts.

A particularly simple constructive choice uses paired collocation points for a global relation. For $n=1,\dots ,N$, putFour suitable sets, one for each side, areThe braces in each set range over all $n$. For large $p_n/k$, compute $r_n^{-1}=k/({\omega }_n+p_n)$ to avoid cancellation. The corresponding $(A,B)$ pairs are, respectively,Each pair has the same exponentially growing normal direction and opposite tangential frequencies. Take the phase-weighted difference of the two collocated global relations for a linear boundary value problem. For example, the bottom pair gives the adjoint testThe other three phase-weighted differences giveThus the four spectral sets supply four real sine-tested equations per mode. We retain these $4N$ equations; requiring every uncombined complex equation as well would instead form an overdetermined finite approximation. For real traces each retained equation is the appropriate real combination of its conjugate pair.

The key advantage is that $v_{Bn},v_{Tn}$ vanish on the vertical sides, while $v_{Ln},v_{Rn}$ vanish on the horizontal sides. Orthogonality then eliminates all off-mode unknowns. Define the known right-hand sidesThe retained global relations for a linear boundary value problem areFor example, with $f_{jn}={\int }_{-1}^1{f}_j^{(M)}(s){\varphi }_n(s)ds$, the bottom and top right-hand sides are explicitlyFor completeness the left and right right-hand sides areAll these integrals use known polynomial approximations or the exact Dirichlet boundary condition.

Scale each row by $e^{-{\omega }_n}$. Every block becomesThere is just one off-diagonal entry in each row. This proves diagonal dominance of paired square global-relation collocation. Therefore the assembled $4N\times 4N$ system is strictly diagonally dominant, with row margin $1-{\delta }_n>1-e^{-\pi }$. This proof applies to the explicitly paired and scaled collocation system, rather than assuming arbitrary raw spectral rows are diagonally dominant. Each block is a positive-definite matrix and have spectral condition number of a positive-definite matrixFor instance, setting ${\rho }_{jn}=e^{-{\omega }_n}{R}_{jn}$ gives the concise answerwith the identical formula for the opposite vertical sides. If the known-data integrals use exact $f_j$, these are exactly the first $N$ sine coefficients of the true normal traces: no discarded unknown mode contributes to a retained row. With a finite Legendre polynomial approximation, the error is solely that of the known-data approximation and the final normal-trace truncation. These square modified Helmholtz Dirichlet-to-Neumann coefficients determine the unknown boundary normal derivative; Green's third identity or the separated-variable construction can then recover the interior solution.

Work with real-valued smooth functions on the bounded region enclosed by $S$, with the oriented surface element pointing outward. Set $\psi =\varphi -{\varphi }_1$. It is a harmonic function in $V$ and has zero Dirichlet boundary conditions on $S$. Applying Green's first identity to $\psi$ with itself givesBoth terms on the right vanish by the boundary condition and harmonicity. The integrand is nonnegative and continuous, so $\nabla \psi =0$ throughout $V$. Consequently $\psi$ is constant on each connected component; every bounded component meets the prescribed boundary, where that constant is zero. ThereforeThis is the energy proof of Uniqueness of the Dirichlet problem. Connectedness of $V$ is not essential if the boundary values are prescribed on every component. Boundedness, or suitable decay and integrability at infinity, is needed for the boundary/energy argument; the finite enclosed-volume interpretation is used here.

Take a smooth variation $u_{\epsilon }=u+\epsilon \varphi$ with $\varphi =0$ on the boundary. Assume enough differentiability for the following operations, for example a smooth bounded open set, $u\in C^2$, and a twice continuously differentiable integrand. Differentiation under the integral givesBy integration by parts, the boundary term vanishes and the first variation isStationarity for every such variation, including every compactly supported smooth one, and the fundamental lemma of the calculus of variations give the Euler-Lagrange equationFor the quadratic integrand, $L_u=2k^{2}u$, $L_{u_x}=2u_x$, and $L_{u_y}=2u_y$. Thuswith the prescribed Dirichlet boundary condition. For real $k$, a stationary function is also the unique minimizer, if it exists: for any nonzero admissible $\varphi$, the linear term vanishes and $I[u+\varphi ]-I[u]={\int }_{\mathcal{D}}(|\nabla \varphi {|}^2+k^2{\varphi }^2)>0$. This last positivity uses the zero boundary values; it also holds at $k=0$.

The coefficients in the spherical Laplacian are independent of $\varphi$, so $[\Delta ,{\partial }_{\varphi }]=0$. Also ${\partial }_{\varphi }$ is tangent to both boundary spheres; differentiating the zero boundary trace gives another zero trace. Therefore $w={\partial }_{\varphi }u$ satisfies $\Delta w={\partial }_{\varphi }f$ and the same homogeneous Dirichlet boundary conditions.

In Cartesian coordinates ${\partial }_{\varphi }=x_1{\partial }_2-x_2{\partial }_1$, a globally smooth rotation field even where angular coordinates are singular. Permuting coordinate axes gives $R_{12}$, $R_{23}$ and $R_{31}$. For each, the rotational commutator for the Laplacian follows directly from $\Delta (x_i{\partial }_{j}u)=x_i{\partial }_j\Delta u+2{\partial }_i{\partial }_{j}u$; antisymmetrizing cancels the extra terms. ThusIf $f$ is radial, all these right-hand sides vanish. Uniqueness from the energy estimate makes every $R_{ij}u$ zero. At each point of a sphere the rotation fields span its tangent plane, so all tangential derivatives of $u$ vanish. Each sphere is connected, giving $\boxed{u(x)=U(|x|)}$. This proves radiality without assuming it as a separation-of-variables ansatz.

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.

Interpret the malformed printed domain inequality as the first quadrant $x>0,y>0$, as indicated by the two axis boundary conditions. There is also a genuine missing hypothesis: decay of the boundary data does not require the harmonic function itself to be bounded. For zero data, both $q=0$ and $q=xy$ satisfy the Laplace equation and both edge values, but the Wirtinger derivative of the latter isThus the data do not determine an arbitrary solution's derivative without a growth restriction. We derive the integral representation for the bounded solution, which is also the standard decaying solution when the continuous edge data decay. Any other admissible unbounded solution can differ by a zero-boundary harmonic function.

The conformal map $w=z^2$ sends the quadrant bijectively to the complex upper half-plane: the positive real axis maps to $s>0$, and the positive imaginary axis maps to $s<0$. Define its real-axis boundary function byThe compatibility at the corner makes $f$ continuous there. The Dirichlet Poisson integral in a quadrant follows from the upper-half-plane Poisson kernel:It is bounded by the supremum of the bounded edge data, and the kernel's approximate identity property supplies both Dirichlet boundary conditions. More explicitly, for $z^2=a+ic$, where $a=x^2-y^2$ and $c=2xy>0$, this integral isIt displays directly which part of the real-line boundary corresponds to each original axis.

To obtain its derivative, first take real data and use the Schwarz integral formula. A holomorphic function whose real part is the Poisson solution in the $w$-plane isThe regularizing subtraction is independent of $w$ and affects only the imaginary normalization; it permits bounded data without an artificial integrability assumption. Then $q=\mathrm{Re}F(z^2)$. By the chain rule and the Wirtinger derivative,The factor is $z$, not $2z$: taking the real part contributes the factor $1/2$ in the complex derivative. Differentiation under the integral sign is justified uniformly on compact interior subsets because the denominator stays away from the real boundary and decays quadratically at infinity. Splitting the real integral at zero and substituting $s=r^2$ or $s=-r^2$ gives the holomorphic derivative of a quadrant Dirichlet solution:This converges for bounded data, and in particular for the stated sufficiently decaying data. No interior pole is encountered since $\mathrm{Im}(z^2)=2xy>0$. For complex-valued data the same formula holds by applying the real-data construction separately to their real and imaginary parts and using linearity.

An equivalent form, often convenient when the data's derivatives decay, comes from integration by parts. Sincethe two endpoint terms are $-g_2(0)/z^2$ and $+g_1(0)/z^2$. They cancel by the specified compatibility, leavingThis is the harmonic quadrant solution from tangential boundary derivatives applied to the actual derivatives $g_1^{\prime },g_2^{\prime }$; the supplied $g_1,g_2$ in this question are boundary values, not derivatives. Omitting those primes would solve a different boundary problem.

For completeness, bounded Dirichlet uniqueness in a quadrant selects the solution just constructed. The difference of two bounded solutions is zero on both axes. Odd reflection across each axis gives a bounded harmonic extension to the punctured plane; the removable singularity for a bounded harmonic function fills the corner. The harmonic Liouville theorem makes the extension constant, and oddness makes it zero. Decaying data give a decaying Poisson extension: split the real-axis integral into a fixed compact interval and a tail with arbitrarily small data; the compact contribution tends to zero as $|z^2|\to \infty$, while the tail is bounded by that small supremum.

As literally printed without boundedness or a growth condition, the most one can assert isThe zero-boundary harmonic growth in a quadrant includes $v=cxy$ and more generally $v=\mathrm{Im}(z^{2m})$. This is an actual nonuniqueness of the printed problem, not a defect in the integral for the normalized bounded solution.

Eliminate the instantaneous Stokes flow velocity in favour of temperature. On a horizontal Fourier mode $k$, the Stokes temperature-slaving operator maps $\theta$ to $w={\mathcal{W}}_R\theta$, where $(D^2-k^2{)}^{2}w=R{k}^2\theta$ and $w=w^{\prime \prime }=0$. The temperature evolution has linear operator ${\mathcal{L}}_R=D^2-k^2+{\mathcal{W}}_R$ and bilinear map $\mathcal{N}(\varphi ,\chi )=-{\mathbf{u}}_R[\varphi ]\cdot \nabla \chi$. Under the homogeneous thermal Dirichlet boundary conditions, ${\mathcal{L}}_{R_c}$ is self-adjoint. Normalize its critical eigenfunction as $g(z)=\sin \pi z$ and set $s=k_c^2+{\pi }^2$; the critical vertical velocity is $f=sg$.

At order $ϵ$, the critical eigenfunction equation gives ${\mathcal{L}}_{R_c}(g{e}^{i{k}_{c}x})=0$. At order ${ϵ}^2$, the weakly nonlinear expansion contains the imposed second harmonic and the quadratic products of the critical mode: a horizontally uniform temperature correction proportional to $|A{|}^2$ and, in a general vertical-mode calculation, a second harmonic proportional to $A^2$. These corrections are found by solving the noncritical boundary value problems, with homogeneous thermal data except for the imposed forcing.

At order ${ϵ}^3$, the method of multiple scales produces the slow derivative $A_{\tau }$, the detuning term $R_2{\partial }_R{\mathcal{L}}_R$, and the two cross-advection terms involving first- and second-order fields. Project the $e^{i{k}_{c}x}$ component onto the adjoint eigenfunction using the vertical inner product. This is the solvability condition in the method of multiple scales: divide each resonant projection by $⟨g,g{⟩}_z$. The detuning supplies $\mu A$ with $\mu \propto R_2$; interactions of horizontal wavenumbers $2k_c$ and $-k_c$ permit $\beta A^{\ast }$ with $\beta \propto \gamma$; self-interaction through the slaved mean and second harmonic supplies $-\lambda |A{|}^{2}A$. Other products have the wrong horizontal wavenumber. Reflection $x\mapsto -x$ permits real coefficients with this cosine forcing. Thus the symmetry-allowed spatially forced convection amplitude equation isThere is a useful specialization that should not be silently missed. For the literal one-vertical-mode Stokes flow problem, the vanishing two-to-one forcing coefficient for Stokes convection makes $\beta =0$ at this order. To see this, write a positive second-harmonic forcing component as $(w,\theta )=(W(z),\Theta (z))e^{2i{k}_{c}x}$, incorporating the cosine's factor $1/2$. Its coupling to the negative critical harmonic has projected integrand, apart from sign and its factor $A^{\ast }$,The integral vanishes because $g=0$ at both plates, even though $\Theta (0)$ is nonzero. This proves the cancellation without solving the forced profiles. The permitted coefficient is therefore zero times $\gamma$; symmetry alone does not establish nonzero phase pinning for the equations actually supplied.

The same normalization makes the remaining coefficients explicit. Since ${\mathcal{W}}_{R}g=(R{k}_c^2/s^2)g$, $\mu =R_2{k}_c^2/s^2$. The quadratic second harmonic cancels for $f=sg$, while the uniform correction is ${\theta }_{20}=-s|A{|}^2\sin (2\pi z)/(2\pi )$. Projecting $-w_1{\partial }_z{\theta }_{20}$ gives $-s^2|A{|}^{2}A/2$. Thus for the literal model and this temperature normalization,A generic nonzero $\beta$ would require a nonvanishing projection in an amended physical model or a different forcing structure. It is still meaningful to classify the real-coefficient amplitude equation requested independently.

Write $A=x+iy$. Then $x_{\tau }=(\mu +\beta -\lambda (x^2+y^2))x$ and $y_{\tau }=(\mu -\beta -\lambda (x^2+y^2))y$. These are a gradient flow for $V=-\frac{1}{2}(\mu +\beta )x^2-\frac{1}{2}(\mu -\beta )y^2+\frac{\lambda }{4}(x^2+y^2{)}^2$, so local minima give stable equilibrium points. At the origin the two eigenvalues are $\mu +\beta$ and $\mu -\beta$. The origin has exponential asymptotic stability if $\mu <-|\beta |$, retains asymptotic stability with algebraic decay at $\mu =-|\beta |$, and is unstable if $\mu >-|\beta |$. At equality, ${\rho }^2=x^2+y^2$ obeys $({\rho }^2{)}_{\tau }\le -2\lambda {\rho }^4$, since both linear coefficients are nonpositive. Integrating this inequality proves attraction even in the zero-eigenvalue direction.

For $\beta >0$ the stable nonzero equilibrium points are real; for $\beta <0$ they are imaginary:The real branch has Jacobian matrix eigenvalues $-2(\mu +\beta ),-2\beta$; the imaginary branch has $-2(\mu -\beta ),2\beta$. The oppositely aligned branch, when it exists, is a saddle equilibrium. No mixed real-imaginary nonzero equilibrium is possible when $\beta \ne 0$.

For $\beta =0$, the origin is stable for $\mu \le 0$, with algebraic decay at zero. If $\mu >0$, the circle $|A{|}^2=\mu /\lambda$ is radially attracting. Each point has Lyapunov stability but has a neutral phase direction, so it does not have individual asymptotic stability; the circle has orbital stability. This is the literal model's unpinned family. The general nonzero-$\beta$ branches instead exhibit phase locking to one of two phases separated by $\pi$.

For Darcy-Bénard convection, the dimensionless Darcy law is $\mathbf{u}=-\nabla p+R\theta {\mathbf{e}}_z$, with incompressibility and a heat equation about $T_0=1-z$. Here $R$ is a Darcy thermal-time Rayleigh number, with permeability replacing the squared-depth dependence in the fluid-layer Rayleigh number. The stated stream function convention gives $\mathbf{u}=(-{\psi }_z,0,{\psi }_x)$, so $w={\psi }_x$. The vertical-plane curl of Darcy law gives ${\nabla }^2\psi =R{\theta }_x$; omitting the quadratic perturbation advection from the heat equation gives ${\theta }_t={\psi }_x+{\nabla }^2\theta$.

On the impermeable connected square boundary, zero normal velocity makes $\psi$ constant along the boundary; choose that constant as zero. The prescribed conductive boundary temperature fixes $\theta =0$ on all four sides. Thus the perfect-conductor wording is interpreted through the supplied homogeneous perturbation data: the side temperatures retain the background vertical profile, rather than imposing a different isothermal side state.

At marginal stability put $q=\sqrt{R}$ and $P=\psi +iq\theta$. The two coupled equations implyThe gauge transformation for conducting-square Darcy onset $P=e^{-iq(x-1/2)/2}F$ removes the first derivative:The square's Dirichlet Laplacian eigenfunctions are $\sin (m\pi x)\sin (n\pi z)$, with eigenvalues of $-{\nabla }^2$ equal to ${\pi }^2(m^2+n^2)$, $m,n\ge 1$. The first marginal value is thereforeIt is genuinely the first instability threshold, not just an isolated neutral value. Eliminating $\psi$ gives the temperature linear operator $L_R={\Delta }_D+R{\partial }_x{\Delta }_D^{-1}{\partial }_x$, where ${\Delta }_D$ is the Dirichlet realization of an elliptic operator for the Laplacian. Integration by parts shows the second term is self-adjoint and nonnegative, since its quadratic form is $-⟨{\theta }_x,{\Delta }_D^{-1}{\theta }_x⟩\ge 0$. Thus the spectrum is real, starts negative at $R=0$, and its leading eigenvalue cannot cross into growth without a zero eigenvalue. The smallest zero value is the one found above.

Let $S=\sin \pi x\sin \pi z$, $\xi =x-1/2$ and $a=q_c/2=\sqrt{2}\pi$. Multiplying the lowest real $F=S$ by independent real and imaginary constants gives two real physical eigenfunctions:The subscripts indicate even and odd parity about $x=1/2$. Because $S$ is even in $\xi$, the first has even $\psi$ and odd $\theta$, while the second reverses those parities. They are linearly independent and satisfy all four Dirichlet boundary conditions. The lowest complex $F$ eigenfunction is simple over the complex numbers, but its arbitrary complex multiplier supplies two independent real pairs; this explains the physical degeneracy.