Lie Derivatives and Killing Vectors

Ask how a tensor field changes as you slide along a flow, and you hit a subtlety ordinary calculus hides: the field at a nearby point lives in a different vector space than the field where you stand, so you cannot subtract the two. The covariant derivative fixes this by importing a connection — an extra rule for carrying tensors between tangent spaces. The Lie derivative fixes it with nothing beyond the flow itself. A vector field XX generates a family of diffeomorphisms ϕt\phi_t that drag every geometric object along; to differentiate TT, you let that drag carry the value from the future point back to where you are, compare, and take the rate at t=0t=0. No connection, no metric, no coordinates enter — which is why LX\mathcal{L}_X is the natural derivative of bare differential topology, while ∇\nabla belongs to Riemannian geometry. Two special cases fix everything. On a function the Lie derivative is the directional derivative Xμ∂μfX^\mu\partial_\mu f; on a vector field it is the commutator [X,Y][X,Y], which measures how badly the flows of XX and YY fail to commute. Demand that LX\mathcal{L}_X obey the Leibniz rule and commute with contraction, and its action on any tensor is forced. The payoff is physical. A Killing vector is a direction along which the drag leaves the metric fixed, Lξg=0\mathcal{L}_\xi g=0 — its flow is an isometry, a continuous symmetry of the geometry. Killing's equation ∇μξν+∇νξμ=0\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=0 turns "the metric has a symmetry" into a solvable first-order system, and every solution buys a conserved quantity ξμuμ\xi_\mu u^\mu along geodesics — the constants that make orbits, light-bending, and black-hole shadows computable.

Key takeaways
  • The Lie derivative LXT\mathcal{L}_X T measures how a tensor field changes as the flow of XX drags it, comparing TT at a point with its push-forward from a neighbouring point — a metric-free, connection-free operation.
  • On a scalar it is the directional derivative Xμ∂μfX^\mu\partial_\mu f; on a vector it is the Lie bracket [X,Y][X,Y]; the Leibniz rule then fixes it on every tensor.
  • The bracket [X,Y][X,Y] is the failure of flows to close: run along XX then YY then back and you miss the start by ϵ2[X,Y]\epsilon^2[X,Y].
  • A Killing vector satisfies Lξgμν=∇μξν+∇νξμ=0\mathcal{L}_\xi g_{\mu\nu}=\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=0 — its flow is an isometry, encoding a continuous symmetry of the metric.
  • Each Killing vector gives a conserved charge ξμuμ\xi_\mu u^\mu along any geodesic; Schwarzschild's ∂t\partial_t and ∂ϕ\partial_\phi yield conserved energy and angular momentum.
  • A maximally symmetric nn-manifold admits exactly 12n(n+1)\tfrac12 n(n+1) Killing vectors — 1010 for 4D Minkowski (the Poincaré algebra), 33 for the sphere S2S^2.
  • Key idea: The Lie derivative differentiates a tensor field by dragging it with the flow of XX — a metric-free, connection-free operation; a Killing vector is one whose drag leaves the metric fixed, so Lξgμν=∇μξν+∇νξμ=0\mathcal{L}_\xi g_{\mu\nu}=\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=0, encoding a continuous isometry and hence a conserved quantity.

Where this lives today

Killing vectors are the daily currency of relativistic astrophysics. Every stationary, axisymmetric black-hole spacetime carries ∂t\partial_t and ∂ϕ\partial_\phi, and the conserved EE and LL they generate reduce the geodesic equations to a one-dimensional radial problem — the machinery behind the Event Horizon Telescope's images of M87 and Sgr A, whose shadow radius is the 33 GM/c23\sqrt3\,GM/c^2 impact parameter computed above. In cosmology, the observed homogeneity and isotropy of the universe is the statement that space is maximally symmetric — six Killing vectors per spatial slice — forcing the FLRW metric of the concordance model, while de Sitter and anti-de Sitter space, each with the full 1010 isometries, anchor inflation and the AdS/CFT correspondence. Closer to Earth, Cartan's formula drives the transport identities of continuum mechanics — the vorticity two-form's vanishing Lie derivative is Helmholtz's frozen-in theorem — and the "3+1" codes of numerical relativity.

The physical picture

Think of XX as the velocity field of a fluid whose flow ϕt\phi_t carries every particle — and every field painted on it — downstream. To Lie-differentiate a tensor you do not compare its values at two points directly; you let the flow carry the value from ϕt(p)\phi_t(p) back to pp, lay it beside the value already there, and read off the rate of mismatch. A quantity advected without distortion has zero Lie derivative; a nonzero LXT\mathcal{L}_X T is the rate at which the flow stretches, shears, or rotates TT. The bracket [X,Y][X,Y] is the sharpest instance: what is left over when two flows refuse to close a little coordinate box.

