Covariant Derivative

In flat space with Cartesian axes the partial derivative of a vector field is already a tensor — the basis vectors never move, so differentiating the components captures the field's entire change. That is a pleasant accident, and the exception rather than the rule. Switch to polar or spherical coordinates, or step onto a genuinely curved manifold, and the coordinate basis vectors turn from point to point. Now the change of a vector field splits in two: the change of its components, which the partial derivative sees, and the change of the basis carrying them, which it does not. So ∂μVν\partial_\mu V^\nu no longer transforms as a tensor — its value depends on the coordinate grid you happened to draw. The covariant derivative ∇μ\nabla_\mu is the repair: it adds one connection term per index, built from the coefficients Γλμν\Gamma^\lambda{}_{\mu\nu} that encode exactly how the basis rotates, and the non-tensorial defect cancels precisely. What survives is a genuine tensor measuring the change of a field relative to parallel transport — the change you cannot blame on the grid. On a metric manifold one connection is distinguished, singled out by demanding it be torsion-free and metric-compatible, and its coefficients are the Christoffel symbols. This single operator is the workhorse of all gravitation: geodesics, conservation laws, the geodesic deviation we feel as tides, and the curvature that is gravity itself all descend from it. Everything downstream in relativity is written with a ∇\nabla.

Key takeaways
  • The covariant derivative ∇μ\nabla_\mu generalizes ∂μ\partial_\mu so that differentiating a tensor field yields another tensor field, compensating for the coordinate-dependence of the partial derivative in curved spaces or curvilinear coordinates.
  • It adds one connection term per index: ∇μVν=∂μVν+ΓνμλVλ\nabla_\mu V^\nu=\partial_\mu V^\nu+\Gamma^\nu{}_{\mu\lambda}V^\lambda for an upper index, with a minus sign −Γλμνωλ-\Gamma^\lambda{}_{\mu\nu}\omega_\lambda for a lower one.
  • On a metric manifold two demands — zero torsion and metric compatibility ∇σgμν=0\nabla_\sigma g_{\mu\nu}=0 — fix the connection uniquely as the Christoffel symbols Γλμν=12gλσ(∂μgσν+∂νgσμ−∂σgμν)\Gamma^\lambda{}_{\mu\nu}=\tfrac12 g^{\lambda\sigma}(\partial_\mu g_{\sigma\nu}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}).
  • A vanishing covariant derivative along a curve is parallel transport, and a curve that parallel-transports its own tangent is a geodesic — this is where the equation of motion comes from.
  • The failure of two covariant derivatives to commute is the Riemann curvature tensor, [∇μ,∇ν]Vρ=RρσμνVσ[\nabla_\mu,\nabla_\nu]V^\rho=R^\rho{}_{\sigma\mu\nu}V^\sigma — the gateway from kinematics to the geometry of gravity.
  • Key idea: to differentiate a tensor covariantly, take the partial derivative of each component and add one connection term +Γνμλ+\Gamma^\nu{}_{\mu\lambda} for every upper index and subtract one −Γλμν-\Gamma^\lambda{}_{\mu\nu} for every lower index; the connection cancels the non-tensorial part of ∂\partial, leaving a true tensor that measures change relative to parallel transport.
The core insight

To differentiate a tensor covariantly, take the ordinary partial derivative of each component and then add one connection term +Γνμλ+\Gamma^{\nu}{}_{\mu\lambda} for every upper index and subtract one −Γλμν-\Gamma^{\lambda}{}_{\mu\nu} for every lower index. The connection terms cancel the non-tensorial part of ∂\partial, leaving a genuine tensor that measures change relative to parallel transport.

Where this lives today

The covariant derivative is the differential operator inside every line of working relativity. Local energy-momentum conservation is written ∇μTμν=0\nabla_\mu T^{\mu\nu}=0, an identity that follows from the contracted Bianchi identity ∇μGμν=0\nabla_\mu G^{\mu\nu}=0 and is the consistency condition of the Einstein field equations — the reason matter cannot be specified independently of the geometry it sources. In numerical relativity, the codes that model merging black holes for LIGO and Virgo evolve the spatial metric and its covariant derivatives on a grid, recomputing the Christoffel symbols at every point and time step; the waveforms compared against detector data are solutions of covariant-derivative equations. The tidal signal a gravitational-wave detector is built to sense — the relative acceleration of nearby freely-falling masses — is the second covariant derivative of their separation, D2Sμdτ2=RμνρσuνuρSσ\tfrac{D^2 S^\mu}{\mathrm{d}\tau^2}=R^\mu{}_{\nu\rho\sigma}u^\nu u^\rho S^\sigma. Symmetries and their conserved charges come from Killing's equation ∇(μξν)=0\nabla_{(\mu}\xi_{\nu)}=0, which turns the time-translation symmetry of Schwarzschild into conserved orbital energy and gives Kerr its angular momentum. And every matter field is placed in a gravitational background by the minimal-coupling rule ∂→∇\partial\to\nabla — electromagnetism as ∇μFμν=μ0Jν\nabla_\mu F^{\mu\nu}=\mu_0 J^\nu, scalars, relativistic fluids alike — so that ∇\nabla is the single operation coupling all of physics to gravity.

The physical picture

Walk across the surface of the Earth carrying an arrow that you keep "as straight as possible," never deliberately twisting it. In your local north–east frame the arrow's components change even though you did nothing to it, because the frame itself swings around as you follow the curved ground. The partial derivative is an over-eager accountant who blames you for that swing; the connection terms are the exchange-rate correction that subtracts it. Once subtracted, "no covariant change" means "genuinely parallel-transported." The most dramatic consequence is that parallel transport around a closed loop does not bring the arrow home unchanged: on a curved surface it returns rotated, and the rotation angle equals the area the loop encloses, in units of the radius squared. That single fact is the local fingerprint of curvature, and it is what the covariant derivative was built to detect.

