Book 12A

Β© 2026 NeckPinch Β· www.neckpinch.com Β· All rights reserved.

Course 12Book 12A: Entropy and NoncollapsingChapter 6

𝓛-Geodesics and the Variational Formulas

Perelman's Β§7.1 derived line by line: the first and second variation of 𝓛-length, 𝓛-Jacobi fields, adapted fields, Hamilton's Harnack expressions, and the differential inequalities for the reduced distance.

51 min read Β· Updated Oct 23, 2026

Read Perelman I, Β§7.1, equations (7.1)–(7.15), with Morgan and Tian, Ricci Flow and the PoincarΓ© Conjecture, chapter 6, Β§Β§1–5 (Lemma 6.4 to Corollary 6.51), which carry out every computation Perelman compresses. Kleiner and Lott's notes, the sections on L\mathcal L-geodesics and reduced length, are a second account. For the Riemannian model, read the second variation of length and the Hessian comparison for the distance function in 9A.7 Jacobi Fields and Curvature versus Topology and 9B.1 Laplacian Comparison first; this chapter follows them step by step.

In this chapter Β· 11 sections
  1. 6.1The setting and the variable s=Ο„s = \sqrt\taus=τ​
  2. 6.2The first variation and the L\mathcal LL-geodesic equation
  3. 6.3The gradient and the time derivative of LLL
  4. 6.4Hamilton's trace Harnack expression appears
  5. 6.5The second variation and L\mathcal LL-Jacobi fields
  6. 6.6The Hessian comparison with adapted fields
  7. 6.7The trace: the Laplacian comparison
  8. 6.8Combining: the inequalities (7.13)–(7.15)
  9. 6.9When the curvature operator is nonnegative
  10. 6.10History
  11. 6.11Exercises
  12. β†—Free resources

12A.5 Reduced Distance and Reduced Volume defined the L\mathcal L-length, the reduced distance β„“\ell and the reduced volume V~\tilde V, stated the differential inequalities that β„“\ell satisfies, and showed how they make V~\tilde V monotone. It took those inequalities on trust. This chapter proves them. The proofs are the content of a single page of Perelman's first preprint, Β§7.1, which he describes as "natural modifications of those in the classical variational theory of geodesics that can be found in any textbook on Riemannian geometry". That description is accurate and also misleading: every step is classical, but there are many steps, the signs are treacherous, and the time dependence of the metric enters every one of them, producing Ricci terms that must be tracked and, at the end, recognised as Hamilton's Harnack expressions. Morgan and Tian spend twenty pages on the same material. A reader who has done the second variation of arc length once (9A.7 Jacobi Fields and Curvature versus Topology) and the Laplacian comparison for the distance function once (9B.1 Laplacian Comparison) has all the tools; what remains is to see how they adapt.

The chapter follows Perelman's equations in order, (7.1) to (7.15), and gives each its derivation. The one new idea, compared with Riemannian geometry, is the choice of test fields in the second variation: Perelman's fields YY solving βˆ‡XY=βˆ’Ric⁑(Y,β‹…)+12Ο„Y\nabla_XY = -\operatorname{Ric}(Y, \cdot) + \frac{1}{2\tau}Y, which play the role that the linear Jacobi fields srE\frac srE play in Euclidean space. With them the second variation of L\mathcal L becomes the integral of Hamilton's matrix Harnack expression, and its trace becomes Hamilton's trace Harnack expression H(X)H(X). Perelman's quantity K=βˆ«Ο„3/2H(X) dΟ„K = \int\tau^{3/2}H(X)\,d\tau is then the only unknown, and it cancels from the combinations that matter. 12A.7 Monotonicity of the Reduced Volume uses the result to prove the monotonicity of the reduced volume and the bound min⁑ℓ≀n2\min\ell \leq \frac n2 rigorously.

By the end of this chapter you will be able to:

  • write the L\mathcal L-length in the variable s=Ο„s = \sqrt\tau, and explain why the L\mathcal L-geodesic equation is regular at Ο„=0\tau = 0;
  • derive the first variation formula (7.1) and the L\mathcal L-geodesic equation (7.2), including the Ricci term that comes from differentiating the metric;
  • prove that βˆ‡L=2Ο„X\nabla L = 2\sqrt\tau X and βˆ‚Ο„L=2Ο„Rβˆ’Ο„(R+∣X∣2)\partial_\tau L = 2\sqrt\tau R - \sqrt\tau(R + |X|^2);
  • derive (7.3)–(7.6), identify H(X)H(X) with Hamilton's trace Harnack expression at t=βˆ’Ο„t = -\tau, and express βˆ‚Ο„β„“\partial_\tau\ell and βˆ£βˆ‡β„“βˆ£2|\nabla\ell|^2 through KK;
  • derive the second variation formula (7.7) and the L\mathcal L-Jacobi equation;
  • prove the Hessian comparison (7.9) with adapted fields, take its trace to get the Laplacian comparison (7.10), and verify the trace identity that turns the matrix Harnack expression into the trace one;
  • combine everything into (7.13)–(7.15), and identify the equality case as a gradient shrinking soliton;
  • state Perelman's Β§7.2 estimates under nonnegative curvature operator, which 12B.3 The Structure of ΞΊ-Solutions uses for ΞΊ\kappa-solutions.

The setting and the variable s=Ο„s = \sqrt\tau

Let g(t)g(t) be a Ricci flow on MM, either with MM closed or with each g(t)g(t) complete of uniformly bounded curvature, and fix a time TT. Write Ο„=Tβˆ’t\tau = T - t for backward time, so that

βˆ‚gβˆ‚Ο„=2Ric⁑.\frac{\partial g}{\partial\tau} = 2\operatorname{Ric}.

A path parameterised by backward time is a curve Ξ³:[Ο„1,Ο„2]β†’M\gamma : [\tau_1, \tau_2] \to M, thought of as the space-time path τ↦(Ξ³(Ο„),Tβˆ’Ο„)\tau \mapsto (\gamma(\tau), T - \tau), which moves backward in tt as Ο„\tau increases. Write X=Ξ³Λ™=dΞ³dΟ„X = \dot\gamma = \frac{d\gamma}{d\tau}, and measure it with g(Ο„)g(\tau). The L\mathcal L-length is

L(Ξ³)=βˆ«Ο„1Ο„2Ο„(R(Ξ³(Ο„),Ο„)+∣X(Ο„)∣g(Ο„)2) dΟ„.\mathcal L(\gamma) = \int_{\tau_1}^{\tau_2}\sqrt\tau\big(R(\gamma(\tau), \tau) + |X(\tau)|^2_{g(\tau)}\big)\,d\tau.

Morgan and Tian work more generally with generalised Ricci flows, space-times that are locally isomorphic to Ricci flows but may change topology, which is what a flow with surgery is (12B.9 Surgery). Every computation in this chapter is local along a path, so it applies verbatim there; that is why they set it up in that generality.

The variable ss. Set s=Ο„s = \sqrt\tau and A(s)=dΞ³dsA(s) = \frac{d\gamma}{ds}. Then dΟ„=2s dsd\tau = 2s\,ds and X=A2sX = \frac{A}{2s}, that is, A=2Ο„XA = 2\sqrt\tau X. Substituting,

Ο„R dΟ„=2s2R ds,Ο„βˆ£X∣2dΟ„=s∣A∣24s22s ds=12∣A∣2 ds,\sqrt\tau R\,d\tau = 2s^2R\,ds, \qquad \sqrt\tau|X|^2d\tau = s\frac{|A|^2}{4s^2}2s\,ds = \frac12|A|^2\,ds,

so

L(Ξ³)=βˆ«Ο„1Ο„2(12∣A(s)∣2+2s2R(Ξ³(s),s2)) ds.\mathcal L(\gamma) = \int_{\sqrt{\tau_1}}^{\sqrt{\tau_2}}\Big(\frac12|A(s)|^2 + 2s^2R(\gamma(s), s^2)\Big)\,ds.

This is an energy plus a potential, and in this form the theory looks like the classical calculus of variations (6A.11 Calculus of Variations and Gradient Flows): the weight Ο„\sqrt\tau, which seemed strange, is exactly what makes L\mathcal L a standard energy in ss. On static flat space, where R=0R = 0, minimisers are straight lines at constant ss-speed, Ξ³(s)=p+2sv\gamma(s) = p + 2sv, that is, Ξ³(Ο„)=p+2τ v\gamma(\tau) = p + 2\sqrt\tau\,v: in Ο„\tau they leave pp with infinite speed and slow down, like the fan of curves in 12A.5 Reduced Distance and Reduced Volume.

The first variation and the L\mathcal L-geodesic equation