Figure
Flow along X = ∂ₓ then Y = x∂_y then back along each does not return to the start. The order-ε² gap drawn in red is the Lie bracket of X and Y, here equal to ∂_y, so the loop misses by ε²∂_y — vector fields commute exactly when this gap vanishes
The Lie bracket is the failure of two flows to close a loop p flow ε along X = ∂ₓ flow ε along Y back −ε along X back −ε along Y gap = ε²∂_y Y = x∂_y bracket of X, Y = ∂_y ≠ 0 loop closes only when the bracket vanishes

A Killing vector is the opposite of a distorting flow: a direction in which the geometry looks exactly the same before and after. Slide the entire metric along ξ\xi and it superimposes on itself perfectly, every length and angle preserved. On a round sphere the rotations do this; on a generic bumpy surface no rigid motion exists and there are no Killing vectors at all. Killing's equation ∇(μξν)=0\nabla_{(\mu}\xi_{\nu)}=0 says the flow has no shear and no expansion, only rotation — the symmetric part of ∇ξ\nabla\xi, which encodes stretching, is forced to vanish. The count of independent Killing vectors is a hard invariant of the space, maxing out at 12n(n+1)\tfrac12 n(n+1).

Figure
A round 2-sphere looks identical after any rotation, so each rotation axis furnishes a Killing vector: n(n+1)/2 = 3 in all, the generators of SO(3). Dragging the metric along ξ = ∂_φ leaves ds² unchanged, so its Lie derivative vanishes; a bumpy surface has none
A Killing vector: a direction the metric does not notice axis of ξ = ∂_φ isometry flow Killing count of S² ½·n(n+1) = 3 rotations about x, y, z the generators of SO(3) drag g along ξ ⇒ ds² fixed

The core equations

The definition is a limit built from the flow. The tensor at ϕt(p)\phi_t(p) cannot be subtracted from the tensor at pp directly, so the flow's own push-forward (ϕ−t)∗(\phi_{-t})_* carries it back to pp first; the rate of the resulting comparison at t=0t=0 is the Lie derivative — using nothing but the manifold's differentiable structure.

LXT  =  lim⁡t→0(ϕ−t)∗Tϕt(p)−Tpt  =  ddt∣t=0(ϕ−t)∗T(1)\mathcal{L}_X T \;=\; \lim_{t\to 0} \frac{(\phi_{-t})_* T_{\phi_t(p)} - T_p}{t} \;=\; \left.\frac{\mathrm{d}}{\mathrm{d}t}\right|_{t=0} (\phi_{-t})_* T \qquad (1)

where: XX — the vector field generating the flow; ϕt\phi_t — its one-parameter group of diffeomorphisms; (ϕ−t)∗(\phi_{-t})_* — the push-forward that drags a tensor from ϕt(p)\phi_t(p) back to pp; TT — an arbitrary (p,q)(p,q) tensor field. Read it as "compare TT with its dragged-back self and take the rate" — the drag-back is what makes the subtraction legal. On a scalar there is nothing to drag but the value itself, so the Lie derivative collapses to the plain directional derivative — the rate of change of ff along XX.

LXf=Xμ∂μf=X(f)(2)\mathcal{L}_X f = X^\mu \partial_\mu f = X(f) \qquad (2)

where: ff — a smooth function; XμX^\mu — the components of XX in a coordinate basis. This is the one case where LX\mathcal{L}_X and ∇X\nabla_X agree identically — a scalar has no indices for a connection to act on. On a vector field the flow drags YY as well as evaluates it, and the extra term −Yν∂νXμ-Y^\nu\partial_\nu X^\mu is exactly the commutator of the two fields — the Lie bracket.

(LXY)μ=Xν∂νYμ−Yν∂νXμ=[X,Y]μ(3)(\mathcal{L}_X Y)^\mu = X^\nu \partial_\nu Y^\mu - Y^\nu \partial_\nu X^\mu = [X,Y]^\mu \qquad (3)

where: YμY^\mu — the components of the field being differentiated; [X,Y][X,Y] — the Lie bracket, the commutator of XX and YY as first-order operators. It is antisymmetric and Jacobi-satisfying, so vector fields form a Lie algebra under it. A one-form is dual to a vector, and the Leibniz rule fixes its Lie derivative once the scalar and vector cases are known — the XX-derivative term switching sign because a one-form carries a lower index.

(LXω)μ=Xν∂νωμ+ων∂μXν(4)(\mathcal{L}_X \omega)_\mu = X^\nu \partial_\nu \omega_\mu + \omega_\nu \partial_\mu X^\nu \qquad (4)

where: ωμ\omega_\mu — the components of the one-form (a lower, covariant index); the second term derives XX rather than ω\omega. The rule generalizes: every upper index contributes a −(∂X)-(\partial X) term and every lower index a +(∂X)+(\partial X) term. Assembling these, the Lie derivative of a general (p,q)(p,q) tensor is a transport term plus one correction per index. This single formula reproduces (2), (3), and (4) as special cases.

