The Li-Yau inequality is a classical and powerful result in geometric analysis that controls how fast solutions to heat-type equations can change in space and time. It has deep consequences, including Harnack inequalities that compare the values of a solution at different points and times. The original result applies to the standard heat equation on curved spaces, but this paper extends it to a much more general family of equations called the normalized parabolic p-Laplacian, which interpolates between the heat equation and equations more sensitive to the gradient of the solution. These equations arise naturally in image processing, optimal transport, and the study of geometric flows.
The main achievement is proving a sharp version of the Li-Yau inequality for this nonlinear family, meaning the constant appearing in the inequality cannot be improved. The authors identify an explicit self-similar solution on flat Euclidean space that achieves equality in the bound, and they show this sharpness carries over to compact curved spaces by taking larger and larger flat tori that approximate Euclidean space. A key technical challenge is that the p-Laplacian equation behaves badly at points where the gradient of the solution vanishes, so many standard tools break down there. The authors sidestep this by working with a smoothed approximating family of equations where this problem does not arise, proving the inequality uniformly in the approximation parameter, and then passing to a limit.
The proof introduces a precise algebraic identity of Bochner type, a technique that encodes geometric and analytic information about how curvature interacts with the equation. This identity turns out to be stable under a classical deformation of the Li-Yau functional, which allows the authors to handle manifolds with negative as well as nonnegative curvature bounds. As a consequence, they derive a sharp global Harnack inequality, which concretely means one can compare the solution's values at widely separated points in a quantitatively optimal way. The results unify and sharpen several earlier partial results and work without any extra assumptions on where the solution's gradient vanishes.