---
title: Li–Yau Inequality Overview
url: https://www.emergentmind.com/topics/li-yau-inequality
type: topic
---

# Li–Yau Inequality Overview

The Li–Yau inequality is a differential Harnack inequality for positive solutions of heat-type equations. In its classical form, if \(u>0\) solves \(\partial_t u=\Delta u\) on a complete \(n\)-dimensional Riemannian manifold with nonnegative Ricci curvature, then
\[
|\nabla \log u|^2-\partial_t \log u \le \frac{n}{2t},
\]
equivalently,
\[
\frac{|\nabla u|^2}{u^2}-\frac{\partial_t u}{u}\le \frac{n}{2t}.
\]
The estimate is sharp in Euclidean space, where equality is achieved by the heat kernel, and it is one of the basic mechanisms behind parabolic Harnack inequalities, heat kernel bounds, and Liouville theorems [1405.0684][2501.12685].

## 1. Classical heat-equation formulation

Li and Yau introduced the inequality in the setting of positive heat solutions on complete manifolds under lower Ricci curvature bounds. In the nonnegative-curvature case, the estimate controls the logarithmic gradient by the parabolic scale \(t^{-1}\), and the right-hand side \(\frac{n}{2t}\) is the canonical dimensional constant [1405.0684]. A standard curved version states that if \(\mathrm{Ric}\ge -K\) with \(K\ge 0\), then for every \(\alpha>1\),
\[
\frac{|\nabla u|^2}{u^2}-\alpha \frac{\partial_t u}{u}
\le
\frac{n\alpha^2}{2t}+\frac{n\alpha^2K}{2(\alpha-1)},
\]
which reduces to the sharp Euclidean formula when \(K=0\) and \(\alpha=1\) [2105.03609].

The estimate is often written in Laplacian form. Since
\[
\Delta \log u=\frac{\Delta u}{u}-\frac{|\nabla u|^2}{u^2},
\]
the classical inequality can also be stated as
\[
-\Delta \log u \le \frac{n}{2t},
\]
or, after substituting \(\partial_t u=\Delta u\), as a control of \(\partial_t \log u\) by \( |\nabla \log u|^2\) and \(t^{-1}\) [2312.16181]. This multiplicity of equivalent formulations is one reason the Li–Yau inequality appears in several distinct analytic guises.

## 2. Analytical structure and proof mechanisms

The original proof is based on differentiating the heat equation, applying Bochner’s identity, and using the parabolic maximum principle on a carefully chosen quantity such as \(F=|\nabla f|^2-\alpha f_t\) with \(f=\log u\) [2105.03609]. In smooth settings this route is naturally encoded by the Bakry–Émery curvature-dimension formalism, where \(CD(K,N)\) or related conditions provide the synthetic replacement for lower Ricci curvature bounds [1412.5165].

A major later development is that Li–Yau theory can be reformulated at semigroup level. Under a Markov diffusion triple satisfying \(CD(\rho,n)\), one obtains a global Li–Yau inequality expressed through an explicit concave function \(\Phi_t\) of \(\frac{LP_t f}{P_t f}\); the authors state that this inequality is stronger than all classical Li–Yau type inequalities known to them, and on Riemannian manifolds it is equivalent to a new parabolic Harnack inequality in both negative and positive curvature [1412.5165]. This replaces the classical pointwise maximum-principle derivation by an entropy and Sturm-comparison argument at the semigroup level.

A second major proof paradigm uses explicit heat-kernel representations rather than curvature-dimension inequalities. For very general non-local diffusion equations, a reduction principle shows that Li–Yau bounds for the heat kernel imply corresponding bounds for all positive solutions represented by that kernel [2012.12974]. The same representation-based viewpoint also yields generalized second- and fourth-derivative Li–Yau type inequalities on \(\mathbb{R}^n\) [2312.16181]. This suggests that Li–Yau theory is not tied to a single differential identity; it can also emerge from precise kernel calculus.

## 3. Consequences for Harnack theory, heat kernels, and rigidity