Let Ξ³u(Ο„)=Ξ³(Ο„,u)\gamma_u(\tau) = \gamma(\tau, u) be a smooth family of paths parameterised by backward time, with Ξ³0=Ξ³\gamma_0 = \gamma, and let Y=βˆ‚uΞ³Y = \partial_u\gamma and X=βˆ‚Ο„Ξ³X = \partial_\tau\gamma on the surface swept out. Because these are coordinate vector fields on that surface, [X,Y]=0[X, Y] = 0, so βˆ‡YX=βˆ‡XY\nabla_YX = \nabla_XY (the connection being torsion-free; the connection of g(Ο„)g(\tau) is used at each Ο„\tau, and its time dependence does not affect this identity, which holds at each fixed Ο„\tau).

One new identity. For vector fields VV, WW along Ξ³\gamma, the derivative of their inner product along the path has an extra term, because the metric changes:

ddΟ„βŸ¨V,W⟩g(Ο„)=βŸ¨βˆ‡XV,W⟩+⟨V,βˆ‡XW⟩+2Ric⁑(V,W).\frac{d}{d\tau}\langle V, W\rangle_{g(\tau)} = \langle\nabla_XV, W\rangle + \langle V, \nabla_XW\rangle + 2\operatorname{Ric}(V, W).

The first two terms are the usual Leibniz rule for the Levi-Civita connection at fixed Ο„\tau; the last is βˆ‚gβˆ‚Ο„(V,W)=2Ric⁑(V,W)\frac{\partial g}{\partial\tau}(V, W) = 2\operatorname{Ric}(V, W). Every Ricci term in this chapter comes from this line.

Now compute, following Perelman's (7.1). Differentiating under the integral,

Ξ΄YL=βˆ«Ο„1Ο„2Ο„(⟨Y,βˆ‡R⟩+2βŸ¨βˆ‡YX,X⟩)dΟ„=βˆ«Ο„1Ο„2Ο„(⟨Y,βˆ‡R⟩+2βŸ¨βˆ‡XY,X⟩)dΟ„.\delta_Y\mathcal L = \int_{\tau_1}^{\tau_2}\sqrt\tau\big(\langle Y, \nabla R\rangle + 2\langle\nabla_YX, X\rangle\big)d\tau = \int_{\tau_1}^{\tau_2}\sqrt\tau\big(\langle Y, \nabla R\rangle + 2\langle\nabla_XY, X\rangle\big)d\tau.

By the identity, βŸ¨βˆ‡XY,X⟩=ddΟ„βŸ¨Y,XβŸ©βˆ’βŸ¨Y,βˆ‡XXβŸ©βˆ’2Ric⁑(Y,X)\langle\nabla_XY, X\rangle = \frac{d}{d\tau}\langle Y, X\rangle - \langle Y, \nabla_XX\rangle - 2\operatorname{Ric}(Y, X). Substituting, and integrating 2Ο„ddΟ„βŸ¨Y,X⟩2\sqrt\tau\frac{d}{d\tau}\langle Y, X\rangle by parts (with ddΟ„2Ο„=1Ο„\frac{d}{d\tau}2\sqrt\tau = \frac{1}{\sqrt\tau}):

Ξ΄YL=2Ο„βŸ¨X,YβŸ©βˆ£Ο„1Ο„2+βˆ«Ο„1Ο„2Ο„βŸ¨Y,Β βˆ‡Rβˆ’2βˆ‡XXβˆ’4Ric⁑(X,β‹…)βˆ’1Ο„X⟩dΟ„.(7.1)\delta_Y\mathcal L = 2\sqrt\tau\langle X, Y\rangle\Big|_{\tau_1}^{\tau_2} + \int_{\tau_1}^{\tau_2}\sqrt\tau\Big\langle Y,\ \nabla R - 2\nabla_XX - 4\operatorname{Ric}(X, \cdot) - \frac1\tau X\Big\rangle d\tau. \tag{7.1}

Here Ric⁑(X,β‹…)\operatorname{Ric}(X, \cdot) denotes the vector dual to the one-form Ric⁑(X,β‹…)\operatorname{Ric}(X, \cdot). For variations fixing the endpoints the boundary term vanishes, and Ξ³\gamma is critical for all of them exactly when the bracket vanishes:

βˆ‡XXβˆ’12βˆ‡R+12Ο„X+2Ric⁑(X,β‹…)=0.(7.2)\nabla_XX - \frac12\nabla R + \frac{1}{2\tau}X + 2\operatorname{Ric}(X, \cdot) = 0. \tag{7.2}

This is the L\mathcal L-geodesic equation. Compare the Riemannian geodesic equation βˆ‡XX=0\nabla_XX = 0: the term 12βˆ‡R\frac12\nabla R is a force pushing the path toward larger scalar curvature (the potential 2s2R2s^2R in the ss-form, with the sign of a Lagrangian), the term 12Ο„X\frac{1}{2\tau}X is the "slowing down" that makes Ο„X\sqrt\tau X rather than XX the natural velocity, and 2Ric⁑(X,β‹…)2\operatorname{Ric}(X, \cdot) is the correction from the moving metric.

Regularity at Ο„=0\tau = 0. In the variable ss, with X=A2sX = \frac{A}{2s} and ddΟ„=12sdds\frac{d}{d\tau} = \frac{1}{2s}\frac{d}{ds}, one finds βˆ‡XX=14s2βˆ‡AAβˆ’14s3A\nabla_XX = \frac{1}{4s^2}\nabla_AA - \frac{1}{4s^3}A, and (7.2) becomes, after multiplying by 4s24s^2 (Exercise 6.1),

βˆ‡AAβˆ’2s2βˆ‡R+4sRic⁑(A,β‹…)=0.\nabla_AA - 2s^2\nabla R + 4s\operatorname{Ric}(A, \cdot) = 0.

This ODE has smooth coefficients at s=0s = 0. So for every v∈TpMv \in T_pM there is a unique L\mathcal L-geodesic with Ξ³(0)=p\gamma(0) = p and A(0)=2vA(0) = 2v, that is, lim⁑τ→0Ο„X(Ο„)=v\lim_{\tau\to0}\sqrt\tau X(\tau) = v (Morgan–Tian, Lemma 6.8). The L\mathcal L-exponential map is Lexp⁑τ(v)=Ξ³v(Ο„)\mathcal L\exp_\tau(v) = \gamma_v(\tau). On a complete flow of bounded curvature, L\mathcal L-geodesics extend to the whole time interval (their Lemma 7.3), and minimising curves exist between any two space-time points with Ο„1<Ο„2\tau_1 < \tau_2 (their Lemma 7.2): a minimising sequence has bounded energy ∫∣A∣2ds\int|A|^2ds in the ss-variable, hence is uniformly HΓΆlder continuous, a subsequence converges weakly in H1H^1, and L\mathcal L is lower semicontinuous, exactly the direct method of 6A.11 Calculus of Variations and Gradient Flows and 4A.7 Weak Convergence and the Direct Method. The minimiser satisfies (7.2) and is smooth.

The gradient and the time derivative of LL

Fix a base point pp at Ο„=0\tau = 0 and write L(q,Ο„Λ‰)L(q, \bar\tau) for the infimum of L\mathcal L over paths from (p,0)(p, 0) to (q,Ο„Λ‰)(q, \bar\tau). Pretend, as Perelman does, that the minimiser Ξ³\gamma is unique and depends smoothly on (q,Ο„Λ‰)(q, \bar\tau); 12A.7 Monotonicity of the Reduced Volume shows that this is true on an open set of full measure (the analogue of the complement of the cut locus), and that elsewhere the inequalities hold in the barrier sense.

The gradient. Vary the endpoint qq with velocity Y(Ο„Λ‰)Y(\bar\tau) through minimisers. The integral term of (7.1) vanishes (Ξ³\gamma is an L\mathcal L-geodesic) and the boundary term at Ο„=0\tau = 0 vanishes (Y(0)=0Y(0) = 0, and Ο„X\sqrt\tau X stays bounded), so Ξ΄YL=2Ο„Λ‰βŸ¨X(Ο„Λ‰),Y(Ο„Λ‰)⟩\delta_YL = 2\sqrt{\bar\tau}\langle X(\bar\tau), Y(\bar\tau)\rangle, that is,

βˆ‡L(q,Ο„Λ‰)=2Ο„Λ‰X(Ο„Λ‰),βˆ£βˆ‡L∣2=4Ο„Λ‰βˆ£X∣2=βˆ’4Ο„Λ‰R+4Ο„Λ‰(R+∣X∣2).\nabla L(q, \bar\tau) = 2\sqrt{\bar\tau}X(\bar\tau), \qquad |\nabla L|^2 = 4\bar\tau|X|^2 = -4\bar\tau R + 4\bar\tau(R + |X|^2).

This is the analogue of βˆ‡d(p,β‹…)=Ξ³Λ™\nabla d(p, \cdot) = \dot\gamma for the Riemannian distance: the gradient of the length function at the endpoint is the velocity of the minimiser, here rescaled.

