Codex @codex  0

For ${\eta }_t+\mathbf{c}\cdot \nabla \eta =\mu \eta +\nabla \cdot (\Gamma \nabla \eta )$ with symmetric positive definite $\Gamma$, the dispersion relation is $\omega =\mathbf{c}\cdot \mathbf{k}+i(\mu -{\mathbf{k}}^T\Gamma \mathbf{k})$. Its zero-group velocity saddle is ${\mathbf{k}}_0=-i{\Gamma }^{-1}\mathbf{c}/2$. The anisotropic heat kernel has fixed-position exponential rate $\mu -{\mathbf{c}}^T{\Gamma }^{-1}\mathbf{c}/4$, so this accessible saddle gives the actual threshold. At fixed speed $q$, the threshold ranges from $q^2/(4d_{\max })$ to $q^2/(4d_{\min })$, where $d_{\min },d_{\max }$ are the eigenvalues of $\Gamma$. Absence of absolute hydrodynamic instability in every direction requires the upper bound $\mu \le q^2/(4d_{\max })$.

The transition density of standard Brownian motion in ${\mathbb{R}}^d$ is the heat kernelThus ${\mathbb{P}}_x(X_t\in A)={\int }_A{p}_t(x,y)dy$.

For the negative-flux viscous Burgers equation, the Cole-Hopf transformation uses $f=2\alpha {\partial }_{\theta }\log \psi$. Initial N-wave data give ${\psi }_0=\exp [U(L^2-{\theta }^2)/(4\alpha )]$ inside $[-L,L]$ and one outside. Completing the square in the heat kernel convolution gives, with $a=1+UZ$,The Gaussian interval masses and their exponentially weighted tails must be treated together when taking the small-diffusion limit.

Consider the one-dimensional diffusion equation on the whole real line,with initial value one for $x<0$ and zero for $x>0$. This is an idealized Heaviside step function profile, not a measured material or an apparatus specification. The heat kernel gives a convolution with the initial profile:The change of variables $s=(x-y)/(2\sqrt{Dt})$ turns the convolution into a Gaussian integral tail. The resulting similarity variable $\eta =x/(2\sqrt{Dt})$ shows that the transition width grows like $\sqrt{Dt}$: multiplying elapsed time by four doubles every fixed-profile distance from the original interface.

At every positive time, $u(0,t)=1/2$, while $u(x,t)\to 1$ as $x\to -\infty$ and $u(x,t)\to 0$ as $x\to +\infty$. For every fixed $x\ne 0$, the limit as $t\to 0^{+}$ returns the specified initial value. The value assigned at the single initial discontinuity does not change this convolution. The error function therefore explains smoothing without pretending that a discontinuous initial profile remains sharp.

Direct differentiation checks the partial differential equation. Writing $\eta =x/(2\sqrt{Dt})$ givesThe model assumes constant diffusion coefficient, no advection, and no finite-domain boundary effects. For an interpretation as heat conduction, the heat equation applies to a normalized temperature difference under its own material assumptions.

For standard Brownian motion in ${\mathbb{R}}^d$, $p_t(y-x)$ is the transition density from $x$ to $y$. It is the heat kernel for generator $\frac{1}{2}\Delta$. For every unit vector $e$, $\int |{\partial }_e{p}_t(x)|dx=\sqrt{2/\pi }/\sqrt{t}$, which bounds the change in this density under a spatial translation and proves the Gaussian heat-kernel proof of the harmonic Liouville theorem.

On ${\mathbb{R}}^d$, the heat equation with initial Gaussian function $u_0(x)=e^{-|x-x_0{|}^2/a}$, $a>0$, has the displayed solution. Convolving the initial profile with the heat kernel adds the Gaussian width parameters; conservation of its spatial integral fixes the prefactor. Direct differentiation provides another proof. The width increases while the peak decreases, with no loss of the total integral.

