Book 8A

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

Course 8Book 8A: Smooth ManifoldsChapter 6

Flows and the Lie Derivative

Vector fields, their flows, brackets, and the diffeomorphism invariance behind solitons and DeTurck’s trick.

22 min read · Updated Oct 3, 2026

Read with Lee, Introduction to Smooth Manifolds (2nd edition), chapter 8 (vector fields and the Lie bracket) and chapter 9 (integral curves and flows, Lie derivatives, commuting vector fields, time-dependent vector fields). Lee's ODE proofs in chapter 9 repeat 2B.10 Ordinary Differential Equations and can be skipped; keep the statements.

In this chapter · 7 sections
  1. 6.1How a fluid strains
  2. 6.2Integral curves and flows
  3. 6.3The Lie bracket
  4. 6.4The Lie derivative
  5. 6.5Naturality and its consequences
  6. 6.6History
  7. 6.7Exercises

A vector field on a manifold is a velocity at every point. Following those velocities moves every point at once: the flow of the vector field, a one-parameter family of diffeomorphisms. Flows are how vector fields act on everything else. Differentiating any geometric object along a flow gives its Lie derivative, the rate at which the object changes as it is dragged along. For a function this is just the directional derivative; for another vector field it is the Lie bracket, which measures how two flows fail to commute; and for a metric it is the strain: how the flow stretches lengths.

This chapter ends with the single most important structural fact about the Ricci flow. The Ricci tensor is natural: pulling a metric back by a diffeomorphism pulls back its Ricci tensor. So if g(t)g(t) is a Ricci flow, so is ϕ∗g(t)\phi^*g(t) for any fixed diffeomorphism ϕ\phi. That freedom, the diffeomorphism invariance of the equation, has two consequences that run through the rest of the guide: the equation cannot be strictly parabolic, which DeTurck's trick repairs (11A.3 Short-Time Existence and Uniqueness); and there are solutions that move only by diffeomorphisms and scaling, the Ricci solitons, whose defining equation Ric⁡+12LXg=λg\operatorname{Ric} + \frac12\mathcal L_Xg = \lambda g is derived here (11B.1 Ricci Solitons).

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

  • find integral curves and flows of vector fields, and state the fundamental theorem on flows;
  • compute Lie brackets in coordinates and explain the bracket as the failure of flows to commute;
  • define the Lie derivative of functions, vector fields and tensors, and compute LXg\mathcal L_Xg;
  • recognise Killing fields and conformal fields;
  • state the naturality of the Ricci tensor and derive the soliton equation from it.

How a fluid strains

In the world Model The rate-of-strain tensor

A fluid moving with velocity field uu carries every small blob of fluid along: the blob is translated, rotated and deformed. The translation and rotation don't change distances between fluid particles; the deformation does. Its rate is measured by the rate-of-strain tensor, in Cartesian coordinates

Dij=12(∂ui∂xj+∂uj∂xi),D_{ij} = \frac12\Big(\frac{\partial u_i}{\partial x^j} + \frac{\partial u_j}{\partial x^i}\Big),

the symmetric part of the velocity gradient. A small disc of fluid becomes, after a short time, an ellipse whose axes point along the eigenvectors of DD, stretched or squashed at rates given by its eigenvalues (Figure 6.1). In a Newtonian fluid such as water or air, the viscous stress is proportional to it: the Navier–Stokes equations use the stress σ=−pI+2μD\sigma = -pI + 2\mu D, with μ\mu the viscosity.

This tensor is exactly 12Lug\frac12\mathcal L_ug, half the Lie derivative of the Euclidean metric along the flow of uu (Exercise 6.7): the rate at which the flow changes the metric. It is an identity, not an analogy. A flow with D=0D = 0 moves the fluid rigidly: its velocity field is a Killing field. And the same object, LXg\mathcal L_Xg for a vector field XX on a Riemannian manifold, is the extra term in the Ricci soliton equation.

Figure 6.1. A disc of fluid under the linear flow u=(x+0.6y,0.6x−0.3y)u = (x + 0.6y, 0.6x - 0.3y) (computed). Its rate-of-strain tensor D=12LugD = \frac12\mathcal L_ug has eigenvectors along which the disc stretches (orange) and squashes (blue); after time 0.40.4 the disc has become the ellipse shown. Thin lines: streamlines of the flow.