Figure
Parallel transport around a closed geodesic triangle on a sphere brings the vector back rotated by the enclosed solid angle. For the octant N→A→B→N, each corner is a right angle, the enclosed area is π/2 of a unit sphere, and the vector returns turned by exactly 90° — while the same loop drawn in flat space returns it unchanged
Parallel transport around a loop rotates the vector by the enclosed area N A B start return 90° loop N → A → B → N each side a geodesic (great circle) three right-angle corners rotation Δα = enclosed area Δα = Ω = A/R² = π/2 sr = 90° = ∮ curvature (Gauss–Bonnet) flat space: same loop, Δα = 0 the holonomy that ∇ measures

That rotation is bookkept, term by term, by the connection Γ\Gamma — how the coordinate basis vectors turn. In flat space with Cartesian axes the basis is rigid, every Γ\Gamma vanishes, and ∇\nabla collapses to ∂\partial: a loop returns the arrow unchanged. The moment the basis rotates — in polar coordinates on a flat plane, or on any curved surface — a moving-basis correction appears, and ∇μ=∂μ+Γ\nabla_\mu=\partial_\mu+\Gamma is exactly the operator that subtracts it. The second figure makes the cancellation concrete: a constant, physically unchanging field acquires changing components in polar coordinates purely because the basis swings, and the connection term erases that spurious change to leave ∇=0\nabla=0.

Figure
On a flat plane in polar coordinates a constant vector field V = ê_x still has coordinate components that change from point to point, because the polar basis rotates. At the point r = 2, ϕ = 90° the partial derivative reads ∂_r V^ϕ = +0.25, but the connection term Γ^ϕ_rϕ V^ϕ = (1/r)(−0.5) = −0.25 cancels it exactly, so ∇_r V^ϕ = 0 — the field is correctly seen as unchanging
∇ = ∂ + Γ : the connection cancels the turning of the basis origin P : r = 2, ϕ = 90° ê_r ê_ϕ V = ê_x same V, larger r ledger at P (r = 2, ϕ = 90°) V = ê_x ⇒ V^r = 0, V^ϕ = −0.5 Γ^ϕ_rϕ = 1/r = 0.5 ∇_r V^ϕ = ∂_r V^ϕ + Γ^ϕ_rϕ V^ϕ ∂_r V^ϕ = +0.25 Γ^ϕ_rϕ V^ϕ = (0.5)(−0.5) = −0.25 ∇_r V^ϕ = 0.25 − 0.25 = 0 ✓

Two covariant derivatives generally fail to commute, and that failure — carrying a vector around an infinitesimal loop and reading off how much it turned — is precisely the Riemann tensor, [∇μ,∇ν]Vρ=RρσμνVσ[\nabla_\mu,\nabla_\nu]V^\rho=R^\rho{}_{\sigma\mu\nu}V^\sigma. Because ∇\nabla manufactures tensors, any law written with it holds in every coordinate system once it holds in one — the operational meaning of general covariance. This is why the equations of motion of special relativity are promoted to curved spacetime by the "comma-goes-to-semicolon" rule ∂→∇\partial\to\nabla: minimal coupling to gravity is nothing but replacing every partial derivative by a covariant one.

The core equations

The covariant derivative of a contravariant vector is the partial derivative plus one connection term contracting the connection's upper index against the vector. Neither piece is a tensor, but their sum is — that cancellation is the whole design brief, established in Derivation 1.

∇μVν=∂μVν+ΓνμλVλ(1)\nabla_\mu V^\nu=\partial_\mu V^\nu+\Gamma^\nu{}_{\mu\lambda}V^\lambda \qquad (1)

where: ∇μ\nabla_\mu — covariant derivative in the xμx^\mu direction; VνV^\nu — contravariant vector components; ∂μ=∂/∂xμ\partial_\mu=\partial/\partial x^\mu; Γνμλ\Gamma^\nu{}_{\mu\lambda} — connection coefficients (Christoffel symbols for Levi-Civita). Reading: the free indices make ∇μVν\nabla_\mu V^\nu a genuine (1,1)(1,1) tensor. For a one-form the connection term flips sign and contracts the lower index. This is not a convention — it is forced by the Leibniz rule on the scalar contraction ωνVν\omega_\nu V^\nu (Derivation 1).

∇μων=∂μων−Γλμνωλ(2)\nabla_\mu \omega_\nu=\partial_\mu \omega_\nu-\Gamma^\lambda{}_{\mu\nu}\omega_\lambda \qquad (2)

where: ων\omega_\nu — one-form components. Reading: an upper index brings +Γ+\Gamma, a lower index brings −Γ-\Gamma. The two rules combine into a single master formula: one +Γ+\Gamma per upper index and one −Γ-\Gamma per lower index. Everything else in the article is a special case of this line.

∇μTν1⋯ρ1⋯=∂μTν1⋯ρ1⋯+∑iΓνiμλT⋯λ⋯ρ1⋯−∑jΓλμρjTν1⋯⋯λ⋯(3)\nabla_\mu T^{\nu_1\cdots}{}_{\rho_1\cdots}=\partial_\mu T^{\nu_1\cdots}{}_{\rho_1\cdots}+\sum_i\Gamma^{\nu_i}{}_{\mu\lambda}T^{\cdots\lambda\cdots}{}_{\rho_1\cdots}-\sum_j\Gamma^{\lambda}{}_{\mu\rho_j}T^{\nu_1\cdots}{}_{\cdots\lambda\cdots} \qquad (3)