The canonical consequence of the Li–Yau inequality is the parabolic Harnack inequality. Integrating the differential bound along space-time curves gives estimates relating \(u(x,s)\) and \(u(y,t)\) for \(0<s<t\), typically with an exponential factor involving \(\frac{d(x,y)^2}{t-s}\) and a time-ratio term such as \((t/s)^{N/2}\) [1405.0684]. This is the standard Li–Yau mechanism in both smooth and synthetic settings.

Heat kernel bounds follow from the same scheme. In smooth and \(RCD^*(K,N)\) settings, Li–Yau and Harnack inequalities are consistent with Gaussian two-sided estimates; in the \(RCD^*(0,N)\) case the heat kernel itself satisfies the Li–Yau gradient estimate
\[
|\nabla_x \log p_t(x,y)|^2-\partial_t \log p_t(x,y)\le \frac{N}{2t}
\]
for \(\mu\)-a.e. \(x,y\) [1405.0684]. In the abstract Markov-semigroup framework, the global Li–Yau inequality also yields new heat-kernel Harnack bounds and, under positive curvature, ultracontractive estimates by a direct method [1412.5165].

Rigidity phenomena are another standard consequence. In classical Li–Yau theory, bounded harmonic functions on complete manifolds with nonnegative Ricci curvature are constant; analogous Liouville properties persist in several extensions, including discrete graphs and metric star graphs [2501.12685][1311.3367]. Spectral applications are also typical: on graphs, Li–Yau inequalities lead to heat-kernel bounds, Buser-type inequalities, and Cheng-type eigenvalue estimates [1306.2561][1801.06021].

## 4. Synthetic, discrete, and non-local generalizations

One major branch of the subject extends Li–Yau from smooth manifolds to synthetic metric-measure geometry. On \(RCD^*(0,N)\) spaces, the heat flow satisfies exactly the classical-form estimate
\[
|\nabla \log u(x,T)|^2-\partial_T\log u(x,T)\le \frac{N}{2T},
\]
for nonnegative \(L^1\cap L^\infty\) initial data, and the corresponding Harnack and heat-kernel inequalities persist for general \(K\in\mathbb{R}\) through a Baudoin–Garofalo-type refinement [1405.0684].

