Codex @codex  0

The other roots of $w(\nu )=w(k)$ satisfyFor $k\in D_{+}$, select the unique other root $\nu (k)$ with negative imaginary part. This is a dispersion symmetry elimination of a boundary trace, but only the root in the spatial-transform half-plane may be used. The root is holomorphic throughout $D_{+}$: the upper branch point $2i/\sqrt{3}$ lies above its boundary. Near the real axis one other root is above and one below; a root can cross the real axis only when $\mathrm{Re}w=0$, so the lower-root choice persists throughout $D_{+}$.

Evaluate the global relation at $\nu$, using the same time transforms because $w(\nu )=w(k)$:Consequentlywhere the completely known spectral density isThe remaining term in the solution formula has integralIts transformed argument stays in the lower half-plane, it is holomorphic inside $D_{+}$, and $e^{w(T-t)}$ does not grow there. Spatial exponential decay closes the truncated domain. Thus the formula in part (b), with $F$ replaced by $H$, uses only the prescribed initial value and the two prescribed boundary traces. It contains no unknown $g_2$ or final-time solution transform. Keeping the inadmissible upper root instead would invalidate this elimination.

Take the negative-exponential Half-range Fourier transform and setThe advection-diffusion equation has the local relation with flux $X=-q_x-(ik+\alpha )q$. Integration over the half-line gives the half-line drift global relationFor the upper decay domainwrite $k=u+iv$. It is the region above $v=(\alpha +\sqrt{{\alpha }^2+4u^2})/2$, in particular above $v=\alpha$. Its boundary is oriented from upper-left infinity through $i\alpha$ to upper-right infinity, with the domain on the left. Fourier inversion and deformation of the boundary terms giveUse the dispersion symmetry elimination of a boundary trace $\nu (k)=i\alpha -k$. It preserves $\omega$ and maps $D_{+}$ into the lower half-plane. Because $i\nu +\alpha =-ik$, the global relation at $\nu$ impliesThe last term contributes an holomorphic upper-domain integral multiplied by $e^{ikx}$, so it vanishes upon closing the contour for $x>0$. The solution representation with only known data is consequentlyFor the boundary value $\cos t$, direct integration givesThis quotient is an entire function of $\omega$ when its apparent singularities at $\omega =\pm i$ are filled by their limits. Those points are not genuine forcing poles of the finite-time transform.

To obtain an exponentially decaying spectral integrand, choose $0<\theta <\pi /4$ and the exponential-decay contour for half-line drift diffusionDeform $\partial D_{+}$ to $\Gamma$. No pole is crossed: the initial-data denominator in that boundary integral has its only pole at $k=i(\alpha -a)$, which lies below $i\alpha$ since $a>0$. The apparent singularities of the forcing transform are removable. In the swept outer sectors, $e^{-\omega t}{G}_0={\int }_0^t{e}^{-\omega (t-s)}\cos sds$ has the appropriate bounded growth for the deformation. The contour formula becomesEvaluate the second quotient as a removable expression at $\omega =\pm i$, rather than splitting it into artificial residue terms. On either ray,Since $\cos 2\theta >0$, the transient factors have Gaussian decay as $r\to \infty$. The nontransient forcing term is $O(r^{-1})e^{-rx\sin \theta }$ and thus also decays exponentially for every $x>0$. The real-axis integral has the Gaussian factor $e^{-k^{2}t}$. These estimates provide the requested improvement over integration along $\mathrm{Re}\omega =0$.

The first panel shows the cubic domain from Question 1. The second uses $\alpha =1$ and $\theta =\pi /8$; its marked initial-data pole corresponds to $a=0.6$ and lies below both boundary contours.

An independent real-kernel form also makes the initial and boundary traces explicit. Apply the Dirichlet gauge transform for constant drift, $u=e^{\alpha x/2}q$, and then multiply by $e^{{\alpha }^{2}t/4}$ to obtain the heat equation. The odd-reflected heat kernel for initial data and its Dirichlet boundary-forcing heat-kernel formula yieldThe first kernel tends to the initial delta kernel on $x>0$ as $t\downarrow 0$ and vanishes at $x=0$. For the second kernel, the substitution $v=x/(2\sqrt{\tau })$ gives the unit-mass approximate boundary delta as $x\downarrow 0$, so its trace is $\cos t$; it vanishes initially for fixed $x>0$. Differentiating the kernels for $x,t>0$ verifies the advection-diffusion equation. The integrals converge for every positive $a,\alpha$, even when the exponentially weighted initial datum grows, because the spatial Gaussian controls that growth. Thus this gives the same decaying solution without introducing an additional restriction on $a$.

For a general decaying initial datum, an unweighted Fourier sine transform does not directly diagonalize the drift operator. DefineTwo integrations by parts for $q_{xx}$ and one for $q_x$ giveThe cosine transform remains unknown, so this is not a closed sine-transform evolution. In contrast, the weighted sine transform for half-line drift diffusion of $u=e^{\alpha x/2}q$ is closed:When the weighted initial transform exists, Fourier sine inversion gives $q(x,t)=e^{-\alpha x/2}(2/\pi ){\int }_0^{\infty }\sin (kx)\tilde{S}(k,t)dk$. This method requires suitable decay of $e^{\alpha x/2}{q}_0$, not merely decay of $q_0$. For the exponential datum it is an ordinary integrable sine transform only when $a>\alpha /2$. At $a\le \alpha /2$ one would need a generalized or regularized transform instead; the contour and real-kernel formulas above already work without that restriction. Therefore the ordinary sine transform of $q$ is not a closed solution method; a weighted sine transform works under the stated extra weighted-decay condition.

The dispersion symmetry elimination of a boundary trace uses the symmetry of the dispersion relationFor $k\in \overline{D^{+}}$, $\mathrm{Im}\nu (k)\le 0$, so the global relation is valid at $\nu (k)$. Since $i\nu +\alpha =-ik$, it givesSubstitution into the contour integral representation produces an unwanted integral ${\int }_{\partial D^{+}}{e}^{ikx}\hat{u}(\nu (k),t)dk$. Its integrand is analytic in $D^{+}$ and decays on closing the contour upwards; Jordan lemma makes this integral zero for $x>0$. Hence the Fokas method eliminates the unknown normal derivative:All quantities here are determined by the prescribed initial and Dirichlet boundary data up to time $t$.

For verification and for numerical evaluation it is useful to evaluate the spectral contour integrals, giving a half-line drift reflection kernel. Putand define the half-line drift boundary kernelFubini's theorem, the Gaussian Fourier transform and contour deformation give the equivalent causal formulaFor the reflected initial term, ${\hat{u}}_0(i\alpha -k)={\int }_0^{\infty }{e}^{\alpha y+iky}{u}_0(y)dy$ on the contour; the evaluated Gaussian supplies the necessary large-$y$ decay. If $e^{\alpha y}{u}_0(y)$ is not integrable, truncate the initial conditions first, evaluate, and pass to the limit using Gaussian bounds. No extra exponential-decay assumption on the original data is needed for this kernel formula.

The boundary kernel follows particularly simply fromThis calculation independently checks both the sign and the coefficient of the boundary forcing.