The heat equation describes how temperature (or other diffusing quantities) spreads over time. A classical result by Li and Yau from the 1980s gives a powerful inequality controlling how fast solutions to the heat equation can vary in space and time. This inequality has been enormously useful in geometry and analysis. However, the classical result applies to settings without boundaries. When you impose a Dirichlet boundary condition, meaning the solution is forced to be zero along some boundary surface, the behavior near that boundary becomes more subtle and the classical estimate no longer captures the full picture.
This paper works in the simplest bounded-boundary setting: the half-space, which is just ordinary three-dimensional (or n-dimensional) space cut in half by a flat wall. The key tool is the explicit formula for the heat kernel in this setting, which can be written using a "reflection trick": the heat kernel with a zero-boundary condition is built by taking the usual Gaussian heat kernel and subtracting a mirror-image Gaussian reflected across the wall. The authors reinterpret logarithmic derivatives of this kernel as statistical averages (expectations) under a natural probability measure the kernel itself defines. This reframing turns a hard analytical problem into one where Jensen's inequality, a basic fact about averages of convex functions, can be applied cleanly.
The main result is a new gradient estimate showing that positive solutions satisfy a specific inequality involving two terms: the usual Li-Yau term that scales like one over time, plus a new correction term proportional to one over the square of the distance to the boundary wall. This correction blows up as you approach the wall, which makes physical sense because the solution must vanish there and therefore varies rapidly near it. The correction term arises mathematically from a hyperbolic cotangent function that appears when analyzing the reflected Gaussian. The approach is transparent about where the boundary correction comes from, and the authors suggest it could be adapted to more complicated domain shapes.