A second branch concerns graphs. Several discrete curvature frameworks support Li–Yau theory: the exponential curvature-dimension condition \(CDE(n,K)\) yields gradient estimates, Harnack inequalities, heat kernel bounds, polynomial volume growth under nonnegative curvature, and Buser-type inequalities [1306.2561]; \(CDE(n,-K)\) supports time-dependent Li–Yau families and Hamilton-type estimates [1311.3367]; \(CDE'(n,K)\) allows analogous results for unbounded graph Laplacians, including Cheng’s eigenvalue estimate [1801.06021]. A distinct nonlinear curvature-dimension condition \(CD_\psi(d,0)\) recovers an exact logarithmic Li–Yau inequality on finite graphs, and Ricci-flat graphs have nonnegative curvature in this sense [1412.3340]. Another graph-theoretic approach replaces \(\log P_t f\) by a modified nonlinear flow \(\partial_t u=\Delta u+\Gamma u\); under \(CD(0,n)\) this yields
\[
-\Delta u_t\le \frac{n}{2t},
\]
from which volume doubling is deduced, solving a major open problem in discrete Ricci curvature, and implying that there exist no expander graphs satisfying \(CD(0,n)\) [1909.10242].

Discrete curvature is not the only route. Cayley graphs of virtually Abelian groups satisfy Li–Yau type gradient estimates even though they do not satisfy any known variant of the curvature-dimension inequality with non-negative curvature [1610.05227]. For non-local diffusions, a heat-kernel reduction principle yields Li–Yau inequalities for the fractional heat equation,
\[
(-\Delta)^{\beta/2}(\log u)\le \frac{C_{\mathrm{LY}(\beta,d)}}{t},
\]
thereby solving the problem of obtaining a \(t^{-1}\)-scale Li–Yau inequality for positive solutions of the fractional heat equation [2012.12974].

## 5. Metric graphs and singular one-dimensional geometry

Metric graphs interpolate between manifolds and discrete graphs: edges carry continuous one-dimensional geometry, while vertices introduce singularities. For a metric star graph \(\mathcal G\) consisting of \(N\) half-lines joined at one central vertex and equipped with weighted Kirchhoff matching
\[
\sum_{j=1}^N a_j \partial_x u_j(0,t)=0,\qquad a_j>0,\quad \sum_{j=1}^N a_j=1,
\]
positive heat solutions admit an explicit kernel representation [2501.12685].

In that setting, a Li–Yau type estimate holds on each edge:
\[
\partial_t \ln u(x,t)-|\partial_x\ln u(x,t)|^2
\ge
-\frac{1}{2t}-(1-2a_i)L_i(x,t),
\]
where \(L_i(x,t)\ge 0\) is an explicit integral term depending on the initial data, the edge, the point, the time, and the coefficient \(a_i\) [2501.12685]. The correction term is the analytic signature of the branching singularity and the weighted transmission law at the vertex.

This formula reduces sharply to the classical one-dimensional Li–Yau inequality on \(\mathbb R\) when the star graph degenerates to two symmetric half-lines with \(a_1=a_2=\tfrac12\); in that case the extra term vanishes [2501.12685]. The same estimate implies a parabolic Harnack inequality on the graph and a Liouville property: every bounded harmonic function satisfying the vertex conditions is constant [2501.12685]. A plausible implication is that, on singular one-dimensional spaces, Li–Yau theory remains valid but generally acquires an explicit defect term encoding vertex scattering.

## 6. Higher-order, semilinear, and geometric variants

The Li–Yau paradigm has also been extended beyond the linear heat equation. On \(\mathbb R^n\), a representation-formula method yields a generalized second-derivative family,
\[
\frac{\Delta u}{u}
-\alpha \sum_{i\ne j}\frac{u_{x_i x_j}}{u}
-\beta \sum_{i\ne j}\frac{u_{x_i}u_{x_j}}{u^2}
-\gamma \frac{|\nabla u|^2}{u^2}
\ge -\frac{n}{2t},
\]
under the parameter constraint \((n-1)(\alpha+\beta)+\gamma\le 1\), and also a fourth-order Li–Yau type inequality with \(t^{-2}\) scaling [2312.16181]. These results show that higher-order Li–Yau structures can be derived directly from Gaussian kernel calculus.

For the semilinear heat equation
\[
u_t=\Delta u+u^p
\]
on complete manifolds with \(\mathrm{Ric}\ge 0\), a semilinear Li–Yau inequality holds in the exponent range \(p\in(0,\overline p_n)\). In particular, for admissible \((\alpha,\beta)\),
\[
u_t-\alpha\frac{|\nabla u|^2}{u}-\beta u^p
\ge -\frac{u}{\varepsilon t},
\]
which leads to a Harnack inequality and to monotonicity, convexity, decay estimates, and triviality statements for ancient and eternal solutions [2201.02530].

The name “Li–Yau inequality” also appears in geometric analysis in a different but related sense. For closed immersed surfaces in \(\mathbb R^3\), the classical Willmore-theoretic Li–Yau inequality states that \(W(f)<8\pi\) implies embeddedness; Simon’s monotonicity formula extends this principle to surfaces with boundary, producing boundary-dependent thresholds \(C_{\mathrm{LY}}\) and \(C_{\mathrm{LY}}^{\mathrm{rot}}\) for embeddedness or exclusion of axis self-intersections [2402.07755]. In dimension one, an analogous embeddedness threshold for closed planar curves is expressed using the scale-invariant elastic energy \(\mathcal E(\gamma)\mathcal L(\gamma)\), with the figure-eight elastica furnishing the sharp non-embedded threshold [2101.08509]. This suggests that “Li–Yau inequality” now denotes a broader multiplicity-control paradigm: in heat-flow settings it controls oscillation through differential Harnack estimates, while in Willmore-type settings it controls self-intersection through curvature-energy thresholds.

Source: https://www.emergentmind.com/topics/li-yau-inequality