where: Tν1⋯ρ1⋯T^{\nu_1\cdots}{}_{\rho_1\cdots} — an arbitrary (k,l)(k,l) tensor; the sums run over its upper and lower indices. Reading: ∇\nabla raises the covariant rank by one and preserves tensor character exactly. On a metric manifold there is one preferred connection, its coefficients the Christoffel symbols — an explicit combination of first derivatives of the metric. This is what "the" covariant derivative means in general relativity.

Γλμν=12gλσ(∂μgσν+∂νgσμ−∂σgμν)(4)\Gamma^\lambda{}_{\mu\nu}=\tfrac{1}{2}g^{\lambda\sigma}\left(\partial_\mu g_{\sigma\nu}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}\right) \qquad (4)

where: gμνg_{\mu\nu} — metric; gλσg^{\lambda\sigma} — inverse metric; Γλμν\Gamma^\lambda{}_{\mu\nu} — Christoffel symbols, symmetric in μν\mu\nu. Reading: the connection is built entirely from the metric, so the metric fixes how vectors are transported. The property that selects this connection is metric compatibility: the covariant derivative of the metric vanishes identically, so ∇\nabla commutes with raising and lowering indices.

∇σgμν=0(5)\nabla_\sigma g_{\mu\nu}=0 \qquad (5)

where: gμνg_{\mu\nu} — metric. Reading: because ∇g=0\nabla g=0, one may slide the metric through a covariant derivative, e.g. ∇μVν=gνρ∇μVρ\nabla_\mu V_\nu=g_{\nu\rho}\nabla_\mu V^\rho. The second defining property is zero torsion — the connection's antisymmetric part vanishes, making Γ\Gamma symmetric in its lower pair. The torsion tensor measures that antisymmetry.

Tλμν=Γλμν−Γλνμ=2Γλ[μν](6)T^\lambda{}_{\mu\nu}=\Gamma^\lambda{}_{\mu\nu}-\Gamma^\lambda{}_{\nu\mu}=2\Gamma^\lambda{}_{[\mu\nu]} \qquad (6)

where: TλμνT^\lambda{}_{\mu\nu} — torsion (zero for Levi-Civita); [μν][\mu\nu] — antisymmetrization. Reading: unlike Γ\Gamma, the torsion is a genuine tensor, since the inhomogeneous term is symmetric and cancels in the difference. Restricting ∇\nabla to a curve gives the directional (absolute) derivative; demanding a curve parallel-transport its own tangent produces the geodesic equation — where the free-fall equation of motion is born.

DVμdτ≡uν∇νVμ=dVμdτ+ΓμνρuνVρ,uν∇νuμ=0 ⇒ d2xμdτ2+Γμνρdxνdτdxρdτ=0(7)\frac{DV^\mu}{\mathrm{d}\tau}\equiv u^\nu\nabla_\nu V^\mu=\frac{\mathrm{d}V^\mu}{\mathrm{d}\tau}+\Gamma^\mu{}_{\nu\rho}u^\nu V^\rho,\qquad u^\nu\nabla_\nu u^\mu=0\ \Rightarrow\ \frac{\mathrm{d}^2 x^\mu}{\mathrm{d}\tau^2}+\Gamma^\mu{}_{\nu\rho}\frac{\mathrm{d}x^\nu}{\mathrm{d}\tau}\frac{\mathrm{d}x^\rho}{\mathrm{d}\tau}=0 \qquad (7)

where: τ\tau — affine parameter; uμ=dxμ/dτu^\mu=\mathrm{d}x^\mu/\mathrm{d}\tau — tangent; DVμ/dτDV^\mu/\mathrm{d}\tau — absolute derivative along the curve. Reading: a geodesic is a curve whose tangent has zero covariant derivative along itself — "as straight as the geometry allows." Finally, the commutator of two covariant derivatives on a vector is the Riemann curvature tensor (plus a torsion term, zero for Levi-Civita). This is the deepest payoff: ∇\nabla turns curvature from a picture into an operator.

[∇μ,∇ν]Vρ=RρσμνVσ−Tλμν∇λVρ(8)[\nabla_\mu,\nabla_\nu]V^\rho=R^\rho{}_{\sigma\mu\nu}V^\sigma-T^\lambda{}_{\mu\nu}\nabla_\lambda V^\rho \qquad (8)

where: RρσμνR^\rho{}_{\sigma\mu\nu} — Riemann curvature tensor; TλμνT^\lambda{}_{\mu\nu} — torsion (zero for Levi-Civita). Reading: on a flat manifold covariant derivatives commute; their failure to commute is exactly the local curvature (Derivation 3).

Derivation AWhy the connection term is forced, and the minus sign for one-forms

The goal is to show that adding Γ\Gamma to the partial derivative is not an optional patch but the unique fix demanded by tensoriality, and that the same requirement — applied through the Leibniz rule to a scalar contraction — forces the one-form connection term to carry the opposite sign. We transform the partial derivative and the connection separately, watch their non-tensorial pieces cancel in the sum, then use ∇\nabla of a scalar to pin the one-form rule with no further freedom.

Step 1. Transform the partial derivative to a primed frame: differentiating Vν′=(∂xν′/∂xν)VνV^{\nu'}=(\partial x^{\nu'}/\partial x^\nu)V^\nu gives ∂μ′Vν′=∂xμ∂xμ′∂xν′∂xν∂μVν+∂xμ∂xμ′∂2xν′∂xμ∂xνVν\partial_{\mu'}V^{\nu'}=\frac{\partial x^\mu}{\partial x^{\mu'}}\frac{\partial x^{\nu'}}{\partial x^\nu}\partial_\mu V^\nu+\frac{\partial x^\mu}{\partial x^{\mu'}}\frac{\partial^2 x^{\nu'}}{\partial x^\mu\partial x^\nu}V^\nu. Why: the first term is the clean tensor transformation; the second — a derivative of the Jacobian — is the non-tensorial defect that spoils ∂μVν\partial_\mu V^\nu.