The time derivative. Along the minimiser, ddΟ„L(Ξ³(Ο„),Ο„)=Ο„(R+∣X∣2)\frac{d}{d\tau}L(\gamma(\tau), \tau) = \sqrt\tau(R + |X|^2) by the fundamental theorem of calculus (an initial segment of a minimiser is a minimiser, as for geodesics). On the other hand, by the chain rule, ddΟ„L(Ξ³(Ο„),Ο„)=βˆ‚Ο„L+βŸ¨βˆ‡L,X⟩\frac{d}{d\tau}L(\gamma(\tau), \tau) = \partial_\tau L + \langle\nabla L, X\rangle. Hence

βˆ‚Ο„Λ‰L=Ο„Λ‰(R+∣X∣2)βˆ’βŸ¨X,βˆ‡L⟩=Ο„Λ‰(R+∣X∣2)βˆ’2Ο„Λ‰βˆ£X∣2=2Ο„Λ‰Rβˆ’Ο„Λ‰(R+∣X∣2).\partial_{\bar\tau}L = \sqrt{\bar\tau}(R + |X|^2) - \langle X, \nabla L\rangle = \sqrt{\bar\tau}(R + |X|^2) - 2\sqrt{\bar\tau}|X|^2 = 2\sqrt{\bar\tau}R - \sqrt{\bar\tau}(R + |X|^2).

(Morgan–Tian, Lemma 6.32.) Both formulas involve the unknown quantity R+∣X∣2R + |X|^2 at the endpoint, and the next step computes it.

Hamilton's trace Harnack expression appears

Differentiate R+∣X∣2R + |X|^2 along an L\mathcal L-geodesic. The scalar curvature contributes ddΟ„R(Ξ³(Ο„),Ο„)=RΟ„+βŸ¨βˆ‡R,X⟩\frac{d}{d\tau}R(\gamma(\tau), \tau) = R_\tau + \langle\nabla R, X\rangle. The speed contributes, by the new identity, ddΟ„βˆ£X∣2=2βŸ¨βˆ‡XX,X⟩+2Ric⁑(X,X)\frac{d}{d\tau}|X|^2 = 2\langle\nabla_XX, X\rangle + 2\operatorname{Ric}(X, X), and by (7.2), 2βŸ¨βˆ‡XX,X⟩=βŸ¨βˆ‡R,XβŸ©βˆ’1Ο„βˆ£X∣2βˆ’4Ric⁑(X,X)2\langle\nabla_XX, X\rangle = \langle\nabla R, X\rangle - \frac1\tau|X|^2 - 4\operatorname{Ric}(X, X). Adding, and adding and subtracting 1Ο„R\frac1\tau R,

ddΟ„(R+∣X∣2)=RΟ„+1Ο„R+2βŸ¨βˆ‡R,XβŸ©βˆ’2Ric⁑(X,X)βˆ’1Ο„(R+∣X∣2)=βˆ’H(X)βˆ’1Ο„(R+∣X∣2),(7.3)\frac{d}{d\tau}(R + |X|^2) = R_\tau + \frac1\tau R + 2\langle\nabla R, X\rangle - 2\operatorname{Ric}(X, X) - \frac1\tau(R + |X|^2) = -H(X) - \frac1\tau(R + |X|^2), \tag{7.3}

where

H(X)=βˆ’RΟ„βˆ’1Ο„Rβˆ’2βŸ¨βˆ‡R,X⟩+2Ric⁑(X,X).H(X) = -R_\tau - \frac1\tau R - 2\langle\nabla R, X\rangle + 2\operatorname{Ric}(X, X).

Perelman identifies H(X)H(X) as Hamilton's expression for the trace Harnack inequality with t=βˆ’Ο„t = -\tau. Indeed, Hamilton's trace Harnack quantity for a forward Ricci flow is Rt+Rt+2βŸ¨βˆ‡R,V⟩+2Ric⁑(V,V)R_t + \frac Rt + 2\langle\nabla R, V\rangle + 2\operatorname{Ric}(V, V) (11B.2 Ancient Solutions and the Harnack Inequality); substituting t=βˆ’Ο„t = -\tau (so Rt=βˆ’RΟ„R_t = -R_\tau and Rt=βˆ’RΟ„\frac Rt = -\frac R\tau) and V=βˆ’XV = -X (the velocity in forward time) gives exactly H(X)H(X). On a Ricci flow with nonnegative curvature operator defined for all earlier times, Hamilton's theorem says this expression is nonnegative for every VV when tt is measured from the initial time; here the base of Ο„\tau is the final time, which is why Perelman's Β§7.2 needs a correction term (below).

Integrating (7.3). Multiply by Ο„3/2\tau^{3/2}: ddΟ„(Ο„3/2(R+∣X∣2))=32Ο„(R+∣X∣2)+Ο„3/2ddΟ„(R+∣X∣2)=12Ο„(R+∣X∣2)βˆ’Ο„3/2H(X)\frac{d}{d\tau}\big(\tau^{3/2}(R + |X|^2)\big) = \frac32\sqrt\tau(R + |X|^2) + \tau^{3/2}\frac{d}{d\tau}(R + |X|^2) = \frac12\sqrt\tau(R + |X|^2) - \tau^{3/2}H(X). Integrate from 00 to Ο„Λ‰\bar\tau; the boundary term at 00 vanishes because Ο„3/2∣X∣2=Ο„βˆ£Ο„X∣2β†’0\tau^{3/2}|X|^2 = \sqrt\tau|\sqrt\tau X|^2 \to 0. With

K=K(Ξ³,Ο„Λ‰)=∫0Ο„Λ‰Ο„3/2H(X) dΟ„,K = K(\gamma, \bar\tau) = \int_0^{\bar\tau}\tau^{3/2}H(X)\,d\tau,

the result is

Ο„Λ‰3/2(R+∣X∣2)(Ο„Λ‰)=βˆ’K+12L(q,Ο„Λ‰).(7.4)\bar\tau^{3/2}\big(R + |X|^2\big)(\bar\tau) = -K + \frac12L(q, \bar\tau). \tag{7.4}

Substituting into the formulas for βˆ‚Ο„Λ‰L\partial_{\bar\tau}L and βˆ£βˆ‡L∣2|\nabla L|^2:

βˆ‚Ο„Λ‰L=2Ο„Λ‰Rβˆ’12Ο„Λ‰L+1Ο„Λ‰K,(7.5)\partial_{\bar\tau}L = 2\sqrt{\bar\tau}R - \frac{1}{2\bar\tau}L + \frac1{\bar\tau}K, \tag{7.5}
βˆ£βˆ‡L∣2=βˆ’4Ο„Λ‰R+2Ο„Λ‰Lβˆ’4Ο„Λ‰K.(7.6)|\nabla L|^2 = -4\bar\tau R + \frac{2}{\sqrt{\bar\tau}}L - \frac{4}{\sqrt{\bar\tau}}K. \tag{7.6}

In terms of the reduced distance β„“=L2Ο„Λ‰\ell = \frac{L}{2\sqrt{\bar\tau}} these become (Exercise 6.3)

βˆ‚Ο„Λ‰β„“=Rβˆ’β„“Ο„Λ‰+K2Ο„Λ‰3/2,βˆ£βˆ‡β„“βˆ£2=βˆ’R+β„“Ο„Λ‰βˆ’KΟ„Λ‰3/2.\partial_{\bar\tau}\ell = R - \frac{\ell}{\bar\tau} + \frac{K}{2\bar\tau^{3/2}}, \qquad |\nabla\ell|^2 = -R + \frac{\ell}{\bar\tau} - \frac{K}{\bar\tau^{3/2}}.

On flat space H≑0H \equiv 0, K=0K = 0, and these are βˆ‚Ο„β„“=βˆ’β„“Ο„\partial_\tau\ell = -\frac{\ell}{\tau} and βˆ£βˆ‡β„“βˆ£2=β„“Ο„|\nabla\ell|^2 = \frac\ell\tau, which β„“=∣x∣24Ο„\ell = \frac{|x|^2}{4\tau} satisfies. These are first-order identities: they involve KK, which depends on the whole geodesic, not just its endpoint. Second-order information comes next.

The second variation and L\mathcal L-Jacobi fields