Integral curves and flows

A vector field XX on MM assigns to each pp a tangent vector Xp∈TpMX_p \in T_pM, smoothly (8A.3 Tangent Vectors and Bundles); in coordinates X=Xi∂iX = X^i\partial_i with smooth functions XiX^i. An integral curve is a curve γ\gamma with γ′(t)=Xγ(t)\gamma'(t) = X_{\gamma(t)}: in coordinates, the ODE system γ˙i=Xi(γ)\dot\gamma^i = X^i(\gamma).

Theorem 6.1 The fundamental theorem on flows

For a smooth vector field XX on MM there is a unique maximal flow: an open set D⊆R×M\mathcal D \subseteq \mathbb{R}\times M containing {0}×M\{0\}\times M and a smooth map θ:D→M\theta : \mathcal D \to M, (t,p)↦θt(p)(t, p) \mapsto \theta_t(p), such that t↦θt(p)t \mapsto \theta_t(p) is the maximal integral curve starting at pp, and θt∘θs=θt+s\theta_t\circ\theta_s = \theta_{t+s} where defined. If XX has compact support (in particular, if MM is compact), the flow is defined for all tt and each θt\theta_t is a diffeomorphism.

The existence, uniqueness and smooth dependence on the starting point come from the ODE theory of 2B.10 Ordinary Differential Equations, applied in charts; maximality patches the local solutions together. If XX has compact support, every integral curve stays in a compact set, so by the blow-up alternative (2B.10 Ordinary Differential Equations) it exists for all time. A field whose flow exists for all time is complete. On R\mathbb{R}, the field x2∂xx^2\partial_x is not: its integral curve x(t)=x01−x0tx(t) = \frac{x_0}{1 - x_0t} escapes to infinity in finite time.

Time-dependent vector fields XtX_t have flows too, solving ddtϕt(p)=Xt(ϕt(p))\frac{d}{dt}\phi_t(p) = X_t(\phi_t(p)), though they no longer form a group. They are the diffeomorphisms of DeTurck's trick (11A.3 Short-Time Existence and Uniqueness) and of the soliton equation below.

The Lie bracket

The bracket of two vector fields is the vector field

[X,Y]f=X(Yf)−Y(Xf),[X,Y]=(Xi∂iYj−Yi∂iXj)∂j[X, Y]f = X(Yf) - Y(Xf), \qquad [X, Y] = \big(X^i\partial_iY^j - Y^i\partial_iX^j\big)\partial_j

(Exercise 6.4). Second derivatives of ff cancel, so [X,Y][X, Y] is again a first-order operator: a vector field. The bracket is antisymmetric and satisfies the Jacobi identity, and for left-invariant fields on a Lie group it is the bracket of the Lie algebra (8A.5 Lie Groups and Group Actions). Coordinate vector fields commute: [∂i,∂j]=0[\partial_i, \partial_j] = 0.

The geometric meaning. Flow along XX for time tt, then along YY for time tt, then back along XX, then back along YY. To second order you arrive at

ψ−tθ−tψtθt(p)=p+t2[X,Y]p+O(t3)\psi_{-t}\theta_{-t}\psi_t\theta_t(p) = p + t^2[X, Y]_p + O(t^3)

in coordinates (Lee, chapter 9): the bracket measures the failure of the flows to commute. In particular two flows commute, θtψs=ψsθt\theta_t\psi_s = \psi_s\theta_t, if and only if [X,Y]=0[X, Y] = 0.

In the world Model Parallel parking

A car's position and heading form a three-dimensional configuration (x,y,ϕ)(x, y, \phi), but it has only two controls: drive (move forward along the heading) and steer (change the heading, for this idealised car, while stationary). The corresponding vector fields are D=cos⁡ϕ ∂x+sin⁡ϕ ∂yD = \cos\phi\,\partial_x + \sin\phi\,\partial_y and S=∂ϕS = \partial_\phi. Neither moves the car sideways. Their bracket does:

[S,D]=−sin⁡ϕ ∂x+cos⁡ϕ ∂y,[S, D] = -\sin\phi\,\partial_x + \cos\phi\,\partial_y,