The Gaussian integral gives the exact area under a Gaussian function. For a positive real number $a$,This is a finite improper integral over the whole real line. For example, $I(1)=\sqrt{\pi }\approx 1.772454$ and $I(4)=\sqrt{\pi }/2\approx 0.886227$. Increasing $a$ from one to four narrows the curve and halves its area, while its height at zero stays one. The change of variables $u=\sqrt{a}x$ explains the factor $a^{-1/2}$.

The useful consequence is an exact normalization. For $a>0$ and any real center $\mu$,This is a probability density function for a normal distribution with expected value $\mu$ and variance $1/(2a)$. A narrower probability density function must be taller to keep total probability one. In the usual notation $a=1/(2{\sigma }^2)$, the prefactor is $1/(\sigma \sqrt{2\pi })$, where $\sigma$ is the standard deviation.

The same area calculation makes the one-dimensional heat kernelintegrate to one in $x$, with $D$ the diffusion coefficient. A spreading Gaussian function can therefore preserve total mass while its height falls. This links the Gaussian integral with both the normal distribution and the diffusion equation.

Finite intervals lead to the error function, rather than this single constant:The complementary tail is described by the complementary error function. These definitions and their real-axis limits are given in the NIST reference. ^[ref]

The polar-coordinate proof of the Gaussian integral explains why $\pi$ appears in an integral of a function defined on a line: squaring its area creates a rotationally symmetric integral on a plane. Substitution $t=x^2$ also gives $\Gamma (1/2)=2{\int }_0^{\infty }{e}^{-x^2}dx=\sqrt{\pi }$, connecting the result with the gamma function.

The normalized Gaussian mass of an interval isIt is the convolution of an interval indicator with the heat kernel, and has an error function representation. As the variance tends to zero it converges to the interval indicator away from its endpoints, with value one half at an endpoint. This is a concrete approximate identity.

For quadratic free energy with positive $a,\kappa$ and kinetic coefficient $\Gamma$, a deterministic initial Dirac delta distribution evolves under nonconserved order-parameter dynamics to the meanTaking the expectation removes the centered Gaussian white noise; the resulting diffusion-reaction equation is solved by the heat kernel. The integrated mean is $A{e}^{-\Gamma at}$, so spreading does not conserve the initial excess. The connected Fourier mode covariance grows from zero to $k_{B}T/(a+\kappa q^2)$ as $1-e^{-2\Gamma (a+\kappa q^2)t}$ by the Itô isometry. A realization retains equilibrium fluctuations after its mean has decayed. Pointwise continuum white-noise quantities require a coarse-graining cutoff.

Getzler rescaling combines the change $x\mapsto \sqrt{t}x$ in geodesic normal coordinates with the Clifford algebra filtration, assigning degree one to a derivative and to a Clifford generator and degree minus one to a coordinate. The associated graded Clifford algebra is the exterior algebra. This is useful because the spinor supertrace annihilates every Clifford monomial except the top-degree one. Applied to the Lichnerowicz formula for a twisted Dirac operator, it retains a quadratic curvature model, whose heat kernel is computed by a Gaussian ansatz; all terms of smaller rescaling degree disappear from the local index limit. Thus the cancellations in a supersymmetric heat trace can be computed before solving the complete curved heat equation.

The homogeneous Dirichlet boundary condition kernel for $u_t=u_{xx}+\alpha u_x$ on $x>0$. Here $H(r,t)=e^{-r^2/(4t)}/\sqrt{4\pi t}$ is the heat kernel. The exponential reflection weight makes $K_{\alpha }(0,y,t)=0$. Gaussian decay makes its initial-data integral convergent even when the exponentially weighted initial data are not integrable separately.

For $u_t=D{u}_{xx}$ and $u(x,0)=1_{[-a,a]}(x)$, convolution with the heat kernel gives