LXTμ1⋯ν1⋯=Xλ∂λTμ1⋯ν1⋯−∑iT⋯λ⋯ν⋯ ∂λXμi+∑jTμ⋯⋯λ⋯ ∂νjXλ(5)\mathcal{L}_X T^{\mu_1\cdots}{}_{\nu_1\cdots} = X^\lambda \partial_\lambda T^{\mu_1\cdots}{}_{\nu_1\cdots} - \sum_i T^{\cdots\lambda\cdots}{}_{\nu\cdots}\,\partial_\lambda X^{\mu_i} + \sum_j T^{\mu\cdots}{}_{\cdots\lambda\cdots}\,\partial_{\nu_j} X^{\lambda} \qquad (5)

where: Tμ⋯ν⋯T^{\mu\cdots}{}_{\nu\cdots} — the tensor field; ∂λXμ\partial_\lambda X^\mu — the derivatives of the generator that source every correction term. Because only ∂X\partial X appears, the Lie derivative depends on XX throughout a neighbourhood, not just at a point. Remarkably, the ordinary derivatives in (5) may be swapped for covariant ones with any torsion-free connection and the answer is unchanged, because the Christoffel contributions cancel by the symmetry Γμνλ=Γνμλ\Gamma^\lambda_{\mu\nu}=\Gamma^\lambda_{\nu\mu}. This is the algebraic proof that LX\mathcal{L}_X needs no connection.

LXT=Xλ∇λT−(∇λXμ) T⋯λ⋯+(∇νXλ) T⋯⋯λ⋯(6)\mathcal{L}_X T = X^\lambda \nabla_\lambda T - (\nabla_\lambda X^{\mu})\,T^{\cdots\lambda\cdots} + (\nabla_{\nu} X^{\lambda})\,T^{\cdots}{}_{\cdots\lambda\cdots} \qquad (6)

where: ∇\nabla — any symmetric (torsion-free) connection; the result is independent of which one, and applied to the metric it yields Killing's equation in a line. A Killing vector is defined by demanding that dragging the metric change nothing: Lξg=0\mathcal{L}_\xi g=0. Feeding gg through (6) and using metric compatibility ∇g=0\nabla g=0 turns this into the vanishing of the symmetrized covariant derivative of ξ\xi.

Lξgμν=∇μξν+∇νξμ=2 ∇(μξν)=0(7)\mathcal{L}_\xi g_{\mu\nu} = \nabla_\mu \xi_\nu + \nabla_\nu \xi_\mu = 2\,\nabla_{(\mu}\xi_{\nu)} = 0 \qquad (7)

where: ξμ=gμνξν\xi_\mu=g_{\mu\nu}\xi^\nu — the covariant components of the Killing vector; gμνg_{\mu\nu} — the metric; ∇(μξν)\nabla_{(\mu}\xi_{\nu)} — the symmetrized derivative. It is an overdetermined first-order system whose solutions form the finite-dimensional Lie algebra of infinitesimal isometries. For differential forms there is a slicker route: Cartan's "magic" formula expresses LX\mathcal{L}_X on any pp-form through the exterior derivative and the interior product, with no coordinates at all.

LXα=d (ιXα)+ιX (dα)(8)\mathcal{L}_X \alpha = \mathrm{d}\,(\iota_X \alpha) + \iota_X\,(\mathrm{d}\alpha) \qquad (8)

where: α\alpha — a differential pp-form; d\mathrm{d} — the exterior derivative; ιX\iota_X — the interior product (contraction of XX into the first slot). It packages "differentiate then contract" and its reverse into the flow derivative. The physical harvest is a conservation law: contract a Killing vector with the tangent to a geodesic and the result is constant along the curve — every symmetry becomes a first integral of the motion.

Q=ξμ dxμdτ=ξμuμ=const along a geodesic(9)Q = \xi_\mu\, \frac{\mathrm{d}x^\mu}{\mathrm{d}\tau} = \xi_\mu u^\mu = \text{const along a geodesic} \qquad (9)

where: uμ=dxμ/dτu^\mu=\mathrm{d}x^\mu/\mathrm{d}\tau — the geodesic tangent (four-velocity); τ\tau — an affine parameter (proper time for timelike geodesics); QQ — the conserved charge tied to ξ\xi.

Derivation AThe Lie bracket as the failure of two flows to close

The goal here is to give the bracket [X,Y][X,Y] its geometric meaning as the obstruction to commuting flows, and thereby to explain why commuting vector fields — and only they — can serve as a coordinate grid. We track a single point through a small loop: flow along XX for parameter ϵ\epsilon, then YY, then back along XX, then back along YY, and show it lands a distance ϵ2[X,Y]\epsilon^2[X,Y] from its start, so the bracket is exactly the order-ϵ2\epsilon^2 non-closure of the flow rectangle.

Step 1. Expand each flow to second order from its defining equation x˙μ=Xμ(x)\dot x^\mu = X^\mu(x), giving ϕϵX(x)μ=xμ+ϵXμ(x)+12ϵ2Xν∂νXμ+O(ϵ3)\phi^X_\epsilon(x)^\mu = x^\mu + \epsilon X^\mu(x) + \tfrac12\epsilon^2 X^\nu\partial_\nu X^\mu + O(\epsilon^3), and likewise for YY. Why: the second-order term, where each field is re-evaluated at the already-displaced point, is what carries the bracket.