the sideways direction (Exercise 6.5). So the sequence "turn, drive forward, turn back, drive back", repeated in small steps, moves the car sideways: parallel parking (Figure 6.2). Control systems whose reachable motions are generated by brackets of the available ones are called nonholonomic; the theory, built on this fact (Chow–Rashevskii theorem), is standard in robotics (Murray, Li and Sastry, A Mathematical Introduction to Robotic Manipulation, 1994, chapter 7).

Figure 6.2. A parallel-parking cycle for the idealised car (computed): rotate by ε\varepsilon, drive forward by ε\varepsilon, rotate back, drive back. With ε=0.5\varepsilon = 0.5 the car returns to its original heading displaced by (−0.061,0.240)(-0.061, 0.240): sideways by about ε2=0.25\varepsilon^2 = 0.25, the direction of [S,D][S, D], plus a smaller correction of order ε3\varepsilon^3.

The Lie derivative

Let θt\theta_t be the flow of XX. The pullback of a covariant tensor field TT (a function, a covector field, a metric) by a diffeomorphism ϕ\phi is (ϕ∗T)p(v1,…,vk)=Tϕ(p)(dϕ(v1),…,dϕ(vk))(\phi^*T)_p(v_1, \dots, v_k) = T_{\phi(p)}(d\phi(v_1), \dots, d\phi(v_k)); a vector field is pulled back by (ϕ∗Y)p=(dϕp)−1Yϕ(p)(\phi^*Y)_p = (d\phi_p)^{-1}Y_{\phi(p)}.

Definition 6.2 Lie derivative

The Lie derivative of a tensor field TT along XX is

LXT=ddt∣t=0θt∗T.\mathcal L_XT = \frac{d}{dt}\Big|_{t=0}\theta_t^*T.

It is a tensor of the same type. The cases that recur:

  • functions: LXf=Xf=df(X)\mathcal L_Xf = Xf = df(X);
  • vector fields: LXY=[X,Y]\mathcal L_XY = [X, Y];
  • covariant 2-tensors, in particular a metric: in coordinates,
(LXg)ij=Xk∂kgij+gkj∂iXk+gik∂jXk,(\mathcal L_Xg)_{ij} = X^k\partial_kg_{ij} + g_{kj}\partial_iX^k + g_{ik}\partial_jX^k,

and in Euclidean space with Cartesian coordinates (LXg)ij=∂iXj+∂jXi(\mathcal L_Xg)_{ij} = \partial_iX_j + \partial_jX_i (Exercise 6.6); on a Riemannian manifold, in terms of the covariant derivative of 9A.2 Connections, (LXg)ij=∇iXj+∇jXi(\mathcal L_Xg)_{ij} = \nabla_iX_j + \nabla_jX_i;

The Lie derivative along the flow gives the derivative at every time, not just t=0t = 0:

ddtθt∗T=θt∗(LXT),\frac{d}{dt}\theta_t^*T = \theta_t^*(\mathcal L_XT),

by the group property θt+s∗=θt∗θs∗\theta_{t+s}^* = \theta_t^*\theta_s^* (Exercise 6.9).

Killing and conformal fields. A vector field with LXg=0\mathcal L_Xg = 0 is a Killing field: its flow consists of isometries. On R2\mathbb{R}^2 the rotation field −y∂x+x∂y-y\partial_x + x\partial_y is Killing; the dilation x∂x+y∂yx\partial_x + y\partial_y is not, but it satisfies LXg=2g\mathcal L_Xg = 2g: a conformal field, whose flow scales all lengths by the same factor (Exercise 6.8).

Naturality and its consequences

Theorem 6.3 Naturality of the Ricci tensor

For every diffeomorphism ϕ\phi and every Riemannian metric gg,

Ric⁡(ϕ∗g)=ϕ∗Ric⁡(g).\operatorname{Ric}(\phi^*g) = \phi^*\operatorname{Ric}(g).

The proof is in 9A.4 Curvature and What It Means; the reason is that the Ricci tensor is built from gg by a coordinate-free construction (the Levi-Civita connection and its curvature), so it cannot tell a metric from its reparametrisation. Every geometric quantity has this property: scalar curvature, volume, the Laplacian. Two consequences follow, and between them they shape the analysis of the Ricci flow (Figure 6.3).