Step 2. Transform the connection, which is defined to carry an inhomogeneous piece Γν′μ′λ′=∂xν′∂xν∂xμ∂xμ′∂xλ∂xλ′Γνμλ−∂xμ∂xμ′∂xλ∂xλ′∂2xν′∂xμ∂xλ\Gamma^{\nu'}{}_{\mu'\lambda'}=\frac{\partial x^{\nu'}}{\partial x^\nu}\frac{\partial x^\mu}{\partial x^{\mu'}}\frac{\partial x^\lambda}{\partial x^{\lambda'}}\Gamma^\nu{}_{\mu\lambda}-\frac{\partial x^\mu}{\partial x^{\mu'}}\frac{\partial x^\lambda}{\partial x^{\lambda'}}\frac{\partial^2 x^{\nu'}}{\partial x^\mu\partial x^\lambda}. Why: its second-derivative term has precisely the opposite sign of the Step 1 defect.

Step 3. Add them. In ∇μ′Vν′=∂μ′Vν′+Γν′μ′λ′Vλ′\nabla_{\mu'}V^{\nu'}=\partial_{\mu'}V^{\nu'}+\Gamma^{\nu'}{}_{\mu'\lambda'}V^{\lambda'} the two inhomogeneous terms cancel identically, leaving ∇μ′Vν′=∂xμ∂xμ′∂xν′∂xν∇μVν\nabla_{\mu'}V^{\nu'}=\frac{\partial x^\mu}{\partial x^{\mu'}}\frac{\partial x^{\nu'}}{\partial x^\nu}\nabla_\mu V^\nu. Why: this exact cancellation is the entire reason the connection exists — it makes ∇μVν\nabla_\mu V^\nu a genuine (1,1)(1,1) tensor.

Step 4. Fix the one-form sign from Leibniz. Since ωνVν\omega_\nu V^\nu is a scalar, ∇μ(ωνVν)=∂μ(ωνVν)\nabla_\mu(\omega_\nu V^\nu)=\partial_\mu(\omega_\nu V^\nu); expanding gives (∇μων)Vν+ων(∂μVν+ΓνμλVλ)=(∂μων)Vν+ων∂μVν(\nabla_\mu\omega_\nu)V^\nu+\omega_\nu(\partial_\mu V^\nu+\Gamma^\nu{}_{\mu\lambda}V^\lambda)=(\partial_\mu\omega_\nu)V^\nu+\omega_\nu\partial_\mu V^\nu, and matching the coefficient of arbitrary VνV^\nu forces ∇μων=∂μων−Γλμνωλ\nabla_\mu\omega_\nu=\partial_\mu\omega_\nu-\Gamma^\lambda{}_{\mu\nu}\omega_\lambda. Why: one scalar identity across all vectors leaves no freedom — the minus sign is mandatory. *Verified: reproduces Carroll (3.1)–(3.6) and Wald §3.1; recompute by symmetry of the second derivatives and by matching coefficients of VνV^\nu in the scalar Leibniz identity.*

Go deeper · Derivation B · 1 minThe unique metric-compatible, torsion-free connection
Derivation BThe unique metric-compatible, torsion-free connection

The goal is to prove that the two structural demands of general relativity — that ∇\nabla annihilate the metric (∇σgμν=0\nabla_\sigma g_{\mu\nu}=0) and that the connection have no torsion (Γλμν=Γλνμ\Gamma^\lambda{}_{\mu\nu}=\Gamma^\lambda{}_{\nu\mu}) — do not merely constrain the connection but determine it completely, giving the Christoffel symbols as the one and only solution. The strategy is a classic index gymnastic: write metric compatibility three times with the indices cyclically permuted, add two and subtract the third, and use symmetry to collapse six connection terms into one.

Step 1. Write metric compatibility in full, ∇μgνσ=∂μgνσ−Γλμνgλσ−Γλμσgνλ=0\nabla_\mu g_{\nu\sigma}=\partial_\mu g_{\nu\sigma}-\Gamma^\lambda{}_{\mu\nu}g_{\lambda\sigma}-\Gamma^\lambda{}_{\mu\sigma}g_{\nu\lambda}=0, and record the two cyclic permutations ∇νgσμ=0\nabla_\nu g_{\sigma\mu}=0 and ∇σgμν=0\nabla_\sigma g_{\mu\nu}=0. Why: three copies with rotated indices give three linear relations among the same coefficients, enough to solve for them.

Step 2. Add the first two and subtract the third. The partial terms assemble into ∂μgνσ+∂νgσμ−∂σgμν\partial_\mu g_{\nu\sigma}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}, and the six connection terms — using Γλμν=Γλνμ\Gamma^\lambda{}_{\mu\nu}=\Gamma^\lambda{}_{\nu\mu} — cancel in four places and reinforce in two, leaving ∂μgνσ+∂νgσμ−∂σgμν=2Γλμνgλσ\partial_\mu g_{\nu\sigma}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}=2\Gamma^\lambda{}_{\mu\nu}g_{\lambda\sigma}. Why: torsion-freeness is exactly what makes the unwanted terms pair off, isolating a single Γg\Gamma g.

Step 3. Contract with gλσg^{\lambda\sigma} to strip the metric factor: Γλμν=12gλσ(∂μgσν+∂νgσμ−∂σgμν)\Gamma^\lambda{}_{\mu\nu}=\tfrac{1}{2}g^{\lambda\sigma}\left(\partial_\mu g_{\sigma\nu}+\partial_\nu g_{\sigma\mu}-\partial_\sigma g_{\mu\nu}\right). Why: this delivers Γ\Gamma alone; every step was forced, so the connection is unique — the fundamental theorem of Riemannian geometry. *Verified: reproduces Carroll (3.27) and Wald (3.1.30); recompute by symmetrizing the three metric-compatibility equations and confirming that the antisymmetric-in-μν\mu\nu parts drop out.*

