Moser's Harnack Inequality

• Moser's Harnack inequality is an interior comparison principle for nonnegative weak solutions in elliptic and parabolic settings, ensuring bounded oscillation and controlled regularity.
• The approach combines local boundedness, weak Harnack estimates, and logarithmic control via methods like John–Nirenberg and Bombieri–Giusti to convert integral bounds into pointwise estimates.
• Its extensions to fractional, quasilinear, anisotropic, and metric-measure frameworks demonstrate its broad applicability in advancing modern PDE analysis.

Moser’s Harnack inequality is the interior comparison principle asserting that a nonnegative weak solution of a uniformly elliptic or parabolic equation cannot oscillate arbitrarily on smaller domains. In the standard elliptic divergence-form setting,

$(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$

with uniformly elliptic coefficient matrix $(a_{ij})_{n\times n}$ and ellipticity constants $\lambda,\Lambda$, the classical statement is

$\sup_{Q_1} u \le C \inf_{Q_1} u,$

where $C$ depends only on $n,\lambda,\Lambda$ (Li et al., 2019). In the parabolic setting, the same principle compares the supremum on an earlier cylinder with the infimum on a later one, reflecting forward-in-time propagation (Niebel et al., 2022). Within De Giorgi–Moser theory, Harnack’s inequality sits alongside local boundedness and weak Harnack estimates as a cornerstone of regularity theory for uniformly elliptic divergence-form equations, and it is one of the central tools behind interior continuity, positivity propagation, and compactness arguments (Li et al., 2019).

1. Historical position and classical meaning

In the chronology recalled by Li and Zhang, De Giorgi proved Hölder continuity for elliptic equations in divergence form in 1957, Moser gave a new proof in 1960, and in the following year Moser obtained the Harnack inequality (Li et al., 2019). That sequence is historically important because the later note argues that the Harnack inequality was already “hidden” in De Giorgi’s 1957 work. For the parabolic theory, the 2022 trajectorial paper states that in 1971 Moser published a simplified version of his proof of the parabolic Harnack inequality, with the fundamental lemma of Bombieri and Giusti as the core new ingredient (Niebel et al., 2022).

The classical content is quantitative rather than merely qualitative. In the elliptic case, the inequality compares the maximum and minimum of a nonnegative solution on a smaller cube. In the parabolic case, it compares values on separated subcylinders, with the supremum taken earlier in time and the infimum later in time. This temporal asymmetry is not an incidental technicality: it is built into the formulation because diffusion propagates forward in time (Niebel et al., 2022).

A common oversimplification is to identify Moser’s contribution only with one final sup–inf estimate. The literature summarized here presents a more structured picture: Harnack’s inequality emerges from a chain consisting of local boundedness, weak Harnack control, and a logarithmic mechanism that links upper and lower estimates (Li et al., 2019, Niebel et al., 2022).

2. Elliptic divergence-form theorem and the “hidden” De Giorgi structure

For the standard linear uniformly elliptic equation in divergence form,

$(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$

Moser’s recalled theorem states that if $u\ge 0$ is a weak solution, then

$\sup_{Q_1} u \le C \inf_{Q_1} u,$

with $C$ depending only on $(a_{ij})_{n\times n}$0 (Li et al., 2019).

Li and Zhang identify two ingredients already present in De Giorgi’s 1957 work. The first is the local maximum estimate for nonnegative weak subsolutions,

$(a_{ij})_{n\times n}$1

and, in the more general form used in the note,

$(a_{ij})_{n\times n}$2

for any $(a_{ij})_{n\times n}$3. The second is a weak Harnack-type measure-to-point estimate for nonnegative weak supersolutions: for any $(a_{ij})_{n\times n}$4, there exists $(a_{ij})_{n\times n}$5, depending only on $(a_{ij})_{n\times n}$6, such that

$(a_{ij})_{n\times n}$7

From this, the note derives a weak Harnack inequality

$(a_{ij})_{n\times n}$8

for some $(a_{ij})_{n\times n}$9 (Li et al., 2019).

The proof architecture is short once these two estimates are available. One normalizes by assuming $\lambda,\Lambda$0, constructs a decreasing sequence $\lambda,\Lambda$1, proves the decay

$\lambda,\Lambda$2

and uses Theorem 1.3 together with a Calderón–Zygmund cube decomposition to propagate the measure decay from level $\lambda,\Lambda$3 to level $\lambda,\Lambda$4. The full Harnack inequality then follows by combining the local maximum principle with the weak Harnack estimate. The note explicitly frames this as showing that De Giorgi’s 1957 argument already contains the Harnack inequality in implicit form, and attributes that historical observation to DiBenedetto (Li et al., 2019).

3. Parabolic formulation and the Bombieri–Giusti mechanism

In Moser’s parabolic setting, the equation is

$\lambda,\Lambda$5

on a space-time cylinder $\lambda,\Lambda$6, where $\lambda,\Lambda$7 is uniformly elliptic: $\lambda,\Lambda$8 For a nonnegative weak solution in a parabolic cylinder $\lambda,\Lambda$9, the Harnack inequality compares the supremum on a backward subcylinder $\sup_{Q_1} u \le C \inf_{Q_1} u,$0 with the infimum on a later forward subcylinder $\sup_{Q_1} u \le C \inf_{Q_1} u,$1: $\sup_{Q_1} u \le C \inf_{Q_1} u,$2 with $\sup_{Q_1} u \le C \inf_{Q_1} u,$3 (Niebel et al., 2022).

The 2022 reinterpretation isolates the two quantitative inputs in Moser’s 1971 streamlined proof. The first is an $\sup_{Q_1} u \le C \inf_{Q_1} u,$4 estimate for positive solutions, obtained by iteration. The second is a weak $\sup_{Q_1} u \le C \inf_{Q_1} u,$5 estimate for $\sup_{Q_1} u \le C \inf_{Q_1} u,$6: for a positive weak supersolution $\sup_{Q_1} u \le C \inf_{Q_1} u,$7,

$\sup_{Q_1} u \le C \inf_{Q_1} u,$8

and similarly on $\sup_{Q_1} u \le C \inf_{Q_1} u,$9 for lower tails of $C$0. These two inputs are then glued by the Bombieri–Giusti lemma, which is the abstract device converting an $C$1 estimate and weak $C$2 logarithmic control into the Harnack comparison (Niebel et al., 2022).

The paper’s new contribution is geometric. It gives a trajectorial proof of the weak $C$3 estimate by following parabolic curves

$C$4

and integrating $C$5 along these paths. This replaces the classical Poincaré-based argument and interprets Moser’s logarithmic estimate as positivity propagation along parabolic trajectories. The authors emphasize that the new estimate is slightly weaker than Moser’s classical one because the cylinders must be separated by a positive time gap, but that this is still enough for Harnack since the Bombieri–Giusti mechanism only needs the estimate on well-separated sets (Niebel et al., 2022).

4. Weak Harnack, logarithmic control, and the structure of the method

Across the papers summarized here, Moser’s Harnack inequality appears not as an isolated endpoint but as the output of a recurring analytic pattern. In the elliptic note, Moser’s 1961 proof is described as an iteration scheme involving $C$6 bounds for $C$7 and $C$8, followed by the John–Nirenberg inequality to combine upper and lower estimates (Li et al., 2019). In the parabolic reinterpretation, the analogous bridge is the Bombieri–Giusti lemma rather than John–Nirenberg (Niebel et al., 2022).

This distinction clarifies the role of weak Harnack inequalities. A weak Harnack estimate is an $C$9-to-inf inequality for supersolutions; it is not yet the full sup–inf comparison. The full Harnack inequality emerges only after weak Harnack control is paired with local boundedness of subsolutions. The logarithmic estimate is the mechanism that allows these two directions to communicate.

A related misconception is to treat the logarithmic step as a technical ornament. In the papers considered here, it is indispensable. Whether formulated through weighted BMO and John–Nirenberg, through Bombieri–Giusti, or through a direct trajectorial estimate for $n,\lambda,\Lambda$0, the logarithmic control is the oscillation estimate that turns integral information into pointwise comparison (Li et al., 2019, Niebel et al., 2022).

5. Nonlocal, fractional, and metric-measure extensions

The same Harnack architecture persists in several nonclassical settings, although the statement often changes form or acquires auxiliary terms.

Setting	Representative statement	Main mechanism
Fractional Laplace with lower order terms	$n,\lambda,\Lambda$1	Caffarelli–Silvestre extension, Moser iteration, John–Nirenberg (Tan et al., 2010)
Mixed local/nonlocal parabolic equation	weak Harnack with tail term $n,\lambda,\Lambda$2	energy estimates, reverse Hölder, logarithmic estimate, Bombieri–Giusti (Garain et al., 2021)
Time-fractional diffusion	weak Harnack $n,\lambda,\Lambda$3 with optimal critical exponent	Moser iteration, Yosida approximation, Bombieri–Giusti (Zacher, 2010)
Parabolic minimizers in metric spaces	$n,\lambda,\Lambda$4	variational method, doubling, weak $n,\lambda,\Lambda$5-Poincaré, Bombieri–Giusti (Marola et al., 2013)

In the fractional Laplace case, the nonlocal equation

$n,\lambda,\Lambda$6

is converted into a degenerate local elliptic extension in one higher dimension, and the authors emphasize that no sign condition is assumed on $n,\lambda,\Lambda$7 (Tan et al., 2010). In the mixed local/nonlocal parabolic setting, the nonlocal tail is intrinsic and cannot be discarded; it records the influence of values of the solution outside the spatial cylinder (Garain et al., 2021). In the time-fractional problem, the inequality is weak rather than full, and the paper states that the full parabolic-style Harnack inequality remains open there (Zacher, 2010). In geodesic metric measure spaces, the argument becomes strictly variational and extends the sufficiency direction of the Grigor’yan–Saloff-Coste theorem from the linear heat equation to $n,\lambda,\Lambda$8 (Marola et al., 2013).

6. Quasilinear, anisotropic, subelliptic, and hypoelliptic variants

Several papers recast Moser’s Harnack philosophy in settings where the operator is quasilinear, anisotropic, degenerate, or only hypoelliptic.

For the quasilinear equation

$n,\lambda,\Lambda$9

the 2016 paper proves a Harnack comparison inequality for positive $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$0 weak solutions by adapting Moser iteration in the style of Trudinger; the conclusion is a comparison estimate for $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$1, not the classical single-function theorem (Merchán et al., 2016). The 2025 anisotropic extension replaces $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$2 by a Finsler norm $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$3 and studies

$(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$4

again via Moser iteration, obtaining a Harnack type comparison inequality and a Harnack type inequality for the linearized operator (Vuono, 11 Jan 2025).

In the Grushin setting, the subelliptic nondivergence-form equation

$(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$5

with $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$6 and $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$7 is treated by a completely different route: weighted ABP, critical density, double ball property, power decay property, and then the axiomatic Harnack theory of Di Fazio–Gutiérrez–Lanconelli (Montanari, 2014). For the hypoelliptic operator

$(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$8

the 2011 paper proves a Moser-style Harnack inequality for bounded nonnegative distributional solutions, but by a probabilistic proof based on the Feynman–Kac formula and hypoelliptic density estimates rather than Moser iteration (Hamel et al., 2011). These examples show that “Moser-style” may refer to the form of the conclusion and the regularity principle, even when the proof mechanism is no longer iterative.

A further variation appears in the one-dimensional singular parabolic $(a_{ij}u_i)_j=0 \quad \text{in } Q_6,$9-Laplacian, $u\ge 0$0, where the paper proves a Harnack-type “sidewise spreading of positivity” estimate. It explicitly states that this is not the classical pointwise sup/inf relation of Moser, but a propagation-of-positivity estimate with spatial decay $u\ge 0$1 (Düzgün et al., 2015).

7. Refined analogues, discrete versions, and the scope of the name

The name “Harnack inequality” also appears in settings that are adjacent to, but not identical with, Moser’s theorem. For positive harmonic functions on the unit disc, the 2025 paper proves a strengthened inequality depending on both $u\ge 0$2 and $u\ge 0$3: $u\ge 0$4 The paper explicitly states that this is not Moser’s theorem in full generality, but a harmonic, unit-disc, gradient-refined Harnack inequality (Svetlik, 16 Jan 2025).

A different analogue appears in the Jordan–Kinderlehrer–Otto scheme for the heat equation on the flat torus. There the paper proves a discrete-time Harnack inequality for the JKO approximation, derived from a matrix lower bound on $u\ge 0$5, and notes that as $u\ge 0$6 one recovers Hamilton’s matrix differential Harnack inequality and then the heat-equation Harnack inequality itself (Lee, 2015). This situates the classical Moser theory alongside Hamilton/Li–Yau differential inequalities and optimal-transport discretizations.

By contrast, the determinantal inequalities for contractive matrices studied in matrix analysis are explicitly described as “Harnack-type determinantal inequalities,” not Moser’s Harnack inequality from PDE and geometry (Lin et al., 2017). That distinction matters terminologically. “Harnack-type” may indicate formal analogy, whereas “Moser’s Harnack inequality” conventionally refers to the sup–inf estimates for positive solutions of elliptic or parabolic equations and to the analytic machinery—iteration, logarithmic control, and positivity propagation—that produces them.