Step 2. Compose the four maps in order ϕ−ϵY∘ϕ−ϵX∘ϕϵY∘ϕϵX\phi^Y_{-\epsilon}\circ\phi^X_{-\epsilon}\circ\phi^Y_{\epsilon}\circ\phi^X_{\epsilon}, at each stage evaluating the next field at the point reached so far and keeping terms through O(ϵ2)O(\epsilon^2). Why: evaluating YY at the XX-displaced point rather than at pp is what breaks the symmetry between the two orders.

Step 3. Collect the surviving terms. Every first-order displacement cancels between the forward and backward legs, and the 12ϵ2Xν∂νXμ\tfrac12\epsilon^2 X^\nu\partial_\nu X^\mu self-terms cancel too, leaving ϕ−ϵY∘ϕ−ϵX∘ϕϵY∘ϕϵX(p)=p+ϵ2(Xν∂νYμ−Yν∂νXμ)+O(ϵ3)\phi^Y_{-\epsilon}\circ\phi^X_{-\epsilon}\circ\phi^Y_{\epsilon}\circ\phi^X_{\epsilon}(p) = p + \epsilon^2\big(X^\nu\partial_\nu Y^\mu - Y^\nu\partial_\nu X^\mu\big) + O(\epsilon^3). Why: the mixed term Xν∂νYμ−Yν∂νXμX^\nu\partial_\nu Y^\mu - Y^\nu\partial_\nu X^\mu is exactly the bracket components [X,Y]μ[X,Y]^\mu, so the gap is ϵ2[X,Y]p\epsilon^2[X,Y]_p.

Step 4. Read off the meaning: the loop closes to order ϵ2\epsilon^2 if and only if [X,Y]=0[X,Y]=0, in which case the two flows commute and X,YX,Y can be completed to a commuting coordinate basis ∂a,∂b\partial_a,\partial_b. Why: this is the geometric content of LXY=[X,Y]\mathcal{L}_X Y=[X,Y] — the Lie derivative of YY along XX is the rate at which the XX-flow drags YY out of alignment, which is the same non-closure seen infinitesimally. *As a check, the explicit example X=∂x, Y=x∂yX=\partial_x,\ Y=x\partial_y gives [X,Y]=∂y[X,Y]=\partial_y, the gap drawn in the first figure.*

Go deeper · Derivation B · 1 minFrom metric invariance to Killing's equation
Derivation BFrom metric invariance to Killing's equation

The aim is to convert the geometric statement "the flow of ξ\xi is an isometry", written compactly as Lξg=0\mathcal{L}_\xi g=0, into the working equation ∇μξν+∇νξμ=0\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=0 that one solves to find symmetries. The bridge is the covariant form of the Lie derivative plus the metric compatibility of the Levi-Civita connection, which collapse three terms into a clean symmetrized derivative and expose Killing's equation as the vanishing of the shear-and-expansion part of ∇ξ\nabla\xi.

Step 1. Apply the covariant Lie-derivative formula (6) to the (0,2)(0,2) tensor gg, giving Lξgμν=ξλ∇λgμν+gλν∇μξλ+gμλ∇νξλ\mathcal{L}_\xi g_{\mu\nu} = \xi^\lambda\nabla_\lambda g_{\mu\nu} + g_{\lambda\nu}\nabla_\mu\xi^\lambda + g_{\mu\lambda}\nabla_\nu\xi^\lambda, one correction term for each lower index of gg. Why: using ∇\nabla rather than ∂\partial is legal because the Lie derivative is connection-independent, and it lets metric compatibility do the next step's work.

Step 2. Kill the transport term. For the Levi-Civita connection ∇λgμν=0\nabla_\lambda g_{\mu\nu}=0 identically, so the first term drops and only the two index-correction terms remain. Why: metric compatibility says the metric is covariantly constant, which is what makes the ξλ∇λg\xi^\lambda\nabla_\lambda g piece vanish.

Step 3. Lower the free index with gg inside the covariant derivative — again permitted by compatibility — to write gλν∇μξλ=∇μξνg_{\lambda\nu}\nabla_\mu\xi^\lambda=\nabla_\mu\xi_\nu and likewise for the other term, yielding Lξgμν=∇μξν+∇νξμ\mathcal{L}_\xi g_{\mu\nu}=\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu. Why: this is the symmetrized covariant derivative 2∇(μξν)2\nabla_{(\mu}\xi_{\nu)}, whose vanishing is Killing's equation.