Go deeper · Derivation C · 1 minCurvature as the failure of covariant derivatives to commute
Derivation CCurvature as the failure of covariant derivatives to commute

The goal is to show that although a single covariant derivative merely repairs tensoriality, the antisymmetrized double covariant derivative uncovers something new and irreducibly geometric — the Riemann curvature tensor — and to exhibit it explicitly. We compute ∇μ∇νVρ\nabla_\mu\nabla_\nu V^\rho by treating the first derivative as a (1,1)(1,1) tensor and applying the master rule again, then antisymmetrize in μν\mu\nu (torsion vanishing) and watch every non-curvature term cancel, leaving a tensor that acts linearly on VV.

Step 1. Take the first covariant derivative, ∇νVρ=∂νVρ+ΓρνσVσ\nabla_\nu V^\rho=\partial_\nu V^\rho+\Gamma^\rho{}_{\nu\sigma}V^\sigma, which is a (1,1)(1,1) tensor with a lower index ν\nu and an upper index ρ\rho. Why: the result must now be differentiated as a tensor, so the second derivative will carry a −Γ-\Gamma for the ν\nu slot and a +Γ+\Gamma for the ρ\rho slot.

Step 2. Apply the master rule again: ∇μ∇νVρ=∂μ(∇νVρ)−Γλμν∇λVρ+Γρμλ∇νVλ\nabla_\mu\nabla_\nu V^\rho=\partial_\mu(\nabla_\nu V^\rho)-\Gamma^\lambda{}_{\mu\nu}\nabla_\lambda V^\rho+\Gamma^\rho{}_{\mu\lambda}\nabla_\nu V^\lambda. Why: the middle term treats the differentiated lower index ν\nu, the last term treats the upper index ρ\rho — this is just equation (3) applied to a rank-2 object.

Step 3. Antisymmetrize in μν\mu\nu. The ∂μ∂νVρ\partial_\mu\partial_\nu V^\rho terms are symmetric and cancel; the Γλμν∇λVρ\Gamma^\lambda{}_{\mu\nu}\nabla_\lambda V^\rho term is symmetric under μ↔ν\mu\leftrightarrow\nu (torsion-free) and cancels; what survives are the terms in which a derivative or a product hits the connection. Why: only the pieces that are genuinely antisymmetric in μν\mu\nu endure, and those are built from ∂Γ\partial\Gamma and ΓΓ\Gamma\Gamma, not from derivatives of VV.

Step 4. Collect the survivors multiplying VσV^\sigma: [∇μ,∇ν]Vρ=(∂μΓρνσ−∂νΓρμσ+ΓρμλΓλνσ−ΓρνλΓλμσ)Vσ≡RρσμνVσ[\nabla_\mu,\nabla_\nu]V^\rho=\left(\partial_\mu\Gamma^\rho{}_{\nu\sigma}-\partial_\nu\Gamma^\rho{}_{\mu\sigma}+\Gamma^\rho{}_{\mu\lambda}\Gamma^\lambda{}_{\nu\sigma}-\Gamma^\rho{}_{\nu\lambda}\Gamma^\lambda{}_{\mu\sigma}\right)V^\sigma\equiv R^\rho{}_{\sigma\mu\nu}V^\sigma. Why: the bracket is a genuine tensor because the commutator kills the inhomogeneous Γ\Gamma pieces — the Riemann tensor, which vanishes if and only if the manifold is flat. *Verified: matches Carroll (3.112) and Wald (3.2.3); recompute by expanding the nested covariant derivatives and discarding the terms symmetric in μν\mu\nu.*

Limiting cases

The covariant derivative reduces to, or specializes in, several regimes worth committing to memory — each one a check on the general formula.

CaseResultMeaning
Flat space, Cartesian coordinatesΓλμν=0⇒∇μ=∂μ\Gamma^\lambda{}_{\mu\nu}=0\Rightarrow\nabla_\mu=\partial_\muAll connection coefficients vanish, so covariant and partial derivatives coincide everywhere; loops return vectors unchanged.
Acting on a scalar∇μϕ=∂μϕ\nabla_\mu\phi=\partial_\mu\phiA scalar carries no index, so there is no connection term; the gradient of a scalar is already a one-form tensor.
Local inertial frame at a point ppΓλμν(p)=0, ∂σΓλμν(p)≠0\Gamma^\lambda{}_{\mu\nu}(p)=0,\ \partial_\sigma\Gamma^\lambda{}_{\mu\nu}(p)\neq0In Riemann normal coordinates the connection vanishes at one point but not its derivatives — first covariant derivatives equal partials there, yet curvature persists.
Metric-compatible (Levi-Civita) connection∇σgμν=0\nabla_\sigma g_{\mu\nu}=0Covariant differentiation commutes with raising and lowering indices; this, with zero torsion, fixes Γ\Gamma to be the Christoffel symbols.
Flat connection, curvilinear coordinatesΓ≠0\Gamma\neq0 but Rρσμν=0R^\rho{}_{\sigma\mu\nu}=0Polar or spherical coordinates on flat space: the connection is nonzero (the basis rotates) yet curvature is identically zero.

Worked examples

1 · Christoffel symbols of the unit two-sphere (confidence builder)