1. The flow is invariant under diffeomorphisms. If g(t)g(t) solves ∂tg=−2Ric⁡(g)\partial_tg = -2\operatorname{Ric}(g), so does ϕ∗g(t)\phi^*g(t) for any fixed ϕ\phi: ∂tϕ∗g=ϕ∗∂tg=−2ϕ∗Ric⁡(g)=−2Ric⁡(ϕ∗g)\partial_t\phi^*g = \phi^*\partial_tg = -2\phi^*\operatorname{Ric}(g) = -2\operatorname{Ric}(\phi^*g). So solutions come in huge families, one for each diffeomorphism, and the linearised equation has a kernel in the directions LXg\mathcal L_Xg. An equation with this degeneracy cannot be strictly parabolic (6A.1 What a PDE Is, 6A.7 Nonlinear Parabolic Equations). DeTurck's trick breaks the symmetry by adding a term LW(g)g\mathcal L_{W(g)}g, solves the resulting strictly parabolic equation, and recovers a Ricci flow by pulling back along the time-dependent flow of WW (11A.3 Short-Time Existence and Uniqueness).

2. Solitons. The most symmetric solutions change only by diffeomorphisms and scaling:

g(t)=σ(t) ϕt∗g0,g(t) = \sigma(t)\,\phi_t^*g_0,

with σ(t)>0\sigma(t) > 0, σ(0)=1\sigma(0) = 1, and ϕt\phi_t the flow of a (time-dependent) vector field with ϕ0=id⁡\phi_0 = \operatorname{id}. Differentiating at t=0t = 0, using ddtϕt∗g0∣t=0=LXg0\frac{d}{dt}\phi_t^*g_0|_{t=0} = \mathcal L_Xg_0:

−2Ric⁡(g0)=σ′(0)g0+LXg0,that isRic⁡(g0)+12LXg0=λg0,-2\operatorname{Ric}(g_0) = \sigma'(0)g_0 + \mathcal L_Xg_0, \qquad\text{that is}\qquad \operatorname{Ric}(g_0) + \frac12\mathcal L_Xg_0 = \lambda g_0,

with λ=−12σ′(0)\lambda = -\frac12\sigma'(0). This is the Ricci soliton equation. The soliton is shrinking, steady or expanding as λ>0\lambda > 0, =0= 0 or <0< 0. When X=∇fX = \nabla f is a gradient, L∇fg=2∇2f\mathcal L_{\nabla f}g = 2\nabla^2f and the equation becomes Ric⁡+∇2f=λg\operatorname{Ric} + \nabla^2f = \lambda g (Exercise 6.10). Solitons are the fixed points of the Ricci flow modulo symmetry, and they are the models for its singularities (11B.1 Ricci Solitons, 12A.3 The 𝓦-Entropy).

Figure 6.3. Naturality of the Ricci tensor and its two consequences: diffeomorphism invariance, which makes the flow only weakly parabolic and is repaired by DeTurck's trick; and self-similar solutions, which satisfy the Ricci soliton equation.

History

Sophus Lie introduced the bracket and the derivative now named after him in his work on transformation groups in the 1870s–80s; the term "Lie derivative" and the modern formulation are due to Władysław Ślebodziński (1931) and David van Dantzig. Wilhelm Killing introduced the fields named after him in 1888. Wei-Liang Chow (1939) and Petr Rashevskii (1938) proved that brackets of the available directions determine the reachable configurations of a nonholonomic system. Richard Hamilton introduced Ricci solitons in the 1980s; Dennis DeTurck's trick appeared in 1983.

Recall Where we stand