Step 4. Demand isometry, Lξg=0\mathcal{L}_\xi g=0, to obtain ∇(μξν)=0\nabla_{(\mu}\xi_{\nu)}=0. Why: the symmetric part of ∇ξ\nabla\xi is the strain of the flow (its expansion and shear); setting it to zero says the flow moves the geometry rigidly, preserving every length and angle — the infinitesimal statement of an isometry. *Note that the antisymmetric part ∇[μξν]\nabla_{[\mu}\xi_{\nu]} is unconstrained, consistent with a Killing flow carrying rotation but no strain.*

Go deeper · Derivation C · 1 minWhy a Killing vector conserves a charge along geodesics
Derivation CWhy a Killing vector conserves a charge along geodesics

The target is the single most useful consequence of Killing symmetry: that the contraction of a Killing vector with the tangent to any geodesic is constant along the curve, so each symmetry supplies one first integral of the motion. The proof is a two-line collaboration between the geodesic equation, which kills the derivative of the tangent, and Killing's equation, which kills the symmetric derivative of ξ\xi — together forcing the total derivative of ξμuμ\xi_\mu u^\mu to vanish.

Step 1. Differentiate the scalar ξμuμ\xi_\mu u^\mu along the curve using the tangent as the directional operator: ddτ(ξμuμ)=uν∇ν(ξμuμ)\dfrac{\mathrm{d}}{\mathrm{d}\tau}(\xi_\mu u^\mu) = u^\nu\nabla_\nu(\xi_\mu u^\mu). Why: along a curve with tangent uu, the ordinary τ\tau-derivative of any scalar equals uν∇νu^\nu\nabla_\nu acting on it, which lets covariant machinery take over.

Step 2. Apply the Leibniz rule to split the derivative into a piece acting on ξ\xi and a piece acting on uu: uν∇ν(ξμuμ)=uνuμ∇νξμ+ξμ uν∇νuμu^\nu\nabla_\nu(\xi_\mu u^\mu) = u^\nu u^\mu\nabla_\nu\xi_\mu + \xi_\mu\,u^\nu\nabla_\nu u^\mu. Why: isolating the two factors lets each defining equation strike its own term.

Step 3. Use the geodesic equation uν∇νuμ=0u^\nu\nabla_\nu u^\mu=0 to eliminate the second term outright. Why: a geodesic parallel-transports its own tangent, so the derivative of uu along the motion is zero — this is what "straightest possible path" means covariantly.

Step 4. Use Killing's equation on the first term: uνuμ∇νξμ=uνuμ∇(νξμ)=0u^\nu u^\mu\nabla_\nu\xi_\mu = u^\nu u^\mu\nabla_{(\nu}\xi_{\mu)} = 0, since the symmetric product uνuμu^\nu u^\mu projects onto the symmetric part of ∇νξμ\nabla_\nu\xi_\mu, which vanishes. Hence d(ξμuμ)/dτ=0\mathrm{d}(\xi_\mu u^\mu)/\mathrm{d}\tau=0. Why: only the symmetric part of ∇ξ\nabla\xi survives contraction with uνuμu^\nu u^\mu, and that is exactly what Killing's equation zeroes — so the charge is conserved. *Specializing to Schwarzschild's ∂t\partial_t and ∂ϕ\partial_\phi returns the conserved energy EE and angular momentum LL of orbit theory.*

Limiting cases

The Lie derivative simplifies sharply in special cases, and the Killing structure constrains the geometry it lives on:

CaseResultMeaning
Scalar fieldLXf=Xμ∂μf=X(f)\mathcal{L}_X f = X^\mu\partial_\mu f = X(f)Just the directional derivative — no connection, no bracket subtlety; it agrees with ∇Xf\nabla_X f.
Vector fieldLXY=[X,Y]\mathcal{L}_X Y = [X,Y]The Lie bracket; antisymmetric and Jacobi-satisfying, making vector fields a Lie algebra.
Commuting fields[X,Y]=0  ⟺  ϕsX∘ϕtY=ϕtY∘ϕsX[X,Y]=0 \iff \phi^X_s\circ\phi^Y_t=\phi^Y_t\circ\phi^X_sFlows commute exactly when the bracket vanishes; such fields complete to a coordinate basis.
Maximally symmetric space# Killing vectors=12n(n+1)\#\,\text{Killing vectors}=\tfrac12 n(n+1)In nn dimensions: nn translations plus 12n(n−1)\tfrac12 n(n-1) rotations; 1010 for 4D Minkowski (Poincaré).
Conformal Killing limitLξgμν=2 ω gμν\mathcal{L}_\xi g_{\mu\nu}=2\,\omega\,g_{\mu\nu}Metric preserved only up to scale, ω=1n∇λξλ\omega=\tfrac1n\nabla_\lambda\xi^\lambda; true Killing needs ω=0\omega=0.

Worked examples

1 · The three Killing vectors of the round 2-sphere (fully worked)

Take the unit-radius-aa sphere with metric ds2=a2(dθ2+sin⁡2θ dϕ2)\mathrm{d}s^2=a^2(\mathrm{d}\theta^2+\sin^2\theta\,\mathrm{d}\phi^2). Verify that ξ=∂ϕ\xi=\partial_\phi solves Killing's equation, and count how many independent Killing vectors the sphere admits.