The unit sphere has metric ds2=dθ2+sin⁡2θ dϕ2\mathrm{d}s^2=\mathrm{d}\theta^2+\sin^2\theta\,\mathrm{d}\phi^2, so gθθ=1g_{\theta\theta}=1, gϕϕ=sin⁡2θg_{\phi\phi}=\sin^2\theta, and gθϕ=0g_{\theta\phi}=0. Find the nonzero Christoffel symbols, evaluate them at θ=45°\theta=45°, and verify metric compatibility there numerically.

Step 1. Only gϕϕ=sin⁡2θg_{\phi\phi}=\sin^2\theta depends on a coordinate, so the sole nonzero metric derivative is ∂θgϕϕ=2sin⁡θcos⁡θ\partial_\theta g_{\phi\phi}=2\sin\theta\cos\theta. Feeding this into the Christoffel formula with gθθ=1g^{\theta\theta}=1, gϕϕ=1/sin⁡2θg^{\phi\phi}=1/\sin^2\theta gives Γθϕϕ=12gθθ(−∂θgϕϕ)=−sin⁡θcos⁡θ\Gamma^\theta{}_{\phi\phi}=\tfrac12 g^{\theta\theta}(-\partial_\theta g_{\phi\phi})=\mathbf{-\sin\theta\cos\theta} and Γϕθϕ=Γϕϕθ=12gϕϕ∂θgϕϕ=cot⁡θ\Gamma^\phi{}_{\theta\phi}=\Gamma^\phi{}_{\phi\theta}=\tfrac12 g^{\phi\phi}\partial_\theta g_{\phi\phi}=\mathbf{\cot\theta}; all others vanish.

Step 2. At θ=45°\theta=45°, sin⁡θ=cos⁡θ=22≈0.7071\sin\theta=\cos\theta=\tfrac{\sqrt2}{2}\approx0.7071, so Γθϕϕ=−(0.7071)(0.7071)=−0.5\Gamma^\theta{}_{\phi\phi}=-(0.7071)(0.7071)=\mathbf{-0.5} and Γϕθϕ=cot⁡45°=1.0\Gamma^\phi{}_{\theta\phi}=\cot45°=\mathbf{1.0} — the entire connection of the round sphere at that latitude.

Step 3. Verify metric compatibility, ∇θgϕϕ=∂θgϕϕ−2Γϕθϕgϕϕ\nabla_\theta g_{\phi\phi}=\partial_\theta g_{\phi\phi}-2\Gamma^\phi{}_{\theta\phi}g_{\phi\phi} (only λ=ϕ\lambda=\phi survives). At θ=45°\theta=45°, ∂θgϕϕ=2(0.7071)(0.7071)=1.0\partial_\theta g_{\phi\phi}=2(0.7071)(0.7071)=\mathbf{1.0} and 2Γϕθϕgϕϕ=2(1.0)(0.5)=1.02\Gamma^\phi{}_{\theta\phi}g_{\phi\phi}=2(1.0)(0.5)=\mathbf{1.0}, so ∇θgϕϕ=1.0−1.0=0\nabla_\theta g_{\phi\phi}=1.0-1.0=\mathbf{0} — covariantly constant, as the Levi-Civita connection demands. Γθϕϕ=−0.5\Gamma^\theta{}_{\phi\phi}=-0.5, Γϕθϕ=1.0\Gamma^\phi{}_{\theta\phi}=1.0 at θ=45°\theta=45°, and ∇θgϕϕ=0\nabla_\theta g_{\phi\phi}=0.

2 · The covariant derivative sees through a rotating basis (fully worked)

On the flat plane in polar coordinates, ds2=dr2+r2dϕ2\mathrm{d}s^2=\mathrm{d}r^2+r^2\mathrm{d}\phi^2, take the constant Cartesian field V=e^xV=\hat e_x — the same arrow everywhere. Show that its polar components change from point to point yet ∇V=0\nabla V=0, evaluating everything at r=2r=2, ϕ=90°\phi=90°.

Step 1. The nonzero Christoffel symbols of the polar plane are Γrϕϕ=−r\Gamma^r{}_{\phi\phi}=-r and Γϕrϕ=Γϕϕr=1/r\Gamma^\phi{}_{r\phi}=\Gamma^\phi{}_{\phi r}=1/r (from gϕϕ=r2g_{\phi\phi}=r^2, as in Example 1). The uniform field e^x=cos⁡ϕ ∂r−(sin⁡ϕ/r) ∂ϕ\hat e_x=\cos\phi\,\partial_r-(\sin\phi/r)\,\partial_\phi has components Vr=cos⁡ϕV^r=\cos\phi, Vϕ=−sin⁡ϕ/rV^\phi=-\sin\phi/r that vary with position even though the physical arrow is fixed.

Step 2. At r=2r=2, ϕ=90°\phi=90°: Vr=cos⁡90°=0V^r=\cos90°=\mathbf{0} and Vϕ=−sin⁡90°/2=−0.5V^\phi=-\sin90°/2=\mathbf{-0.5}. The offending partial is ∂rVϕ=∂r(−sin⁡ϕ/r)=sin⁡ϕ/r2=sin⁡90°/4=+0.25\partial_r V^\phi=\partial_r(-\sin\phi/r)=\sin\phi/r^2=\sin90°/4=\mathbf{+0.25} — a reported change that is spurious: the arrow did not move, the basis did.

Step 3. Apply ∇rVϕ=∂rVϕ+ΓϕrϕVϕ\nabla_r V^\phi=\partial_r V^\phi+\Gamma^\phi{}_{r\phi}V^\phi (only λ=ϕ\lambda=\phi survives). With Γϕrϕ=1/r=0.5\Gamma^\phi{}_{r\phi}=1/r=0.5 and Vϕ=−0.5V^\phi=-0.5, the connection term is (0.5)(−0.5)=−0.25(0.5)(-0.5)=\mathbf{-0.25}, so ∇rVϕ=0.25−0.25=0\nabla_r V^\phi=0.25-0.25=\mathbf{0}; the other components check similarly, giving ∇V=0\nabla V=0 everywhere. The partial wrongly reports +0.25+0.25; the connection subtracts exactly −0.25-0.25, and ∇\nabla correctly certifies the field as unchanging.