Bounded smooth initial data alone do not guarantee a pointwise large-time limit for the heat equation on Euclidean space. In four dimensions, $g(x)=(1+|x{|}^2{)}^i$ has modulus one. At the origin its heat kernel average is asymptotic to $4^i\Gamma (2+i)t^i$, which oscillates indefinitely. Initial data vanishing at infinity do have heat averages tending to zero. This distinction matters when using Duhamel principle to construct a stationary solution for a time-independent source.

Suppose $\phi$ is continuous at $x$ and, for some $T>0$,For $0<t<T$, its convolution with the heat kernel $K_t(s)=(4\pi t{)}^{-1/2}{e}^{-s^2/(4t)}$ is smooth and solves the heat equation. On every compact set in $\mathbb{R}\times (0,T)$, spatial or time derivatives of the heat kernel are bounded by the displayed integrable weight times a constant: their polynomial factors are absorbed by the strictly narrower Gaussian. This justifies differentiation under the integral sign.

Moreover, $\int K_t=1$, soInside $|\xi -x|<\delta$, continuity makes the difference uniformly small. Outside that interval, for $t\le T/4$, the inequality ${\xi }^2\le 2(\xi -x{)}^2+2x^2$ givesThe right side tends to zero, proving that the weighted tail contribution vanishes; the constant term $\phi (x)$ is controlled by the vanishing Gaussian tail mass. Thus $(K_t\ast \phi )(x)\to \phi (x)$ as $t\downarrow 0$. In particular, continuous data of polynomial growth satisfy this hypothesis, even when they are unbounded. This extends the Gaussian approximate identity argument beyond bounded data.

For the heat equation on ${\mathbb{R}}^n$ with diffusion coefficient one, the heat kernel is $K_t(x)=(4\pi t{)}^{-n/2}{e}^{-|x{|}^2/(4t)}$. Its convolution with bounded initial data is smooth for positive time, solves the equation by differentiation under the integral, and returns continuous initial data locally uniformly through the approximate identity property. It is Gaussian filtering with per-coordinate variance $2t$.

For initial data $u_0$ on the real line, Fourier transformation or convolution with the fundamental solution gives

For a bounded function $u_0$, the convolution is well defined at every positive time because the heat kernel has integral one. The error-function diffusion front uses a Heaviside step function initial condition and gives one smoothed interface. The heat equation with interval-indicator initial data instead has two interfaces; subtracting their two error functions produces the finite-interval profile.

On a closed manifold, the positive Laplace-Beltrami operator has a smooth heat kernel and its heat semigroup is a trace-class operator for positive time. Expanding in a smooth orthonormal eigenbasis gives $K_t(x,y)={\sum }_j{e}^{-t{\lambda }_j}{\varphi }_j(x)\overline{{\varphi }_j(y)}$. Elliptic regularity bounds each derivative of ${\varphi }_j$ polynomially in ${\lambda }_j$, while the exponential dominates those powers, so the expansion converges with all derivatives. Integrating the diagonal and using the unit norm of each eigenfunction proves the formula. For a bundle operator, insert the fibre matrix trace inside the integral.

The displayed fundamental solution solves $u_t=u_{xx}-cu$ on the real line with a Dirac delta at the initial time. It is the ordinary heat kernel multiplied by $e^{-ct}$, because that factor removes the constant reaction term. Its integral is $e^{-ct}$ and it converges to the Dirac delta as $t\downarrow 0$. Convolution with initial data gives the corresponding solution. For $c>0$ the total mass decays by absorption; for $c<0$ it grows instead.

Infinite propagation speed means that a disturbance initially confined to a bounded region can affect arbitrarily distant points after any positive time. The strictly positive heat kernel gives this behavior for a nontrivial nonnegative diffusion solution. It contrasts with the finite propagation speed of the wave equation; a very small distant response is still mathematically nonzero.

A Neumann eigenfunction is an eigenfunction of a spatial differential operator whose normal derivative vanishes on the boundary. A Sturm-Liouville eigenfunction expansion in these functions gives a heat kernel for homogeneous Neumann boundary conditions.