Codex @codex  0

The heat equation is linear and homogeneous in $u$: the operator $L={\partial }_t-\lambda {\partial }_x^2$ satisfies $L(A{u}_1+B{u}_2)=AL{u}_1+BL{u}_2$. Thus the superposition principle applies to its solutions. Spatial translations are also allowed because the coefficients are constant.

Choose $a=1/(4\lambda )$ in part (b) and multiply by $\sqrt{a}$ to make the initial amplitude one. Taking the two shifts $b=1$ and $b=-1$ then givesAt $t=0$ this reproduces both initial peaks. Each term is an instance of Gaussian initial data for the heat equation, so their sum satisfies the equation. Its extension to negative times is a real smooth solution precisely on the connected intervalAt the lower endpoint the peaks narrow without bound and their heights diverge; there is no finite classical continuation there. For positive times each peak spreads and decreases in height while preserving its integral.

For the Gaussian initial data for the heat equation, substitute the initial profile into the heat-kernel solution. Completing the square in $y$ givesThe remaining Gaussian integral yieldsThe centre stays at $x_0$ and the squared width grows linearly with time.

For the interval pulse, the convolution reduces to $B{\int }_a^{b}G(x-y,t)dy$. Using $u=(y-x)/(2\sqrt{t})$ and the error function givesThis is the translated and scaled heat equation with interval-indicator initial data. At its jump points the $t\downarrow 0$ limits are $B/2$; away from them the initial data are recovered pointwise.

For fixed $a<x<b$ at late time, both error function arguments are small. Since $\mathrm{erf}u=2u/\sqrt{\pi }+O(u^3)$,The estimate is uniform for $x\in [a,b]$. Equivalently, on this fixed interval the Gaussian in the convolution is $1+O((b-a{)}^2/t)$.

The spreading distance is more apparent on the growing spatial scale. For the Gaussian pulse the exponent is order one when $|x-x_0|=O(\sqrt{t})$. For the interval pulse, set $x=c+\xi \sqrt{t}$ with $c=(a+b)/2$ and fixed $\xi$. Uniformly for bounded $\xi$, its convolution givesThus both nonzero pulses spread over a distance $O(\sqrt{t})$, while their heights are $O(t^{-1/2})$ and their total integrals remain constant. More quantitatively, the variances of their normalized profiles are $1/(2a_0)+2t$ and $(b-a{)}^2/12+2t$, respectively, confirming the same width law.

Complete the square in the Fourier transform exponent:The supplied shifted Gaussian integral therefore givesThis is the Fourier transform of a Gaussian, with angular-wavenumber convention $e^{-ikx}$.

Take the spatial Fourier transform of the heat equation. Two integrations by parts, with decaying fields and first derivatives at infinity, give $\mathcal{F}[{\psi }_{xx}]=-k^2\tilde{\psi }$. Thus each mode satisfies the ordinary differential equationInserting the transformed initial condition and integrating the time equation givesFinally set $A=t+1/(4\lambda )>0$. The inverse Fourier transform, evaluated by the same Gaussian identity, givesEquivalently the Gaussian initial data for the heat equation evolve asThe profile broadens and its peak decreases, while its total integral remains $\sqrt{\pi /\lambda }$, as required by diffusion without loss on the whole line. Direct differentiation verifies the equation and the initial limit.