Take a variation with Y(0)=0Y(0) = 0 and differentiate (7.1) once more. Perelman's (7.7) is the result after the same manipulations as in Riemannian geometry: commute βˆ‡Y\nabla_Y and βˆ‡X\nabla_X using the curvature tensor, move derivatives of YY along the curve using the new identity (which now produces derivatives of Ric⁑\operatorname{Ric}, since the metric's derivative is 2Ric⁑2\operatorname{Ric} and that changes along the curve), and discard terms proportional to the geodesic equation. The two ingredients are:

  • ddΟ„βŸ¨βˆ‡YY,X⟩=βŸ¨βˆ‡Xβˆ‡YY,X⟩+βŸ¨βˆ‡YY,βˆ‡XX⟩+2Yβ‹…Ric⁑(Y,X)βˆ’Xβ‹…Ric⁑(Y,Y)\frac{d}{d\tau}\langle\nabla_YY, X\rangle = \langle\nabla_X\nabla_YY, X\rangle + \langle\nabla_YY, \nabla_XX\rangle + 2Y\cdot\operatorname{Ric}(Y, X) - X\cdot\operatorname{Ric}(Y, Y), the identity for the derivative of βŸ¨βˆ‡YY,X⟩\langle\nabla_YY, X\rangle along the curve, in which the last two terms come from the Ο„\tau-derivative of the Christoffel symbols: βˆ‚Ο„Ξ“ijk=gkl(βˆ‡iRic⁑jl+βˆ‡jRic⁑ilβˆ’βˆ‡lRic⁑ij)\partial_\tau\Gamma^k_{ij} = g^{kl}(\nabla_i\operatorname{Ric}_{jl} + \nabla_j\operatorname{Ric}_{il} - \nabla_l\operatorname{Ric}_{ij}) (11A.2 How Curvature Evolves; Morgan–Tian, Claim 6.34, do the computation in coordinates);
  • the commutation βˆ‡Yβˆ‡XY=βˆ‡Xβˆ‡YY+R(Y,X)Y\nabla_Y\nabla_XY = \nabla_X\nabla_YY + R(Y, X)Y, with the curvature convention of Perelman's text.

The outcome, for any variation with Y(0)=0Y(0) = 0, is

Ξ΄Y2L=2βŸ¨βˆ‡YY,XβŸ©Ο„Λ‰+∫0Ο„Λ‰Ο„(βˆ‡Yβˆ‡YR+2⟨R(Y,X)Y,X⟩+2βˆ£βˆ‡XY∣2+2βˆ‡XRic⁑(Y,Y)βˆ’4βˆ‡YRic⁑(Y,X))dΟ„.(7.7)\delta^2_Y\mathcal L = 2\langle\nabla_YY, X\rangle\sqrt{\bar\tau} + \int_0^{\bar\tau}\sqrt\tau\Big(\nabla_Y\nabla_YR + 2\langle R(Y, X)Y, X\rangle + 2|\nabla_XY|^2 + 2\nabla_X\operatorname{Ric}(Y, Y) - 4\nabla_Y\operatorname{Ric}(Y, X)\Big)d\tau. \tag{7.7}

Compare the second variation of Riemannian energy, ∫(βˆ£βˆ‡XY∣2βˆ’βŸ¨R(Y,X)X,Y⟩)\int\big(|\nabla_XY|^2 - \langle R(Y, X)X, Y\rangle\big) (9A.7 Jacobi Fields and Curvature versus Topology): the curvature term and the βˆ£βˆ‡XY∣2|\nabla_XY|^2 term are the same up to factors; the Hessian of RR is the second variation of the potential; and the two derivatives of Ric⁑\operatorname{Ric} are new, the price of a moving metric.

L\mathcal L-Jacobi fields. The variation field of a family of L\mathcal L-geodesics satisfies the linearisation of (7.2), the L\mathcal L-Jacobi equation (Morgan–Tian, Lemma 6.10):

βˆ‡Xβˆ‡XY+R(Y,X)Xβˆ’12βˆ‡Y(βˆ‡R)+12Ο„βˆ‡XY+2(βˆ‡YRic⁑)(X,β‹…)+2Ric⁑(βˆ‡XY,β‹…)=0.\nabla_X\nabla_XY + R(Y, X)X - \frac12\nabla_Y(\nabla R) + \frac{1}{2\tau}\nabla_XY + 2(\nabla_Y\operatorname{Ric})(X, \cdot) + 2\operatorname{Ric}(\nabla_XY, \cdot) = 0.

As for Riemannian Jacobi fields: the second variation of a minimising L\mathcal L-geodesic is a nonnegative quadratic form on fields vanishing at both ends, it vanishes exactly on L\mathcal L-Jacobi fields (their Proposition 6.13), and the differential of Lexp⁑τˉ\mathcal L\exp_{\bar\tau} at vv sends ww to the value at Ο„Λ‰\bar\tau of the L\mathcal L-Jacobi field with Y(0)=0Y(0) = 0 and limβ‘Ο„βˆ‡XY=w\lim\sqrt\tau\nabla_XY = w (their Lemma 6.19). So Lexp⁑τˉ\mathcal L\exp_{\bar\tau} is singular at vv exactly when there is a nonzero L\mathcal L-Jacobi field vanishing at both ends: the L\mathcal L-conjugate points. An L\mathcal L-Jacobi field along a minimiser also computes the Hessian of LL exactly (their Lemma 6.40):

Hess⁑L(Y(Ο„Λ‰),Y(Ο„Λ‰))=2Ο„Λ‰βŸ¨βˆ‡XY,Y⟩(Ο„Λ‰).\operatorname{Hess}L(Y(\bar\tau), Y(\bar\tau)) = 2\sqrt{\bar\tau}\langle\nabla_XY, Y\rangle(\bar\tau).

The Hessian comparison with adapted fields

To bound Hess⁑L\operatorname{Hess}L from above, compare LL with the length of a family of curves that are not L\mathcal L-geodesics but are explicit, exactly as one bounds Hess⁑d\operatorname{Hess}d by using variations through linear fields. Let Ξ³\gamma be the minimiser to (q,Ο„Λ‰)(q, \bar\tau), and Y(Ο„Λ‰)Y(\bar\tau) a unit vector at qq. Any variation of Ξ³\gamma through curves from pp with variation field YY has length at least LL at its endpoints, with equality at u=0u = 0; so Hess⁑L(Y,Y)≀δY2Lβˆ’2Ο„Λ‰βŸ¨βˆ‡YY,X⟩\operatorname{Hess}L(Y, Y) \leq \delta^2_Y\mathcal L - 2\sqrt{\bar\tau}\langle\nabla_YY, X\rangle, the subtracted term being the correction for the endpoint moving along a curve rather than a geodesic of g(Ο„Λ‰)g(\bar\tau) (it cancels the boundary term of (7.7)). So

Hess⁑L(Y,Y)β‰€βˆ«0Ο„Λ‰Ο„(βˆ‡Yβˆ‡YR+2⟨R(Y,X)Y,X⟩+2βˆ£βˆ‡XY∣2+2βˆ‡XRic⁑(Y,Y)βˆ’4βˆ‡YRic⁑(Y,X))dΟ„\operatorname{Hess}L(Y, Y) \leq \int_0^{\bar\tau}\sqrt\tau\Big(\nabla_Y\nabla_YR + 2\langle R(Y, X)Y, X\rangle + 2|\nabla_XY|^2 + 2\nabla_X\operatorname{Ric}(Y, Y) - 4\nabla_Y\operatorname{Ric}(Y, X)\Big)d\tau

for every field YY along Ξ³\gamma with Y(0)=0Y(0) = 0 and the given value at Ο„Λ‰\bar\tau. Perelman now chooses YY to solve

βˆ‡XY=βˆ’Ric⁑(Y,β‹…)+12Ο„Y.(7.8)\nabla_XY = -\operatorname{Ric}(Y, \cdot) + \frac{1}{2\tau}Y. \tag{7.8}