Step 1. Write the candidate in components. With coordinates (θ,ϕ)(\theta,\phi), ξμ=(0,1)\xi^\mu=(0,1) and the metric gθθ=a2, gϕϕ=a2sin⁡2θg_{\theta\theta}=a^2,\ g_{\phi\phi}=a^2\sin^2\theta give the lowered field ξμ=gμνξν=(0, a2sin⁡2θ)\xi_\mu=g_{\mu\nu}\xi^\nu=(0,\ a^2\sin^2\theta). The nonzero Christoffel symbols are Γϕϕθ=−sin⁡θcos⁡θ\Gamma^\theta_{\phi\phi}=-\sin\theta\cos\theta and Γθϕϕ=Γϕθϕ=cot⁡θ\Gamma^\phi_{\theta\phi}=\Gamma^\phi_{\phi\theta}=\cot\theta.

Step 2. Test the three independent components of ∇μξν+∇νξμ=∂μξν+∂νξμ−2Γμνλξλ\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=\partial_\mu\xi_\nu+\partial_\nu\xi_\mu-2\Gamma^\lambda_{\mu\nu}\xi_\lambda. The (θθ)(\theta\theta) component is 2∂θξθ−2Γθθλξλ=02\partial_\theta\xi_\theta-2\Gamma^\lambda_{\theta\theta}\xi_\lambda=\mathbf{0} (since ξθ=0, Γθθϕ=0\xi_\theta=0,\ \Gamma^\phi_{\theta\theta}=0), and (ϕϕ)(\phi\phi) is 2∂ϕξϕ−2Γϕϕθξθ=02\partial_\phi\xi_\phi-2\Gamma^\theta_{\phi\phi}\xi_\theta=\mathbf{0} (since ξθ=0\xi_\theta=0).

Step 3. The mixed component is the real test: (θϕ)(\theta\phi) gives ∂θξϕ+∂ϕξθ−2Γθϕλξλ=∂θ(a2sin⁡2θ)−2Γθϕϕξϕ=2a2sin⁡θcos⁡θ−2cot⁡θ⋅a2sin⁡2θ=2a2sin⁡θcos⁡θ−2a2sin⁡θcos⁡θ=0\partial_\theta\xi_\phi+\partial_\phi\xi_\theta-2\Gamma^\lambda_{\theta\phi}\xi_\lambda=\partial_\theta(a^2\sin^2\theta)-2\Gamma^\phi_{\theta\phi}\xi_\phi=2a^2\sin\theta\cos\theta-2\cot\theta\cdot a^2\sin^2\theta=2a^2\sin\theta\cos\theta-2a^2\sin\theta\cos\theta=\mathbf{0}. All three components vanish, so ∂ϕ\partial_\phi is a Killing vector — the manifest one from the metric's ϕ\phi-independence.

Step 4. Count the rest. As a maximally symmetric 2-space the sphere admits 12n(n+1)=12⋅2⋅3=3\tfrac12 n(n+1)=\tfrac12\cdot 2\cdot 3=\mathbf{3} Killing vectors: the generators of SO(3)SO(3), rotations about the xx, yy, and zz axes. Only ∂ϕ\partial_\phi is aligned with a coordinate; the other two are "hidden" rotations, and there are no translational Killing vectors because the sphere is curved. Result: ∂ϕ\partial_\phi is Killing, and S2S^2 has exactly 33 Killing vectors.

2 · Conserved charges and the photon capture radius in Schwarzschild

In Schwarzschild spacetime ds2=−(1−2GM/r) dt2+(1−2GM/r)−1dr2+r2(dθ2+sin⁡2θ dϕ2)\mathrm{d}s^2=-(1-2GM/r)\,\mathrm{d}t^2+(1-2GM/r)^{-1}\mathrm{d}r^2+r^2(\mathrm{d}\theta^2+\sin^2\theta\,\mathrm{d}\phi^2), use the Killing vectors to build the conserved quantities and find the critical impact parameter below which light is captured.

Step 1. The metric is independent of tt and ϕ\phi, so ξ(t)=∂t\xi_{(t)}=\partial_t and ξ(ϕ)=∂ϕ\xi_{(\phi)}=\partial_\phi are Killing. Their conserved charges for a geodesic of tangent uμ=x˙μu^\mu=\dot x^\mu are E=−ξ(t)μuμ=(1−2GM/r) t˙E=-\xi_{(t)\mu}u^\mu=(1-2GM/r)\,\dot t and L=ξ(ϕ)μuμ=r2sin⁡2θ ϕ˙L=\xi_{(\phi)\mu}u^\mu=r^2\sin^2\theta\,\dot\phi. Working in the equatorial plane θ=π/2\theta=\pi/2, these are the conserved energy and angular momentum per unit mass.