3 · Parallel-transport holonomy equals the enclosed area

A vector is parallel-transported around the geodesic triangle N→A→B→N on the unit sphere, where N is the north pole and A, B are two equatorial points 90°90° apart in longitude. By how much does it rotate, and how does that relate to curvature?

Step 1. Each side is a great-circle arc, hence a geodesic, and the three interior angles are each 90°90° (meridians meet the equator at right angles and A, B are a quarter-turn apart). The angular excess is 3π2−π=π2\tfrac{3\pi}{2}-\pi=\mathbf{\tfrac{\pi}{2}}.

Step 2. The holonomy around any closed loop equals the integral of the Gaussian curvature over the enclosed region, Δα=∬K dA\Delta\alpha=\iint K\,\mathrm{d}A. On the unit sphere K=1/R2=1K=1/R^2=1, so Δα\Delta\alpha equals the enclosed area — one octant of 4π4\pi: A=18(4π)=π2≈1.5708A=\tfrac18(4\pi)=\tfrac{\pi}{2}\approx\mathbf{1.5708} sr =90°=\mathbf{90°}. This matches the angular excess of Step 1 — Gauss–Bonnet in miniature.

Step 3. Interpret through the covariant derivative. Transporting around an infinitesimal loop of area δAμν\delta A^{\mu\nu} rotates a vector by δVρ=−RρσμνVσ 12δAμν\delta V^\rho=-R^\rho{}_{\sigma\mu\nu}V^\sigma\,\tfrac12\delta A^{\mu\nu}, and summing these over the triangle reproduces the finite 90°90°; in flat space R=0R=0 and every loop closes. The vector returns rotated by π/2=90°\pi/2=90°, equal to the enclosed area and to ∬K dA\iint K\,\mathrm{d}A — the finite fingerprint of the curvature that [∇μ,∇ν][\nabla_\mu,\nabla_\nu] measures locally.

Physics of the portrait

The covariant derivative is the central operation of tensor calculus, which three people built across half a century: Christoffel supplied the symbols in 1869, Ricci-Curbastro forged covariant differentiation into a systematic "absolute differential calculus," and Levi-Civita gave it its geometric soul through parallel transport.

Origins · 2 minTullio Levi-Civita and how this was discovered
TL
Tullio Levi-Civita
1873–1941 · Italian geometer who created the absolute differential calculus with Ricci and gave the covariant derivative its meaning through parallel transport

Tullio Levi-Civita, born in Padua in 1873, was the student and then collaborator of Gregorio Ricci-Curbastro, and together they turned a set of formal index manipulations into a coherent theory. Ricci had spent the 1880s and 1890s developing "covariant differentiation" — the operation this article is about — as the backbone of what he called the absolute differential calculus, a way of writing geometric relations that hold in every coordinate system. In 1900 the two published the founding memoir of the subject, "Méthodes de calcul différentiel absolu et leurs applications" in Mathematische Annalen, laying out tensors, the covariant derivative, and their transformation laws in essentially the form used today. Its importance was not obvious for fifteen years, until Einstein, groping toward a relativistic theory of gravity, found in the Ricci calculus exactly the language he needed; the 1915 field equations are written in it, and Einstein and Levi-Civita carried on a warm correspondence over signs and conventions. Levi-Civita's own decisive contribution came in 1917 with the notion of parallel transport on an arbitrary manifold. Before it the Christoffel symbols were a computational device; after it the covariant derivative had a picture — it measures how a vector changes relative to being carried "as parallel as possible" along a curve, and curvature is the failure of that transport to close up around a loop. That geometric reading, threaded through this article, is what made the connection intelligible. Levi-Civita moved to Rome in 1918 and became one of the great figures of mathematical physics, but in 1938 the Fascist racial laws stripped him, as a Jew, of his chair and his academy memberships; isolated and in failing health, he died in Rome in 1941. The connection that bears his name outlived the regime that silenced him.

History
1869
— Christoffel introduces the symbols Γ that bear his name while studying quadratic differential forms
1887
— Ricci-Curbastro develops systematic covariant differentiation, the core of the absolute differential calculus
1900
— Ricci and Levi-Civita publish "Méthodes de calcul différentiel absolu", the founding text of tensor calculus
1915
— Einstein writes the field equations of general relativity in the Ricci calculus
1917
— Levi-Civita defines parallel transport, giving the covariant derivative its geometric meaning

Going further

Texts. Misner, Thorne and Wheeler, Gravitation (1973), Ch. 8–10 · Wald, General Relativity (1984), §3.1–3.2 · Carroll, Spacetime and Geometry (2004), Ch. 3 · Schutz, A First Course in General Relativity (2009), Ch. 5–6 · Hartle, Gravity (2003), Ch. 20.

Provenance. Equations and conventions follow Carroll, Spacetime and Geometry, Ch. 3 (covariant derivative 3.1–3.6, Christoffel 3.27, Riemann 3.112) and Wald, General Relativity, §3.1–3.2. A verifier should recompute: (1) the cancellation of the inhomogeneous transformation terms and the forced minus sign for one-forms; (2) the Christoffel formula by symmetrizing the three metric-compatibility equations; (3) the Riemann tensor from the nested covariant derivatives; (4) the two-sphere connection Γθϕϕ=−0.5\Gamma^\theta{}_{\phi\phi}=-0.5, Γϕθϕ=1.0\Gamma^\phi{}_{\theta\phi}=1.0 at θ=45°\theta=45° with ∇θgϕϕ=0\nabla_\theta g_{\phi\phi}=0; (5) the polar-plane check ∂rVϕ=+0.25\partial_r V^\phi=+0.25, ΓϕrϕVϕ=−0.25\Gamma^\phi{}_{r\phi}V^\phi=-0.25, ∇rVϕ=0\nabla_r V^\phi=0 at r=2r=2, ϕ=90°\phi=90°; and (6) the spherical-octant holonomy =π/2=90°=\pi/2=90° equal to the enclosed area. All structural results are standard textbook content; the numbers are elementary and reproducible.