A vector field has integral curves and a flow θt\theta_t, defined for all time when the field has compact support. The bracket [X,Y]=(Xi∂iYj−Yi∂iXj)∂j[X, Y] = (X^i\partial_iY^j - Y^i\partial_iX^j)\partial_j measures how the flows fail to commute; brackets of drive and steer give parallel parking. The Lie derivative LXT=ddtθt∗T∣0\mathcal L_XT = \frac{d}{dt}\theta_t^*T|_0 is XfXf on functions, [X,Y][X, Y] on vector fields, and ∂iXj+∂jXi\partial_iX_j + \partial_jX_i on the Euclidean metric: twice the rate of strain. Killing fields have LXg=0\mathcal L_Xg = 0. The Ricci tensor is natural, Ric⁡(ϕ∗g)=ϕ∗Ric⁡(g)\operatorname{Ric}(\phi^*g) = \phi^*\operatorname{Ric}(g), so the Ricci flow is diffeomorphism invariant (hence weakly parabolic) and has self-similar solutions satisfying Ric⁡+12LXg=λg\operatorname{Ric} + \frac12\mathcal L_Xg = \lambda g. 8A.7 Tensors and Index Notation sets up the tensor notation in which all of this is computed.

Exercises

Exercise 6.4 The bracket in coordinates

Show that [X,Y]f=(Xi∂iYj−Yi∂iXj)∂jf[X, Y]f = (X^i\partial_iY^j - Y^i\partial_iX^j)\partial_jf, and that the second derivatives cancel. Compute [x∂y,y∂x][x\partial_y, y\partial_x] on R2\mathbb{R}^2.

Solution

X(Yf)=Xi∂i(Yj∂jf)=Xi∂iYj∂jf+XiYj∂i∂jfX(Yf) = X^i\partial_i(Y^j\partial_jf) = X^i\partial_iY^j\partial_jf + X^iY^j\partial_i\partial_jf; subtract the same with XX and YY exchanged, and the symmetric second-derivative terms cancel. [x∂y,y∂x]=x∂y(y)∂x−y∂x(x)∂y=x∂x−y∂y[x\partial_y, y\partial_x] = x\partial_y(y)\partial_x - y\partial_x(x)\partial_y = x\partial_x - y\partial_y.

Exercise 6.5 Drive and steer

For D=cos⁡ϕ ∂x+sin⁡ϕ ∂yD = \cos\phi\,\partial_x + \sin\phi\,\partial_y and S=∂ϕS = \partial_\phi on R2×S1\mathbb{R}^2\times S^1, compute [S,D][S, D], and show that DD, SS, [S,D][S, D] span the tangent space at every point. Then compute the result of the four-step manoeuvre "steer by ε\varepsilon, drive ε\varepsilon, steer by −ε-\varepsilon, drive −ε-\varepsilon" from (0,0,0)(0, 0, 0) exactly, and compare with ε2[S,D]\varepsilon^2[S, D].

Solution

[S,D]=∂ϕ(cos⁡ϕ)∂x+∂ϕ(sin⁡ϕ)∂y=−sin⁡ϕ ∂x+cos⁡ϕ ∂y[S, D] = \partial_\phi(\cos\phi)\partial_x + \partial_\phi(\sin\phi)\partial_y = -\sin\phi\,\partial_x + \cos\phi\,\partial_y. The three are orthonormal in the (x,y,ϕ)(x, y, \phi) coordinates. The manoeuvre: after steering, heading ε\varepsilon; driving moves to (εcos⁡ε,εsin⁡ε,ε)(\varepsilon\cos\varepsilon, \varepsilon\sin\varepsilon, \varepsilon); steering back gives heading 00; driving back moves to (εcos⁡ε−ε,εsin⁡ε,0)=(−ε32+…,ε2+…,0)(\varepsilon\cos\varepsilon - \varepsilon, \varepsilon\sin\varepsilon, 0) = (-\frac{\varepsilon^3}{2} + \dots, \varepsilon^2 + \dots, 0). To leading order the displacement is ε2∂y=ε2[S,D]\varepsilon^2\partial_y = \varepsilon^2[S, D] at ϕ=0\phi = 0.

Exercise 6.6 The Lie derivative of a metric

Derive (LXg)ij=Xk∂kgij+gkj∂iXk+gik∂jXk(\mathcal L_Xg)_{ij} = X^k\partial_kg_{ij} + g_{kj}\partial_iX^k + g_{ik}\partial_jX^k by differentiating (θt∗g)ij(p)=gkl(θt(p)) ∂iθtk ∂jθtl(\theta_t^*g)_{ij}(p) = g_{kl}(\theta_t(p))\,\partial_i\theta_t^k\,\partial_j\theta_t^l at t=0t = 0, using θt(x)=x+tX(x)+O(t2)\theta_t(x) = x + tX(x) + O(t^2). Specialise to gij=δijg_{ij} = \delta_{ij}.