Step 2. For a photon the null condition gμνuμuν=0g_{\mu\nu}u^\mu u^\nu=0 combines with the two constants to give r˙2=E2−L2r2(1−2GMr)\dot r^2=E^2-\dfrac{L^2}{r^2}\Big(1-\dfrac{2GM}{r}\Big). The Killing constants have removed the tt and ϕ\phi equations entirely, leaving a one-dimensional radial problem with effective potential V(r)=L2r2(1−2GMr)V(r)=\dfrac{L^2}{r^2}\Big(1-\dfrac{2GM}{r}\Big).

Step 3. Capture versus escape is decided by the peak of V(r)V(r). Setting V′(r)=0V'(r)=0 gives the photon sphere rph=3GM/c2r_{\rm ph}=3GM/c^2, where the maximum equals E2E^2 when the impact parameter b≡L/Eb\equiv L/E reaches bcrit=rph/1−2GM/(c2rph)b_{\rm crit}=r_{\rm ph}/\sqrt{1-2GM/(c^2 r_{\rm ph})}.

Step 4. Plug in rph=3GM/c2r_{\rm ph}=3GM/c^2: the factor 1−2GM/(c2rph)=1−23=131-2GM/(c^2 r_{\rm ph})=1-\tfrac23=\tfrac13, so bcrit=3GM/c21/3=33 GMc2≈5.196 GMc2b_{\rm crit}=\dfrac{3GM/c^2}{\sqrt{1/3}}=3\sqrt{3}\,\dfrac{GM}{c^2}\approx \mathbf{5.196}\,\dfrac{GM}{c^2}. Light with b<bcritb<b_{\rm crit} spirals in and is captured; b>bcritb>b_{\rm crit} escapes. For Sgr A* (GM/c2≈6.0×109 mGM/c^2\approx 6.0\times10^9\ \mathrm{m}) this is bcrit≈3.1×1010 mb_{\rm crit}\approx 3.1\times10^{10}\ \mathrm{m} — the shadow radius the Event Horizon Telescope images, set by exactly this Killing-derived number. Result: bcrit=33 GM/c2b_{\rm crit}=3\sqrt{3}\,GM/c^2, straight from the two Killing vectors.

Physics of the portrait

The Lie derivative and its bracket are the working end of an idea Sophus Lie spent his life building — that continuous symmetry is itself a mathematical object, generated by vector fields and organized by their commutators.

Origins · 1 minSophus Lie and how this was discovered
SL
Sophus Lie
1842–1899 · founder of continuous transformation groups and the geometry of flows

Marius Sophus Lie was a Norwegian mathematician who set out to do for differential equations what Galois had done for polynomials: understand them through their symmetries. Working alongside Felix Klein around 1870, he conceived of continuous groups of transformations — now called Lie groups — captured infinitesimally by the vector fields generating their flows, with the commutator of two generators encoding the group's structure. That commutator is the Lie bracket, the flow-derivative it controls is the Lie derivative, and the algebra of generators is the Lie algebra; his three-volume Theorie der Transformationsgruppen (1888–1893, with Friedrich Engel) founded the whole subject. The physical half of this article rests on a contemporary, Wilhelm Killing, who in the early 1890s wrote down the equation for the vector fields whose flows preserve a metric — the isometry generators now bearing his name. When Emmy Noether tied continuous symmetry to conservation in 1918, Lie's flows and Killing's equation became the language in which "why is energy conserved?" finally had an answer.

History
1842
— Sophus Lie is born in Nordfjordeid, Norway
1870
— Lie, working with Felix Klein, begins the theory of continuous transformation groups
1888
— Lie and Engel publish the first volume of Theorie der Transformationsgruppen
1892
— Wilhelm Killing writes the equation for the vector fields generating isometries
1918
— Emmy Noether links continuous symmetries to conserved quantities, giving Killing vectors their dynamical meaning

Going further

Texts. Nakahara, Geometry, Topology and Physics, 2nd ed., Ch. 5 (Lie derivatives, flows, isometries) · Stone & Goldbart, Mathematics for Physics (CUP 2009), differential-geometry chapters · Arfken, Weber & Harris, Mathematical Methods for Physicists, 7th ed. · Wald, General Relativity (1984), App. C and the symmetries chapter · Weinberg, Gravitation and Cosmology, Ch. 13 (Killing vectors and maximally symmetric spaces) · Riley, Hobson & Bence, Mathematical Methods for Physics and Engineering, 3rd ed. Key results to keep in mind: the coordinate formula (5) remains valid when ∂\partial is replaced by ∇\nabla for any torsion-free connection; Killing's equation ∇(μξν)=0\nabla_{(\mu}\xi_{\nu)}=0 follows from Lξg=0\mathcal{L}_\xi g=0 by metric compatibility; the maximal Killing-vector count is 12n(n+1)\tfrac12 n(n+1) — 1010 for 4D Minkowski, 33 for S2S^2, 66 for 3-space; and the Schwarzschild conserved quantities EE, LL and the photon capture radius bcrit=33 GM/c2b_{\rm crit}=3\sqrt3\,GM/c^2 are both worked in full above.