Builds on: The Affine Connection · Christoffel Symbols · Parallel Transport. Leads to: The Geodesic Equation · Metric Compatibility · Geodesic Deviation Equation · The Levi-Civita Connection.

Ask the tutor: Ask me to compute the Christoffel symbols of any metric you name — sphere, Schwarzschild, FLRW — and then covariantly differentiate a field or parallel-transport a vector around a loop, step by numbered step.

Is the covariant derivative just the partial derivative plus a correction?

Yes, and the correction is the whole point. It is literally ∇μ=∂μ+(connection terms)\nabla_\mu=\partial_\mu+(\text{connection terms}), one Γ\Gamma term per index. But do not read "correction" as "small adjustment": neither ∂μVν\partial_\mu V^\nu nor the connection term is a tensor on its own — each depends on the coordinate grid — and only their sum transforms cleanly. The connection is engineered so its non-tensorial transformation piece exactly cancels the non-tensorial piece of ∂\partial. So ∇\nabla is the smallest modification of the partial derivative that yields a genuine tensor, and in flat Cartesian coordinates, where the basis is rigid and every Γ\Gamma vanishes, the correction is zero and ∇\nabla collapses back to ∂\partial.

Why isn't the ordinary partial derivative good enough?

Because in general coordinates it does not produce a tensor, so equations written with it change form under a change of coordinates — they are artifacts of a chosen grid, not physical laws. The failure is geometric: the change of a vector field has two parts, the change of its components and the change of the basis vectors carrying them, and ∂μ\partial_\mu sees only the first. On a curved manifold, or even on a flat plane in polar coordinates, the basis rotates from point to point, so the partial derivative attributes spurious change to a field that is actually constant. The covariant derivative adds the connection precisely to subtract that basis rotation, leaving the true, coordinate-independent rate of change.

Are the Christoffel symbols a tensor?

No, and this is the trap that catches everyone once. Under a coordinate change Γλμν\Gamma^\lambda{}_{\mu\nu} acquires an inhomogeneous term in second derivatives of the transformation — the very term that lets it cancel the non-tensorial part of ∂\partial. Because of it you can make all the Christoffel symbols vanish at any single point (Riemann normal coordinates, the mathematics of a freely falling frame), even in curved spacetime. Were Γ\Gamma a tensor, "zero in one frame" would mean "zero in all frames," and gravity could be transformed away globally — which the equivalence principle forbids beyond a point. The tensorial objects built from Γ\Gamma are its antisymmetric part, the torsion, and the curvature combination ∂Γ+ΓΓ\partial\Gamma+\Gamma\Gamma, the Riemann tensor.

What does metric compatibility actually buy us?

Metric compatibility, ∇σgμν=0\nabla_\sigma g_{\mu\nu}=0, says the metric is covariantly constant, so parallel transport preserves lengths and angles. Practically, it means ∇\nabla commutes with raising and lowering indices: you can slide the metric through a covariant derivative, ∇μVν=gνρ∇μVρ\nabla_\mu V_\nu=g_{\nu\rho}\nabla_\mu V^\rho, with no extra terms. Combined with vanishing torsion it does more still — it fixes the connection uniquely as the Christoffel symbols of the metric. Without metric compatibility there are infinitely many connections; with it and zero torsion there is exactly one, which is why "the" covariant derivative of general relativity is unambiguous.

How is the covariant derivative related to parallel transport?

They are two views of one structure. Parallel transport carries a vector along a curve "without turning it," and a vector is parallel-transported exactly when its covariant derivative along the curve vanishes, uν∇νVμ=0u^\nu\nabla_\nu V^\mu=0. So ∇\nabla measures the rate at which a field departs from being parallel-transported — a nonzero ∇\nabla is the failure to stay parallel. A geodesic is the special curve that parallel-transports its own tangent, which is why uν∇νuμ=0u^\nu\nabla_\nu u^\mu=0 yields the geodesic equation directly. And the curvature that [∇μ,∇ν][\nabla_\mu,\nabla_\nu] detects is precisely the rotation a vector suffers under parallel transport around a closed loop, as the spherical-triangle example makes numerically vivid.

What is the comma-to-semicolon rule and when does it fail?

In old index notation a partial derivative is a comma (Vν,μV^\nu{}_{,\mu}) and a covariant derivative a semicolon (Vν;μV^\nu{}_{;\mu}), so the minimal-coupling recipe "replace every ∂\partial by ∇\nabla" is nicknamed comma-to-semicolon. It promotes a special-relativistic law to curved spacetime: ∂μTμν=0\partial_\mu T^{\mu\nu}=0 becomes ∇μTμν=0\nabla_\mu T^{\mu\nu}=0, and Maxwell's equations become ∇μFμν=μ0Jν\nabla_\mu F^{\mu\nu}=\mu_0 J^\nu. The subtlety is ordering: because covariant derivatives do not commute, a law with two derivatives is ambiguous — different orderings differ by curvature terms, and the rule cannot say which is right. For first-derivative laws it is unambiguous; for higher derivatives it needs physical input about non-minimal curvature couplings.

Before this

Where this leads