Exercise 6.7 Strain is a Lie derivative

For a velocity field uu on R3\mathbb{R}^3, show that 12(Lug)ij=Dij\frac12(\mathcal L_ug)_{ij} = D_{ij}, the rate-of-strain tensor. Show that a rigid motion u(x)=ω×x+bu(x) = \omega\times x + b has D=0D = 0: rotations and translations are Killing fields of Euclidean space.

Exercise 6.8 Dilations and rotations

On R2\mathbb{R}^2, show that X=x∂x+y∂yX = x\partial_x + y\partial_y has LXg=2g\mathcal L_Xg = 2g and that R=−y∂x+x∂yR = -y\partial_x + x\partial_y has LRg=0\mathcal L_Rg = 0. What are their flows, and how do they act on the length of a curve?

Exercise 6.9 The derivative at every time

Using θt+s=θt∘θs\theta_{t+s} = \theta_t\circ\theta_s, show that ddtθt∗T=θt∗(LXT)\frac{d}{dt}\theta_t^*T = \theta_t^*(\mathcal L_XT).

Solution

ddtθt∗T=dds∣s=0θt+s∗T=dds∣s=0θt∗θs∗T=θt∗dds∣s=0θs∗T\frac{d}{dt}\theta_t^*T = \frac{d}{ds}\Big|_{s=0}\theta_{t+s}^*T = \frac{d}{ds}\Big|_{s=0}\theta_t^*\theta_s^*T = \theta_t^*\frac{d}{ds}\Big|_{s=0}\theta_s^*T, since θt∗\theta_t^* is linear and doesn't depend on ss.

Exercise 6.10 Rehearsal: flat space is a shrinking soliton

On Rn\mathbb{R}^n with the Euclidean metric gg, let τ>0\tau > 0 and f=∣x∣24τf = \frac{|x|^2}{4\tau}. (a) Show ∇f=x2τ\nabla f = \frac{x}{2\tau} and L∇fg=2∇2f=1τg\mathcal L_{\nabla f}g = 2\nabla^2f = \frac1\tau g. (b) Since Ric⁡(g)=0\operatorname{Ric}(g) = 0, deduce Ric⁡+∇2f=12τg\operatorname{Ric} + \nabla^2f = \frac{1}{2\tau}g: flat space with this ff is a shrinking soliton, the Gaussian soliton, with λ=12τ\lambda = \frac{1}{2\tau}. (c) Reconcile with the fact that flat space doesn't change under Ricci flow: write the solution as g(t)=σ(t)ϕt∗gg(t) = \sigma(t)\phi_t^*g with σ(t)=1−tτ\sigma(t) = 1 - \frac t\tau and ϕt\phi_t a scaling of Rn\mathbb{R}^n that exactly undoes it. This ff is the one in Perelman's Gaussian (6A.3 The Heat Equation on ℝⁿ, 12A.3 The 𝓦-Entropy).

Solution

(a) ∂if=xi2τ\partial_if = \frac{x_i}{2\tau}, ∂i∂jf=δij2τ\partial_i\partial_jf = \frac{\delta_{ij}}{2\tau}, and in Euclidean coordinates (LYg)ij=∂iYj+∂jYi(\mathcal L_Yg)_{ij} = \partial_iY_j + \partial_jY_i, so L∇fg=2∇2f=1τg\mathcal L_{\nabla f}g = 2\nabla^2f = \frac1\tau g. (b) Immediate. (c) ϕt(x)=x1−t/τ\phi_t(x) = \frac{x}{\sqrt{1 - t/\tau}} multiplies the Euclidean metric by 11−t/τ\frac{1}{1 - t/\tau}, so σ(t)ϕt∗g=g\sigma(t)\phi_t^*g = g for all tt: the flat metric is constant, as it must be, but it can also be read as shrinking by σ\sigma while being stretched back by ϕt\phi_t. Differentiating at 00: σ′(0)=−1τ\sigma'(0) = -\frac1\tau, so λ=12τ\lambda = \frac{1}{2\tau}.

© 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.