These are adapted fields (Morgan–Tian's term). By the new identity, ddΟ„βˆ£Y∣2=2βŸ¨βˆ‡XY,Y⟩+2Ric⁑(Y,Y)=1Ο„βˆ£Y∣2\frac{d}{d\tau}|Y|^2 = 2\langle\nabla_XY, Y\rangle + 2\operatorname{Ric}(Y, Y) = \frac1\tau|Y|^2, so

∣Y(Ο„)∣2=Ο„Ο„Λ‰βˆ£Y(Ο„Λ‰)∣2,|Y(\tau)|^2 = \frac{\tau}{\bar\tau}|Y(\bar\tau)|^2,

and in particular Y(0)=0Y(0) = 0. In flat space they are Y(Ο„)=Ο„/τˉ EY(\tau) = \sqrt{\tau/\bar\tau}\,E with EE parallel: the variation fields of the fan of L\mathcal L-geodesics 2Ο„v2\sqrt\tau v, which are linear in s=Ο„s = \sqrt\tau. The choice is designed so that, in flat space, the comparison is an equality.

Substitute (7.8). The term 2βˆ£βˆ‡XY∣22|\nabla_XY|^2 becomes 2∣Ric⁑(Y,β‹…)∣2βˆ’2Ο„Ric⁑(Y,Y)+12Ο„2∣Y∣22|\operatorname{Ric}(Y, \cdot)|^2 - \frac2\tau\operatorname{Ric}(Y, Y) + \frac{1}{2\tau^2}|Y|^2, and with ∣Y∣2=ττˉ|Y|^2 = \frac\tau{\bar\tau} (normalising ∣Y(Ο„Λ‰)∣=1|Y(\bar\tau)| = 1), ∫0Ο„Λ‰Ο„12Ο„2ττˉdΟ„=1Ο„Λ‰\int_0^{\bar\tau}\sqrt\tau\frac{1}{2\tau^2}\frac{\tau}{\bar\tau}d\tau = \frac{1}{\sqrt{\bar\tau}}. The remaining terms are organised by the derivative of Ric⁑(Y,Y)\operatorname{Ric}(Y, Y) along the curve,

ddΟ„Ric⁑(Y,Y)=Ric⁑τ(Y,Y)+βˆ‡XRic⁑(Y,Y)+2Ric⁑(βˆ‡XY,Y)=Ric⁑τ(Y,Y)+βˆ‡XRic⁑(Y,Y)+1Ο„Ric⁑(Y,Y)βˆ’2∣Ric⁑(Y,β‹…)∣2,\frac{d}{d\tau}\operatorname{Ric}(Y, Y) = \operatorname{Ric}_\tau(Y, Y) + \nabla_X\operatorname{Ric}(Y, Y) + 2\operatorname{Ric}(\nabla_XY, Y) = \operatorname{Ric}_\tau(Y, Y) + \nabla_X\operatorname{Ric}(Y, Y) + \frac1\tau\operatorname{Ric}(Y, Y) - 2|\operatorname{Ric}(Y, \cdot)|^2,

which lets one replace βˆ‡XRic⁑(Y,Y)\nabla_X\operatorname{Ric}(Y, Y) by a total derivative and lower-order terms. Integrating the total derivative ddΟ„(2Ο„Ric⁑(Y,Y))\frac{d}{d\tau}(2\sqrt\tau\operatorname{Ric}(Y, Y)) gives the boundary term 2Ο„Λ‰Ric⁑(Y(Ο„Λ‰),Y(Ο„Λ‰))2\sqrt{\bar\tau}\operatorname{Ric}(Y(\bar\tau), Y(\bar\tau)) (Morgan–Tian carry this out as their equation (6.13)). The result is Perelman's (7.9):

Hess⁑L(Y,Y)≀1Ο„Λ‰βˆ’2Ο„Λ‰Ric⁑(Y,Y)βˆ’βˆ«0Ο„Λ‰Ο„H(X,Y) dΟ„,(7.9)\operatorname{Hess}L(Y, Y) \leq \frac{1}{\sqrt{\bar\tau}} - 2\sqrt{\bar\tau}\operatorname{Ric}(Y, Y) - \int_0^{\bar\tau}\sqrt\tau H(X, Y)\,d\tau, \tag{7.9}

where

H(X,Y)=βˆ’βˆ‡Yβˆ‡YRβˆ’2⟨R(Y,X)Y,XβŸ©βˆ’4(βˆ‡XRic⁑(Y,Y)βˆ’βˆ‡YRic⁑(Y,X))βˆ’2Ric⁑τ(Y,Y)+2∣Ric⁑(Y,β‹…)∣2βˆ’1Ο„Ric⁑(Y,Y)H(X, Y) = -\nabla_Y\nabla_YR - 2\langle R(Y, X)Y, X\rangle - 4\big(\nabla_X\operatorname{Ric}(Y, Y) - \nabla_Y\operatorname{Ric}(Y, X)\big) - 2\operatorname{Ric}_\tau(Y, Y) + 2|\operatorname{Ric}(Y, \cdot)|^2 - \frac1\tau\operatorname{Ric}(Y, Y)

is Hamilton's expression for the matrix Harnack inequality, at t=βˆ’Ο„t = -\tau. It is quadratic in YY (despite the notation, it depends on XX and YY). Equality in (7.9) holds exactly when the adapted field YY is also an L\mathcal L-Jacobi field (Morgan–Tian, Proposition 6.37).

Compare once more with Riemannian geometry. For the distance function rr from a point in a manifold with Ric⁑β‰₯0\operatorname{Ric} \geq 0, the Hessian comparison along a unit-speed geodesic with the linear fields srE\frac srE gives Hess⁑r(E,E)≀1rβˆ’1r2∫0rs2Rm⁑(E,Ξ³Λ™,Ξ³Λ™,E) ds\operatorname{Hess}r(E, E) \leq \frac1r - \frac1{r^2}\int_0^rs^2\operatorname{Rm}(E, \dot\gamma, \dot\gamma, E)\,ds, and tracing gives Ξ”r≀nβˆ’1r\Delta r \leq \frac{n - 1}{r} (9B.1 Laplacian Comparison). Here the role of the linear fields is taken by the adapted fields, the role of the sectional curvature integral by the integral of the matrix Harnack expression, and the curvature hypothesis is replaced by nothing: the Harnack expression is simply carried along.

The trace: the Laplacian comparison

Choose an orthonormal basis {YΞ±(Ο„Λ‰)}\{Y_\alpha(\bar\tau)\} at qq and extend each by (7.8). Since (7.8) is linear and ddΟ„βŸ¨YΞ±,Yβ⟩=1Ο„βŸ¨YΞ±,Yβ⟩\frac{d}{d\tau}\langle Y_\alpha, Y_\beta\rangle = \frac1\tau\langle Y_\alpha, Y_\beta\rangle, the fields Ο„Λ‰/τ YΞ±(Ο„)\sqrt{\bar\tau/\tau}\,Y_\alpha(\tau) are orthonormal at every Ο„\tau. Summing (7.9) over Ξ±\alpha:

Ξ”L≀nΟ„Λ‰βˆ’2Ο„Λ‰Rβˆ’βˆ«0Ο„Λ‰Ο„βˆ‘Ξ±H(X,YΞ±) dΟ„.\Delta L \leq \frac{n}{\sqrt{\bar\tau}} - 2\sqrt{\bar\tau}R - \int_0^{\bar\tau}\sqrt\tau\sum_\alpha H(X, Y_\alpha)\,d\tau.

The trace identity says that for an orthonormal basis {WΞ±}\{W_\alpha\},

βˆ‘Ξ±H(X,WΞ±)=H(X),\sum_\alpha H(X, W_\alpha) = H(X),

the matrix Harnack expression traces to the trace Harnack expression (Exercise 6.4; Morgan–Tian, Claim 6.44). It uses three facts: the contracted second Bianchi identity βˆ‘Ξ±(βˆ‡WΞ±Ric⁑)(WΞ±,X)=12βŸ¨βˆ‡R,X⟩\sum_\alpha(\nabla_{W_\alpha}\operatorname{Ric})(W_\alpha, X) = \frac12\langle\nabla R, X\rangle; the evolution RΟ„=βˆ’Ξ”Rβˆ’2∣Ric⁑∣2R_\tau = -\Delta R - 2|\operatorname{Ric}|^2 of the scalar curvature in backward time (11A.2 How Curvature Evolves); and RΟ„=βˆ‚Ο„(gijRij)=βˆ’2∣Ric⁑∣2+βˆ‘Ξ±Ric⁑τ(WΞ±,WΞ±)R_\tau = \partial_\tau(g^{ij}R_{ij}) = -2|\operatorname{Ric}|^2 + \sum_\alpha\operatorname{Ric}_\tau(W_\alpha, W_\alpha), so that βˆ‘Ξ±Ric⁑τ(WΞ±,WΞ±)=βˆ’Ξ”R\sum_\alpha\operatorname{Ric}_\tau(W_\alpha, W_\alpha) = -\Delta R. Since ∣Yα∣2=ττˉ|Y_\alpha|^2 = \frac\tau{\bar\tau}, the sum is ττˉH(X)\frac{\tau}{\bar\tau}H(X), and the integral is 1Ο„Λ‰βˆ«0Ο„Λ‰Ο„3/2H(X) dΟ„=KΟ„Λ‰\frac1{\bar\tau}\int_0^{\bar\tau}\tau^{3/2}H(X)\,d\tau = \frac{K}{\bar\tau}. So

Ξ”Lβ‰€βˆ’2Ο„Λ‰R+nΟ„Λ‰βˆ’1Ο„Λ‰K.(7.10)\Delta L \leq -2\sqrt{\bar\tau}R + \frac{n}{\sqrt{\bar\tau}} - \frac1{\bar\tau}K. \tag{7.10}

This is the space-time Laplacian comparison, the analogue of Ξ”r2≀2n\Delta r^2 \leq 2n in a space of nonnegative Ricci curvature. In terms of β„“\ell, Ξ”β„“=Ξ”L2Ο„Λ‰β‰€βˆ’R+n2Ο„Λ‰βˆ’K2Ο„Λ‰3/2\Delta\ell = \frac{\Delta L}{2\sqrt{\bar\tau}} \leq -R + \frac{n}{2\bar\tau} - \frac{K}{2\bar\tau^{3/2}}, which is the third inequality quoted in 12A.5 Reduced Distance and Reduced Volume.

Combining: the inequalities (7.13)–(7.15)

Now the unknown KK can be eliminated. With the two identities for βˆ‚Ο„β„“\partial_\tau\ell and βˆ£βˆ‡β„“βˆ£2|\nabla\ell|^2 and the inequality for Ξ”β„“\Delta\ell:

βˆ‚Ο„β„“βˆ’Ξ”β„“+βˆ£βˆ‡β„“βˆ£2βˆ’R+n2Ο„β‰₯(Rβˆ’β„“Ο„+K2Ο„3/2)βˆ’(βˆ’R+n2Ο„βˆ’K2Ο„3/2)+(βˆ’R+β„“Ο„βˆ’KΟ„3/2)βˆ’R+n2Ο„=0.(7.13)\partial_\tau\ell - \Delta\ell + |\nabla\ell|^2 - R + \frac{n}{2\tau} \geq \Big(R - \frac\ell\tau + \frac{K}{2\tau^{3/2}}\Big) - \Big(-R + \frac{n}{2\tau} - \frac{K}{2\tau^{3/2}}\Big) + \Big(-R + \frac\ell\tau - \frac{K}{\tau^{3/2}}\Big) - R + \frac{n}{2\tau} = 0. \tag{7.13}

The KK terms cancel: K2+K2βˆ’K=0\frac{K}{2} + \frac K2 - K = 0. Similarly

2Ξ”β„“βˆ’βˆ£βˆ‡β„“βˆ£2+R+β„“βˆ’nτ≀2(βˆ’R+n2Ο„βˆ’K2Ο„3/2)βˆ’(βˆ’R+β„“Ο„βˆ’KΟ„3/2)+R+β„“βˆ’nΟ„=0.(7.14)2\Delta\ell - |\nabla\ell|^2 + R + \frac{\ell - n}{\tau} \leq 2\Big(-R + \frac{n}{2\tau} - \frac{K}{2\tau^{3/2}}\Big) - \Big(-R + \frac\ell\tau - \frac{K}{\tau^{3/2}}\Big) + R + \frac{\ell - n}{\tau} = 0. \tag{7.14}

And for Lˉ=2τL=4τℓ\bar L = 2\sqrt\tau L = 4\tau\ell, from (7.5) and (7.10) (12A.5 Reduced Distance and Reduced Volume, exercise on the minimum of ℓ\ell),

βˆ‚Ο„LΛ‰+Ξ”Lˉ≀2n.(7.15)\partial_\tau\bar L + \Delta\bar L \leq 2n. \tag{7.15}

In Morgan and Tian's formulation (their Corollary 6.51), if one sets Ξ΄=n2Ο„βˆ’Rβˆ’K2Ο„3/2βˆ’Ξ”β„“β‰₯0\delta = \frac{n}{2\tau} - R - \frac{K}{2\tau^{3/2}} - \Delta\ell \geq 0, the "defect" in the Laplacian comparison, then (7.13) and (7.14) are the equalities βˆ‚Ο„β„“βˆ’Ξ”β„“+βˆ£βˆ‡β„“βˆ£2βˆ’R+n2Ο„=Ξ΄\partial_\tau\ell - \Delta\ell + |\nabla\ell|^2 - R + \frac{n}{2\tau} = \delta and 2Ξ”β„“βˆ’βˆ£βˆ‡β„“βˆ£2+R+β„“βˆ’nΟ„=βˆ’2Ξ΄2\Delta\ell - |\nabla\ell|^2 + R + \frac{\ell - n}{\tau} = -2\delta. All the inequality in the theory comes from the single Hessian comparison.

The equality case. Equality in (7.13) at a point forces equality in (7.10), hence in (7.9) for every unit YY, hence every adapted field is an L\mathcal L-Jacobi field, and then (Morgan–Tian, Proposition 6.43)

Ric⁑+12Ο„Hess⁑L=12Ο„g,thatΒ isRic⁑+βˆ‡2β„“=12Ο„g:\operatorname{Ric} + \frac{1}{2\sqrt\tau}\operatorname{Hess}L = \frac{1}{2\tau}g, \qquad\text{that is}\qquad \operatorname{Ric} + \nabla^2\ell = \frac{1}{2\tau}g:

the gradient shrinking soliton equation with potential β„“\ell. So if equality holds everywhere, the flow is a gradient shrinking soliton, and the Gaussian soliton on Rn\mathbb{R}^n with β„“=∣x∣24Ο„\ell = \frac{|x|^2}{4\tau} is the model.

In the world Analogy Why the Harnack expressions appear

In Riemannian geometry, the second variation of length is controlled by curvature, and curvature hypotheses (such as Ric⁑β‰₯0\operatorname{Ric} \geq 0) turn it into comparison theorems. For the Ricci flow, the analogous role is played by Hamilton's Harnack expressions, which are combinations of curvature and its space and time derivatives. Perelman's heuristic of Β§6 (12A.5 Reduced Distance and Reduced Volume) explains why: on the high-dimensional space-time MΓ—SNΓ—R+M\times S^N\times\mathbb{R}_+ with the metric g(Ο„)+Ο„gSN+(N2Ο„+R)dΟ„2g(\tau) + \tau g_{S^N} + (\frac{N}{2\tau} + R)d\tau^2, the curvature components, to leading order in 1N\frac1N, are the Harnack expressions, and the space is Ricci-flat to that order. The L\mathcal L-geodesics are its geodesics, the second variation of L\mathcal L is its second variation of length, and the Harnack expressions appear because they are its curvature. The analogy is exact in the limit and is the reason the computations work; it is not a proof, and the proof is the computation of this chapter.

When the curvature operator is nonnegative

Perelman's Β§7.2 adds one more estimate, used for ΞΊ\kappa-solutions (12B.3 The Structure of ΞΊ-Solutions). Suppose the backward flow, defined for Ο„βˆˆ[0,Ο„0]\tau \in [0, \tau_0], has nonnegative curvature operator. Then Hamilton's Harnack inequality holds, with time measured from the initial time of the forward flow, tt-time Tβˆ’Ο„0T - \tau_0; translated into Ο„\tau, it gives

H(X,Y)β‰₯βˆ’Ric⁑(Y,Y)(1Ο„+1Ο„0βˆ’Ο„)β‰₯βˆ’R(1Ο„+1Ο„0βˆ’Ο„)∣Y∣2,H(X)β‰₯βˆ’R(1Ο„+1Ο„0βˆ’Ο„).H(X, Y) \geq -\operatorname{Ric}(Y, Y)\Big(\frac1\tau + \frac{1}{\tau_0 - \tau}\Big) \geq -R\Big(\frac1\tau + \frac{1}{\tau_0 - \tau}\Big)|Y|^2, \qquad H(X) \geq -R\Big(\frac1\tau + \frac{1}{\tau_0 - \tau}\Big).

The correction 1Ο„0βˆ’Ο„\frac{1}{\tau_0 - \tau} is the price of measuring from the final rather than the initial time; it is harmless as long as τ≀(1βˆ’c)Ο„0\tau \leq (1 - c)\tau_0. With these lower bounds on HH, the identity for βˆ£βˆ‡β„“βˆ£2|\nabla\ell|^2 and (7.11) give, for τ≀(1βˆ’c)Ο„0\tau \leq (1 - c)\tau_0,

βˆ£βˆ‡β„“βˆ£2+R≀Cβ„“Ο„,ddΟ„log⁑∣Y∣2≀1Ο„(Cβ„“+1)(7.16,Β 7.17)|\nabla\ell|^2 + R \leq \frac{C\ell}{\tau}, \qquad \frac{d}{d\tau}\log|Y|^2 \leq \frac1\tau(C\ell + 1) \tag{7.16, 7.17}

for L\mathcal L-Jacobi fields YY, with CC depending on cc (statement; the derivation uses R≀R+∣X∣2R \leq R + |X|^2 and (7.4) to bound βˆ«Ο„3/2R\int\tau^{3/2}R in terms of LL). For an ancient solution, Ο„0=∞\tau_0 = \infty and the correction disappears entirely. These bounds say that on a ΞΊ\kappa-solution, β„“\ell controls both its own gradient and the curvature, at the scale Ο„\sqrt\tau, which is what makes the rescaled reduced distances compact as Ο„β†’βˆž\tau \to \infty and produces the asymptotic soliton (12B.1 ΞΊ-Solutions, 12B.3 The Structure of ΞΊ-Solutions).

Where this goes Where this goes

12A.7 Monotonicity of the Reduced Volume makes the computations rigorous: it studies the L\mathcal L-exponential map, its Jacobian and its injectivity domain, proves that Ο„βˆ’n/2eβˆ’β„“J\tau^{-n/2}e^{-\ell}J is nonincreasing along each L\mathcal L-geodesic (Perelman's (7.12)), and turns (7.13) and (7.15) into barrier inequalities on complete flows, with min⁑ℓ≀n2\min\ell \leq \frac n2 as a consequence. 12A.8 Local Noncollapsing from the Reduced Volume uses the result for local noncollapsing. 12B.1 ΞΊ-Solutions and 12B.3 The Structure of ΞΊ-Solutions use (7.16) for ΞΊ\kappa-solutions.

History

Perelman's Β§7.1 appeared in November 2002. He named Li and Yau's 1986 paper as the closest reference: Li and Yau integrated their gradient estimate along paths using a "length" associated with a linear parabolic equation, the same structure with a potential. Hamilton's trace and matrix Harnack inequalities for the Ricci flow date from 1993, and Bennett Chow and Sun-Chin Chu's geometric interpretation of them on space-time from 1995. The first complete written accounts of Β§7.1 are in Kleiner and Lott's notes (posted from 2003, published 2008), Morgan and Tian's book (2007) and Cao and Zhu's article (2006); Rugang Ye's notes on the reduced distance and reduced volume (2004) also treat it.

Recall Where we stand

In s=Ο„s = \sqrt\tau the L\mathcal L-length is the energy ∫(12∣A∣2+2s2R) ds\int(\frac12|A|^2 + 2s^2R)\,ds, and L\mathcal L-geodesics are regular at Ο„=0\tau = 0, determined by lim⁑τX\lim\sqrt\tau X. The only new rule is ddΟ„βŸ¨V,W⟩=βŸ¨βˆ‡XV,W⟩+⟨V,βˆ‡XW⟩+2Ric⁑(V,W)\frac{d}{d\tau}\langle V, W\rangle = \langle\nabla_XV, W\rangle + \langle V, \nabla_XW\rangle + 2\operatorname{Ric}(V, W). It gives the first variation (7.1), the geodesic equation (7.2), βˆ‡L=2Ο„X\nabla L = 2\sqrt\tau X and βˆ‚Ο„L=2Ο„Rβˆ’Ο„(R+∣X∣2)\partial_\tau L = 2\sqrt\tau R - \sqrt\tau(R + |X|^2). Along a geodesic ddΟ„(R+∣X∣2)=βˆ’H(X)βˆ’1Ο„(R+∣X∣2)\frac{d}{d\tau}(R + |X|^2) = -H(X) - \frac1\tau(R + |X|^2) with HH Hamilton's trace Harnack expression, so Ο„3/2(R+∣X∣2)=12Lβˆ’K\tau^{3/2}(R + |X|^2) = \frac12L - K and (7.5)–(7.6). The second variation (7.7), with adapted fields βˆ‡XY=βˆ’Ric⁑(Y)+12Ο„Y\nabla_XY = -\operatorname{Ric}(Y) + \frac{1}{2\tau}Y (so ∣Y∣2=ττˉ|Y|^2 = \frac\tau{\bar\tau}), gives the Hessian comparison (7.9) with Hamilton's matrix Harnack expression, which traces to H(X)H(X) and gives Ξ”Lβ‰€βˆ’2Ο„R+nΟ„βˆ’KΟ„\Delta L \leq -2\sqrt\tau R + \frac{n}{\sqrt\tau} - \frac K\tau (7.10). Eliminating KK gives (7.13)–(7.15), with equality only on gradient shrinking solitons. Under nonnegative curvature operator, Hamilton's Harnack gives βˆ£βˆ‡β„“βˆ£2+R≀Cβ„“Ο„|\nabla\ell|^2 + R \leq \frac{C\ell}{\tau}. 12A.7 Monotonicity of the Reduced Volume turns this into the monotonicity of the reduced volume.

Exercises

Exercise 6.1 The geodesic equation in ss

With s=Ο„s = \sqrt\tau, A=dΞ³ds=2sXA = \frac{d\gamma}{ds} = 2sX and ddΟ„=12sdds\frac{d}{d\tau} = \frac{1}{2s}\frac{d}{ds}: (a) show βˆ‡XX=14s2βˆ‡AAβˆ’14s3A\nabla_XX = \frac{1}{4s^2}\nabla_AA - \frac{1}{4s^3}A; (b) deduce that (7.2) is equivalent to βˆ‡AAβˆ’2s2βˆ‡R+4sRic⁑(A,β‹…)=0\nabla_AA - 2s^2\nabla R + 4s\operatorname{Ric}(A, \cdot) = 0; (c) on static flat space, solve it and recover Ξ³(Ο„)=p+2Ο„v\gamma(\tau) = p + 2\sqrt\tau v.

Solution

(a) βˆ‡XX=12sβˆ‡A(A2s)=12s(12sβˆ‡AAβˆ’12s2A)\nabla_XX = \frac{1}{2s}\nabla_A\big(\frac{A}{2s}\big) = \frac{1}{2s}\big(\frac{1}{2s}\nabla_AA - \frac{1}{2s^2}A\big). (b) Substituting into (7.2): 14s2βˆ‡AAβˆ’14s3Aβˆ’12βˆ‡R+12s2β‹…A2s+2Ric⁑(A2s,β‹…)=0\frac{1}{4s^2}\nabla_AA - \frac{1}{4s^3}A - \frac12\nabla R + \frac{1}{2s^2}\cdot\frac{A}{2s} + 2\operatorname{Ric}\big(\frac{A}{2s}, \cdot\big) = 0. The terms βˆ“14s3A\mp\frac{1}{4s^3}A cancel, and multiplying by 4s24s^2 gives the stated equation. (c) With R=0R = 0 and Ric⁑=0\operatorname{Ric} = 0, d2Ξ³ds2=0\frac{d^2\gamma}{ds^2} = 0, so Ξ³=p+sA0\gamma = p + sA_0, and A0=2vA_0 = 2v with v=lim⁑τXv = \lim\sqrt\tau X.

Exercise 6.2 The first variation on flat space

On static flat Rn\mathbb{R}^n, compute Ξ΄YL\delta_Y\mathcal L directly for Ξ³(Ο„)=2Ο„v\gamma(\tau) = 2\sqrt\tau v and a variation field YY with Y(0)=0Y(0) = 0, and check that it equals 2Ο„Λ‰βŸ¨X(Ο„Λ‰),Y(Ο„Λ‰)⟩2\sqrt{\bar\tau}\langle X(\bar\tau), Y(\bar\tau)\rangle. Deduce βˆ‡L(x,Ο„Λ‰)=xΟ„Λ‰\nabla L(x, \bar\tau) = \frac{x}{\sqrt{\bar\tau}}, and check against L=∣x∣22Ο„Λ‰L = \frac{|x|^2}{2\sqrt{\bar\tau}}.

Solution

X=vΟ„βˆ’1/2X = v\tau^{-1/2}, and the integrand of (7.1) is Ο„βŸ¨Y,βˆ’2βˆ‡XXβˆ’1Ο„X⟩\sqrt\tau\langle Y, -2\nabla_XX - \frac1\tau X\rangle with βˆ‡XX=XΛ™=βˆ’12vΟ„βˆ’3/2\nabla_XX = \dot X = -\frac12v\tau^{-3/2}, so βˆ’2XΛ™βˆ’1Ο„X=vΟ„βˆ’3/2βˆ’vΟ„βˆ’3/2=0-2\dot X - \frac1\tau X = v\tau^{-3/2} - v\tau^{-3/2} = 0. Only the boundary term remains: 2Ο„Λ‰βŸ¨X(Ο„Λ‰),Y(Ο„Λ‰)⟩2\sqrt{\bar\tau}\langle X(\bar\tau), Y(\bar\tau)\rangle, since at 00, 2Ο„βŸ¨X,Y⟩=2⟨v,Y(0)⟩=02\sqrt\tau\langle X, Y\rangle = 2\langle v, Y(0)\rangle = 0. So βˆ‡L=2Ο„Λ‰X(Ο„Λ‰)=2v\nabla L = 2\sqrt{\bar\tau}X(\bar\tau) = 2v, and with x=2Ο„Λ‰vx = 2\sqrt{\bar\tau}v, βˆ‡L=xΟ„Λ‰\nabla L = \frac{x}{\sqrt{\bar\tau}}. Directly, βˆ‡βˆ£x∣22Ο„Λ‰=xΟ„Λ‰\nabla\frac{|x|^2}{2\sqrt{\bar\tau}} = \frac{x}{\sqrt{\bar\tau}}.

Exercise 6.3 From LL to β„“\ell

From (7.5) and (7.6), derive βˆ‚Ο„β„“=Rβˆ’β„“Ο„+K2Ο„3/2\partial_\tau\ell = R - \frac\ell\tau + \frac{K}{2\tau^{3/2}} and βˆ£βˆ‡β„“βˆ£2=βˆ’R+β„“Ο„βˆ’KΟ„3/2|\nabla\ell|^2 = -R + \frac\ell\tau - \frac{K}{\tau^{3/2}}, where β„“=L2Ο„\ell = \frac{L}{2\sqrt\tau}.

Solution

βˆ‚Ο„β„“=βˆ‚Ο„L2Ο„βˆ’L4Ο„3/2=Rβˆ’L4Ο„3/2+K2Ο„3/2βˆ’L4Ο„3/2=Rβˆ’L2Ο„3/2+K2Ο„3/2\partial_\tau\ell = \frac{\partial_\tau L}{2\sqrt\tau} - \frac{L}{4\tau^{3/2}} = R - \frac{L}{4\tau^{3/2}} + \frac{K}{2\tau^{3/2}} - \frac{L}{4\tau^{3/2}} = R - \frac{L}{2\tau^{3/2}} + \frac{K}{2\tau^{3/2}}, and L2Ο„3/2=β„“Ο„\frac{L}{2\tau^{3/2}} = \frac\ell\tau. Also βˆ£βˆ‡β„“βˆ£2=βˆ£βˆ‡L∣24Ο„=βˆ’R+L2Ο„3/2βˆ’KΟ„3/2=βˆ’R+β„“Ο„βˆ’KΟ„3/2|\nabla\ell|^2 = \frac{|\nabla L|^2}{4\tau} = -R + \frac{L}{2\tau^{3/2}} - \frac{K}{\tau^{3/2}} = -R + \frac\ell\tau - \frac{K}{\tau^{3/2}}.

Exercise 6.4 The trace identity

Using the conventions of the text (in particular βˆ‘Ξ±βŸ¨R(WΞ±,X)WΞ±,X⟩=βˆ’Ric⁑(X,X)\sum_\alpha\langle R(W_\alpha, X)W_\alpha, X\rangle = -\operatorname{Ric}(X, X) for an orthonormal basis), the contracted Bianchi identity, RΟ„=βˆ’Ξ”Rβˆ’2∣Ric⁑∣2R_\tau = -\Delta R - 2|\operatorname{Ric}|^2 and βˆ‘Ξ±Ric⁑τ(WΞ±,WΞ±)=βˆ’Ξ”R\sum_\alpha\operatorname{Ric}_\tau(W_\alpha, W_\alpha) = -\Delta R, show that βˆ‘Ξ±H(X,WΞ±)=H(X)\sum_\alpha H(X, W_\alpha) = H(X).

Solution

Sum the terms of H(X,WΞ±)H(X, W_\alpha): βˆ’βˆ‘βˆ‡WΞ±βˆ‡WΞ±R=βˆ’Ξ”R-\sum\nabla_{W_\alpha}\nabla_{W_\alpha}R = -\Delta R; βˆ’2βˆ‘βŸ¨R(WΞ±,X)WΞ±,X⟩=2Ric⁑(X,X)-2\sum\langle R(W_\alpha, X)W_\alpha, X\rangle = 2\operatorname{Ric}(X, X); βˆ’4βˆ‘βˆ‡XRic⁑(WΞ±,WΞ±)=βˆ’4βŸ¨βˆ‡R,X⟩-4\sum\nabla_X\operatorname{Ric}(W_\alpha, W_\alpha) = -4\langle\nabla R, X\rangle (the trace of βˆ‡XRic⁑\nabla_X\operatorname{Ric} is βˆ‡XR\nabla_XR); +4βˆ‘βˆ‡WΞ±Ric⁑(WΞ±,X)=2βŸ¨βˆ‡R,X⟩+4\sum\nabla_{W_\alpha}\operatorname{Ric}(W_\alpha, X) = 2\langle\nabla R, X\rangle (Bianchi); βˆ’2βˆ‘Ric⁑τ(WΞ±,WΞ±)=2Ξ”R-2\sum\operatorname{Ric}_\tau(W_\alpha, W_\alpha) = 2\Delta R; +2βˆ‘βˆ£Ric⁑(WΞ±,β‹…)∣2=2∣Ric⁑∣2+2\sum|\operatorname{Ric}(W_\alpha, \cdot)|^2 = 2|\operatorname{Ric}|^2; βˆ’1Ο„βˆ‘Ric⁑(WΞ±,WΞ±)=βˆ’RΟ„-\frac1\tau\sum\operatorname{Ric}(W_\alpha, W_\alpha) = -\frac R\tau. Total: Ξ”R+2∣Ric⁑∣2+2Ric⁑(X,X)βˆ’2βŸ¨βˆ‡R,XβŸ©βˆ’RΟ„=βˆ’RΟ„βˆ’RΟ„βˆ’2βŸ¨βˆ‡R,X⟩+2Ric⁑(X,X)=H(X)\Delta R + 2|\operatorname{Ric}|^2 + 2\operatorname{Ric}(X, X) - 2\langle\nabla R, X\rangle - \frac R\tau = -R_\tau - \frac R\tau - 2\langle\nabla R, X\rangle + 2\operatorname{Ric}(X, X) = H(X).

Exercise 6.5 Adapted fields in flat space

On static flat space along Ξ³(Ο„)=2Ο„v\gamma(\tau) = 2\sqrt\tau v, show that the adapted fields are Y(Ο„)=Ο„/τˉ EY(\tau) = \sqrt{\tau/\bar\tau}\,E with EE constant, that they are L\mathcal L-Jacobi fields, and that (7.9) is an equality: Hess⁑L(E,E)=1Ο„Λ‰βˆ£E∣2\operatorname{Hess}L(E, E) = \frac{1}{\sqrt{\bar\tau}}|E|^2. Check against L=∣x∣22Ο„Λ‰L = \frac{|x|^2}{2\sqrt{\bar\tau}}.

Solution

(7.8) with Ric⁑=0\operatorname{Ric} = 0 is YΛ™=12Ο„Y\dot Y = \frac{1}{2\tau}Y, so Y=cΟ„EY = c\sqrt\tau E, and Y(Ο„Λ‰)=EY(\bar\tau) = E gives c=Ο„Λ‰βˆ’1/2c = \bar\tau^{-1/2}. The variation Ξ³u(Ο„)=2Ο„(v+uE/(2Ο„Λ‰))\gamma_u(\tau) = 2\sqrt\tau(v + uE/(2\sqrt{\bar\tau})) is a family of L\mathcal L-geodesics with variation field Ο„/Ο„Λ‰E\sqrt{\tau/\bar\tau}E, so YY is L\mathcal L-Jacobi. All curvature terms vanish and H=0H = 0, so (7.9) reads Hess⁑L≀1Ο„Λ‰\operatorname{Hess}L \leq \frac{1}{\sqrt{\bar\tau}}, with equality since YY is L\mathcal L-Jacobi. Directly, Hess⁑∣x∣22Ο„Λ‰=1Ο„Λ‰I\operatorname{Hess}\frac{|x|^2}{2\sqrt{\bar\tau}} = \frac{1}{\sqrt{\bar\tau}}I.

Exercise 6.6 The Harnack expression on a shrinking sphere

For the round sphere in backward time, g(Ο„)=(1+2(nβˆ’1)Ο„)gSng(\tau) = (1 + 2(n - 1)\tau)g_{S^n} (the unit sphere at Ο„=0\tau = 0; in forward time this is the ancient shrinking sphere), compute R(Ο„)R(\tau) and H(X)H(X) for the constant path Ξ³(Ο„)=p\gamma(\tau) = p, and show that H(0)H(0) is negative. Show that nevertheless H(0)β‰₯βˆ’RΟ„H(0) \geq -\frac R\tau, which is Perelman's lower bound of Β§7.2 with Ο„0=∞\tau_0 = \infty, and explain how it follows from Hamilton's Harnack inequality for this ancient solution.

Solution

βˆ‚Ο„g=2Ric⁑=2(nβˆ’1)gSn\partial_\tau g = 2\operatorname{Ric} = 2(n - 1)g_{S^n}, so the factor is 1+2(nβˆ’1)Ο„1 + 2(n - 1)\tau and R=n(nβˆ’1)1+2(nβˆ’1)Ο„R = \frac{n(n - 1)}{1 + 2(n - 1)\tau}, RΟ„=βˆ’2n(nβˆ’1)2(1+2(nβˆ’1)Ο„)2R_\tau = -\frac{2n(n - 1)^2}{(1 + 2(n - 1)\tau)^2}. With X=0X = 0, H(0)=βˆ’RΟ„βˆ’RΟ„=n(nβˆ’1)1+2(nβˆ’1)Ο„(2(nβˆ’1)1+2(nβˆ’1)Ο„βˆ’1Ο„)=βˆ’n(nβˆ’1)Ο„(1+2(nβˆ’1)Ο„)2<0H(0) = -R_\tau - \frac R\tau = \frac{n(n - 1)}{1 + 2(n - 1)\tau}\Big(\frac{2(n - 1)}{1 + 2(n - 1)\tau} - \frac1\tau\Big) = -\frac{n(n - 1)}{\tau(1 + 2(n - 1)\tau)^2} < 0. But H(0)+RΟ„=βˆ’RΟ„>0H(0) + \frac R\tau = -R_\tau > 0. For an ancient solution with nonnegative curvature operator, Hamilton's trace Harnack inequality (with the initial time sent to βˆ’βˆž-\infty) gives Rtβ‰₯0R_t \geq 0, that is βˆ’RΟ„β‰₯0-R_\tau \geq 0, which is exactly H(0)β‰₯βˆ’RΟ„H(0) \geq -\frac R\tau. The negative term βˆ’RΟ„-\frac R\tau in HH is there because HH measures time from the base point Ο„=0\tau = 0, not from the beginning of the flow.

Free resources for this chapter

Chosen for what this chapter covers, and all free to read. Links checked Oct 16, 2026.

Β© 2026 NeckPinch (www.neckpinch.com). All content in the guidebook (text, mathematics, figures and exercises) is protected by copyright. All rights reserved. No part may be copied, republished or redistributed without written permission.