Prerequisites. Manifolds and Tangent Spaces · Covariant and Contravariant Tensors · The Metric Tensor · The Covariant Derivative · Christoffel Symbols and Affine Connections. Leads to. Parallel Transport and Geodesics · Exterior Calculus and the Exterior Derivative · Differential Forms.

Ask the tutor: Ask me to solve Killing's equation for a metric of your choice, or to trace how a single Killing vector turns into a conserved energy, momentum, or shadow radius.

What is the difference between the Lie derivative and the covariant derivative?

Both differentiate tensor fields, but they need different inputs. The covariant derivative ∇XT\nabla_X T asks how TT changes relative to parallel transport, which requires a connection — an extra rule for carrying a tensor between neighbouring tangent spaces. The Lie derivative LXT\mathcal{L}_X T asks how TT changes as the flow of XX drags it, using only the flow itself, so it needs no connection and no metric. A practical tell: ∇XT\nabla_X T at a point depends on XX only at that point, while LXT\mathcal{L}_X T depends on XX throughout a neighbourhood, because it differentiates XX too.

Why is the Lie derivative not tensorial in X?

A tensorial (pointwise) dependence would mean LfXT=f LXT\mathcal{L}_{fX}T=f\,\mathcal{L}_X T for every function ff — the value at a point set by XX at that point alone. But the coordinate formula contains derivatives ∂λXμ\partial_\lambda X^\mu, so multiplying XX by ff also differentiates ff. The cleanest instance is LfXY=[fX,Y]=f[X,Y]−(Yf)X\mathcal{L}_{fX}Y=[fX,Y]=f[X,Y]-(Yf)X: the extra −(Yf)X-(Yf)X is the non-tensorial term. This is exactly why the Lie derivative is not a "directional derivative of TT" in the way ∇XT\nabla_X T is — it measures a flow effect, which is inherently non-local in XX.

How do I actually find the Killing vectors of a given metric?

Two routes. The quick one: scan the metric for coordinates it does not depend on. If gμνg_{\mu\nu} is independent of xax^a, then ∂a\partial_a is immediately a Killing vector — this catches ∂t\partial_t and ∂ϕ\partial_\phi in Schwarzschild. The complete one: write ξμ\xi_\mu with unknown component functions and solve ∇μξν+∇νξμ=0\nabla_\mu\xi_\nu+\nabla_\nu\xi_\mu=0 as coupled first-order PDEs. This also finds the "hidden" symmetries, like the two non-ϕ\phi rotations of the sphere, and its integrability conditions cap the count at 12n(n+1)\tfrac12 n(n+1).

What is a conformal Killing vector, and how does it differ?

A conformal Killing vector preserves the metric only up to a local rescaling: Lξgμν=2ω gμν\mathcal{L}_\xi g_{\mu\nu}=2\omega\,g_{\mu\nu} with ω=1n∇λξλ\omega=\tfrac1n\nabla_\lambda\xi^\lambda, rather than the strict Lξg=0\mathcal{L}_\xi g=0 of a true Killing vector. Its flow preserves angles and the light-cone structure but not lengths. These generate the conformal symmetries central to massless field theory, 2D CFT, and null geodesics — a photon's path depends only on the conformal class of the metric, so conformal Killing vectors give conserved quantities for light even when full isometries are absent. A true Killing vector is the case ω=0\omega=0.

Why does each Killing vector produce a conserved quantity?

Because a Killing vector is a symmetry of the geometry, and symmetry means something is not being pushed by the dynamics. Concretely, take the scalar ξμuμ\xi_\mu u^\mu along a geodesic and differentiate: one term dies by the geodesic equation (the tangent is parallel-transported), and the other dies by Killing's equation (the symmetric uμuνu^\mu u^\nu contracts against a derivative whose symmetric part vanishes). The result is d(ξμuμ)/dτ=0\mathrm{d}(\xi_\mu u^\mu)/\mathrm{d}\tau=0. This is the geometric face of Noether's theorem: the isometry generated by ξ\xi is a continuous symmetry, and its Noether charge is precisely ξμuμ\xi_\mu u^\mu.

What is Cartan's magic formula good for?

Cartan's formula LXα=d(ιXα)+ιX(dα)\mathcal{L}_X\alpha=\mathrm{d}(\iota_X\alpha)+\iota_X(\mathrm{d}\alpha) computes the Lie derivative of any differential form using only the exterior derivative and interior product — no coordinates, no index gymnastics. It is the fastest route to transport theorems: on a volume form it gives Lu(dV)=(∇⋅u) dV\mathcal{L}_u(\mathrm{d}V)=(\nabla\cdot u)\,\mathrm{d}V, the divergence as the fractional rate of volume change under the flow. On a closed form (dα=0\mathrm{d}\alpha=0) it shows the Lie derivative is exact — the mechanism behind conserved fluxes and Helmholtz's frozen-in vorticity. It is the workhorse of geometric mechanics and the Reynolds transport theorem.

Before this

Where this leads