Gaussian heat kernel solves the 1D heat equation

← All problems

heat_kernel_solves_heat_equation

Submitter: Kim Morrison.

Notes: The Gaussian convolution u(t,x) = (4 pi t)^{-1/2} integral exp(-(x-y)^2/(4t)) f(y) dy satisfies the heat equation on (0, infty) x R, with initial datum f recovered as t -> 0+.

Source: Classical heat-equation theory.

Informal solution: Differentiate under the integral sign in t and twice in x; the kernel itself satisfies the heat equation, so does its convolution with f. As t -> 0+, the Gaussian is an approximate identity with total mass 1. After the change of variables y = x + sqrt(t) z, the integral becomes an average of f(x + sqrt(t) z) against a fixed integrable Gaussian weight; continuity of f at x and boundedness of f then give pointwise convergence to f(x) by dominated convergence in z.

theorem declaration uses `sorry`heat_kernel_solves_heat_equation (f : ) (hf_cont : Continuous f) (hf_bdd : M : , x, |f x| M) : -- The PDE on (0, ∞) × ℝ. ( t : , 0 < t x : , ux : , uxx : , ( y : , HasDerivAt (fun z => heatSolution f t z) (ux y) y) HasDerivAt ux uxx x HasDerivAt (fun s => heatSolution f s x) uxx t) -- Initial condition recovered as a one-sided limit at t = 0. ( x : , Filter.Tendsto (fun t : => heatSolution f t x) (nhdsWithin (0 : ) (Set.Ioi 0)) (nhds (f x))) := f: hf_cont:Continuous fhf_bdd: M, (x : ), |f x| M(∀ (t : ), 0 < t (x : ), ux uxx, (∀ (y : ), HasDerivAt (fun z => heatSolution f t z) (ux y) y) HasDerivAt ux uxx x HasDerivAt (fun s => heatSolution f s x) uxx t) (x : ), Filter.Tendsto (fun t => heatSolution f t x) (nhdsWithin 0 (Set.Ioi 0)) (nhds (f x)) All goals completed! 🐙

Solved by

@mayorov-m-a with Aleph Prover(logicalintelligence.com) on May 8, 2026

@LorenzoLuccioli with Aristotle (Harmonic) on May 8, 2026 (proof)

@daouid with Antigravity (Multi-Model Ensemble: Gemini 3.1 Pro, Gemini 3 Flash, Claude 4.6 Sonnet/Opus) on May 13, 2026 (proof)

@GanjinZero with Seed Prover (ByteDance) on May 20, 2026