SUCURI
U X

tutorial

Differential geometry
through exercises

index notation and index-free notation,
solved — or not — by Sucuri

\nabla_U \nabla_U X = R(U,X)U

85 exercises from ten universities every answer checked by the test suite
SUCURI notebook manual português tutorial

Differential geometry through exercises. This tutorial works through 85 exercises from problem sets used in courses on general relativity and Riemannian geometry — Cambridge, MIT, Caltech, Berlin, Oregon, Nagoya, Utrecht, Leiden —, grouped by the notation they are naturally stated in: index notation, with components, indices and the Einstein summation; and index-free, with vector fields, connections and forms. For each one, what Sucuri does with it: solves it, solves it partly, or does not — and why.

Every exercise solved here is run by the test suite on every change to the program: the printed output is the output. Before each solution comes the full problem statement, quoted from the original with the exact reference and link, and a note on what Sucuri does with it.

How to read

Each exercise gives the source, the problem statement in brief, and the cells that solve it in the notebook, with the answer. Three badges:

solves Sucuri reaches the result of the exercise — by computation, or by prove, which returns the checked chain of steps.
partly it solves one part, or the deciding step, or only with a hypothesis that the exercise asks you to derive.
not solved the exercise asks for something Sucuri cannot yet write or derive — and section V says what.

Two things hold for everything that follows. Declaring decides the reading: an index is only an index, a vector only a vector, once declared — the manual explains why. And prove uses only the named hypotheses: what the engine knows on its own is what holds in any textbook (linearity, Leibniz, antisymmetry of the bracket, d² = 0); the rest, including every sign convention, goes in as a written hypothesis.

The scoreboard

“Solves” means: everything the problem statement asks for, not just the computation at its centre. The count was redone on that basis when the full statements went into the tutorial, ahead of each solution: 23 exercises that came out in the central computation also ask for an embedding, geodesics, a volume, parallel transport, or a proof in the opposite direction — and moved to “partly”. What a symbolic system does not do by nature — draw a graph, sketch, give an interpretation — does not count against it: it is stated in the note on each problem.

solvespartlynot solved
index (30)2460
components (25)2140
index-free (30)2190
total (85)66190

Exercises 1 to 45 came from a first search (Reall, Tong, MIT 2006, Richard, Part III, Cline, Carroll); 46 to 85 from a second, with other sources (Wendl, Ali-Haïmoud, Dray, Evans, Hunt, MIT 2018, Hirata, van Baal, 't Hooft, Hartman). The numbering follows the order of the searches, and each exercise appears in the section for its notation.

Index-free notation went further, and not by chance: that was where Sucuri had a proof engine; now index notation has its own. Two exercises stated in indices (9 and 11) only come out because they were rewritten without indices — which is, after all, what abstract notation exists to allow. In index notation, going from ∇ to Γ and from Γ to ∂g already works (2, 4, 17, 65), and so does the determinant of the metric (5, 64, 66), Ricci as a declared contraction (13, 63), prove with hypotheses written with indices (10, 14, 15, 68, 69), counting independent components (12, 13, 61), coordinates with indices (62), linearized theory (71), and the transformation law of the Christoffel symbols, with both charts on the same index type (3). In components, Sucuri computes Christoffel, Riemann, Ricci and the scalar of any metric — given by its diagonal, by the line element, or induced by an embedding or a change of coordinates —, geodesics (the equations, what is conserved, and the orbits by quadrature), the orthonormal basis with Cartan's forms, and fields by components — ∇A, the Laplacian, ℒXg and the bracket —; and in_chart evaluates any indexed expression component by component (36, 78). What is missing is argument, or finding a transformation.

I. Index notation

1. Symmetric against antisymmetric solves

Problem 1 · Contraction of a symmetrized pair with an antisymmetrized pair

H. S. Reall, Part III General Relativity (2022), §1.7 “Tensors”, exercise right after eq. (1.42), p. 20 — damtp.cam.ac.uk.

Show that T(ab) X[a|cd|b] = 0.

In Sucuri: all of it: the whole identity (section I).

Reall, Part III GR, §1.7: the contraction of a symmetrized pair with an antisymmetrized pair vanishes, T(ab)X[a|cd|b] = 0. The bars exclude c and d from the antisymmetrization.

a, b, c, d = indices T = tensor(2, 0) X = tensor(0, 4) T^{(ab)} X_{[a|cd|b]} simplify(eq1) 0

without the parentheses on T, the same computation does not vanish. This was the first exercise tried, and it found a bug: before the fix, the answer was 2·T·X + 2·T·X (section VI).

15. Conserved current of a Killing field solves

Problem 15 · Conserved current associated with a Killing vector

H. S. Reall, Part III General Relativity (2022), §5.3 “Lie derivative”, last exercise of the section, p. 57 — damtp.cam.ac.uk; and D. Tong, General Relativity: Example Sheet 3 (2019), question 5, p. 2 — davidtong.org.

Reall: Let Ja = TabXb where Tab is symmetric and satisfies ∇aTab = 0 (e.g. the energy-momentum tensor: see later) and Xb is a Killing vector field. Show that ∇aJa = 0, i.e., Ja is a conserved current.

Tong: Let Kμ be a Killing vector field and Tμν the energy momentum tensor. Let Jμ = TμνKν. Show that Jμ is a conserved current, meaning ∇μJμ = 0.

In Sucuri: all of it: both versions: conservation from the two hypotheses.

Reall §5.3 (also Tong, sheet 3, Q5): with Tab symmetric and ∇aTab = 0, and Xb Killing (∇(aXb) = 0), Ja = TabXb is conserved. Both conditions are hypotheses, and prove chains them: the goal must be a combination of them, of ∇ of them, and of products with what appears in the problem.

a, b = indices T = tensor(2, 0, symmetric) X = tensor(0, 1) \nabla = levi-civita g = metric \nabla_a T^{ab} = 0 \nabla_a X_b + \nabla_b X_a = 0 \nabla_a (T^{ab} X_b) = 0 prove(eq3, eq1, eq2) 1/2 · eq2 [a→L_0, b→L_1] × T(-L_0, -L_1)

the certificate: eq1 times X, plus half of eq2 times T — and the sum is checked again before it says proved. Without eq2, it does not come out; without simétrico on T, neither.

16. Contraction of two Levi-Civita tensors partly

Problem 16 · Contraction of two Levi-Civita tensors

H. S. Reall, Part III General Relativity (2022), §8.5 “Volume form”, lemma with eq. (8.54) (“Proof. Optional exercise.”), p. 107 — damtp.cam.ac.uk.

Lemma.

εa₁…ap cp+1…cn εb₁…bp cp+1…cn = ± p!(n − p)! δ[a₁b₁ … δap]bp   (8.54)

where the upper (lower) sign holds for Riemannian (Lorentzian) signature.

Proof. Optional exercise.

[Context: (8.52) ε12…n = √|g| and (8.53) ε12…n = ±1/√|g| in a right-handed chart. The special case εabcdεabcd = −24 (p = 0, n = 4, Lorentzian signature) is not written out in the text.]

In Sucuri, partly: the cases p = 0 and p = 2 in 4D, with both signatures; the lemma for general n and p, no.

Reall §8.5, eq. (8.54): in n dimensions, εa…c…εb…c… = ±p!(n−p)! δ[ab…δ…]…, with the − sign in Lorentzian signature; in particular εabcdεabcd = −24.

\mu, \nu, \rho, \sigma, \alpha, \beta = indices \delta = kronecker g = metric(-,+,+,+) \epsilon = levi-civita(tensor) \epsilon^{\mu\nu\rho\sigma} \epsilon_{\mu\nu\rho\sigma} simplify(eq1) -24
\mu, \nu, \rho, \sigma, \alpha, \beta = indices \delta = kronecker g = metric(-,+,+,+) \epsilon = levi-civita(tensor) \epsilon^{\mu\nu\rho\sigma} \epsilon_{\mu\nu\alpha\beta} simplify(eq1) -2*delta(rho, -alpha)*delta(sigma, -beta) + 2*delta(sigma, -alpha)*delta(rho, -beta)

p = 2: −2! · 2!/2! · (δδ − δδ). The sign comes from the declared signature; with metric(+,+,+,+) the first computation gives +24.

2. ∇ of a (1,1) tensor in Christoffel symbols solves

Problem 2 · Covariant derivative of a (1,1) tensor

H. S. Reall, Part III General Relativity (2022), §3.1 “Introduction”, exercise right after eq. (3.17), p. 34 — damtp.cam.ac.uk.

Now the Leibniz rule can be used to obtain the formula for the coordinate basis components of ∇T where T is a (r, s) tensor:

Tμ₁…μrν₁…νs;ρ = Tμ₁…μrν₁…νs,ρ + Γμ₁σρ Tσμ₂…μrν₁…νs + … + Γμrσρ Tμ₁…μr−1σν₁…νs − Γσν₁ρ Tμ₁…μrσν₂…νs − … − Γσνsρ Tμ₁…μrν₁…νs−1σ   (3.17)

Exercise. Prove this result for a (1, 1) tensor.

[Context: eqs. (3.9) and (3.16) of the same text give, in the coordinate basis, Yμ;ν = Yμ,ν + ΓμρνYρ and ημ;ν = ημ,ν − Γρμνηρ. The result asked for is therefore Tμν;ρ = Tμν,ρ + ΓμσρTσν − ΓσνρTμσ. This is a specialization of (3.17) and does not appear literally in the text.]

In Sucuri: all of it: the (1,1) tensor, in Reall's convention, read from the written definition.

Reall §3.1: starting from the action of ∇ on a vector, write ∇ρTμν in terms of Γ. Reall puts the derivative index last in Γ; Carroll, first. With torsion the difference matters, so the convention is not assumed: it is read from the written definition, and \Gamma = christoffel(eq1) declares it. expand opens each ∇ into ∂ and one Γ per index.

\mu, \nu, \rho, \sigma = indices V = tensor(1, 0) T = tensor(1, 1) \nabla_\rho V^\mu = \partial_\rho V^\mu + \Gamma^\mu{}_{\sigma\rho} V^\sigma \Gamma = christoffel(eq1) \nabla_\rho T^\mu{}_\nu = \partial_\rho T^\mu{}_\nu + \Gamma^\mu{}_{\sigma\rho} T^\sigma{}_\nu - \Gamma^\sigma{}_{\nu\rho} T^\mu{}_\sigma expand(eq2) True

without the −Γ term for the lower index, the answer is not True.

3. The transformation of the Christoffel symbols solves

Problem 3 · The transformation of the Christoffel symbols

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 3, question 4, parts (a)–(b), p. 2 — ocw.mit.edu.

(a) Show that, under a coordinate transformation, the components of the Christoffel symbol transform as follows:

Γα′β′γ′ = (∂xα′/∂xα)(∂xβ/∂xβ′)(∂xγ/∂xγ′) Γαβγ − (∂²xα′/∂xβ∂xγ)(∂xβ/∂xβ′)(∂xγ/∂xγ′).

Do this by considering the form of the Christoffel symbol in terms of derivatives of the metric.

(b) Show that, using this rule, the components of the covariant derivative of a vector transform as tensors should: ∇α′Aβ′ = (∂xα/∂xα′)(∂xβ′/∂xβ) ∇αAβ.

In Sucuri: all of it: (a) and (b), via the form of the Christoffel symbols in terms of derivatives of the metric, as the problem asks. The ∂′ of a product — Leibniz with the chain rule, ∂′λ = Jσλ∂σ — is written as a hypothesis: Sucuri has a single set of coordinates, and the rest is index algebra, which it does.

Both charts on the same index type — both run from 1 to n —, and what links them, named: Jαμ = ∂xα/∂x′μ, its inverse Kμα = ∂x′μ/∂xα, and Hαμν = ∂²xα/∂x′μ∂x′ν, symmetric. Qσαβ = ∂σgαβ; Pλμν = ∂′λg′μν; k = g′−1; G and S, the Christoffel symbols in the two charts. P comes from g′μν = JαμJβνgαβ by Leibniz and the chain rule — it is hypothesis eq3.

\alpha, \beta, \gamma, \kappa, \lambda, \mu, \nu, \rho, \sigma, \tau = indices \delta = kronecker g = metric J = tensor(1, 1) K = tensor(1, 1) H = tensor(1, 2) Q = tensor(0, 3) P = tensor(0, 3) k = tensor(2, 0, symmetric) S = tensor(1, 2) G = tensor(1, 2) H^\alpha{}_{\mu\nu} = H^\alpha{}_{\nu\mu} k^{\mu\nu} = K^\mu{}_\alpha K^\nu{}_\beta g^{\alpha\beta} P_{\lambda\mu\nu} = H^\alpha{}_{\lambda\mu} J^\beta{}_\nu g_{\alpha\beta} + J^\alpha{}_\mu H^\beta{}_{\lambda\nu} g_{\alpha\beta} + J^\sigma{}_\lambda J^\alpha{}_\mu J^\beta{}_\nu Q_{\sigma\alpha\beta} S^\mu{}_{\nu\rho} = \frac{1}{2} k^{\mu\kappa} (P_{\nu\kappa\rho} + P_{\rho\kappa\nu} - P_{\kappa\nu\rho}) G^\tau{}_{\lambda\sigma} = \frac{1}{2} g^{\tau\kappa} (Q_{\lambda\kappa\sigma} + Q_{\sigma\kappa\lambda} - Q_{\kappa\lambda\sigma}) J^\alpha{}_\mu K^\mu{}_\beta = \delta^\alpha_\beta P_{\nu\kappa\rho} + P_{\rho\kappa\nu} - P_{\kappa\nu\rho} = 2 J^\alpha{}_\kappa H^\beta{}_{\nu\rho} g_{\alpha\beta} + J^\lambda{}_\nu J^\alpha{}_\kappa J^\sigma{}_\rho (Q_{\lambda\alpha\sigma} + Q_{\sigma\alpha\lambda} - Q_{\alpha\lambda\sigma}) prove(eq7, eq1, eq3) proved from eq1, eq3

the Christoffel combination: the H terms combine into 2JH, by the symmetry of H.

\alpha, \beta, \gamma, \kappa, \lambda, \mu, \nu, \rho, \sigma, \tau = indices \delta = kronecker g = metric J = tensor(1, 1) K = tensor(1, 1) H = tensor(1, 2) Q = tensor(0, 3) P = tensor(0, 3) k = tensor(2, 0, symmetric) S = tensor(1, 2) G = tensor(1, 2) H^\alpha{}_{\mu\nu} = H^\alpha{}_{\nu\mu} k^{\mu\nu} = K^\mu{}_\alpha K^\nu{}_\beta g^{\alpha\beta} P_{\lambda\mu\nu} = H^\alpha{}_{\lambda\mu} J^\beta{}_\nu g_{\alpha\beta} + J^\alpha{}_\mu H^\beta{}_{\lambda\nu} g_{\alpha\beta} + J^\sigma{}_\lambda J^\alpha{}_\mu J^\beta{}_\nu Q_{\sigma\alpha\beta} S^\mu{}_{\nu\rho} = \frac{1}{2} k^{\mu\kappa} (P_{\nu\kappa\rho} + P_{\rho\kappa\nu} - P_{\kappa\nu\rho}) G^\tau{}_{\lambda\sigma} = \frac{1}{2} g^{\tau\kappa} (Q_{\lambda\kappa\sigma} + Q_{\sigma\kappa\lambda} - Q_{\kappa\lambda\sigma}) J^\alpha{}_\mu K^\mu{}_\beta = \delta^\alpha_\beta k^{\mu\kappa} J^\alpha{}_\kappa g_{\alpha\beta} = K^\mu{}_\beta prove(eq7, eq2, eq6) proved from eq2, eq6

g′−1 contracted with J and g is K.

\alpha, \beta, \gamma, \kappa, \lambda, \mu, \nu, \rho, \sigma, \tau = indices \delta = kronecker g = metric J = tensor(1, 1) K = tensor(1, 1) H = tensor(1, 2) Q = tensor(0, 3) P = tensor(0, 3) k = tensor(2, 0, symmetric) S = tensor(1, 2) G = tensor(1, 2) H^\alpha{}_{\mu\nu} = H^\alpha{}_{\nu\mu} k^{\mu\nu} = K^\mu{}_\alpha K^\nu{}_\beta g^{\alpha\beta} P_{\lambda\mu\nu} = H^\alpha{}_{\lambda\mu} J^\beta{}_\nu g_{\alpha\beta} + J^\alpha{}_\mu H^\beta{}_{\lambda\nu} g_{\alpha\beta} + J^\sigma{}_\lambda J^\alpha{}_\mu J^\beta{}_\nu Q_{\sigma\alpha\beta} S^\mu{}_{\nu\rho} = \frac{1}{2} k^{\mu\kappa} (P_{\nu\kappa\rho} + P_{\rho\kappa\nu} - P_{\kappa\nu\rho}) G^\tau{}_{\lambda\sigma} = \frac{1}{2} g^{\tau\kappa} (Q_{\lambda\kappa\sigma} + Q_{\sigma\kappa\lambda} - Q_{\kappa\lambda\sigma}) J^\alpha{}_\mu K^\mu{}_\beta = \delta^\alpha_\beta P_{\nu\kappa\rho} + P_{\rho\kappa\nu} - P_{\kappa\nu\rho} = 2 J^\alpha{}_\kappa H^\beta{}_{\nu\rho} g_{\alpha\beta} + J^\lambda{}_\nu J^\alpha{}_\kappa J^\sigma{}_\rho (Q_{\lambda\alpha\sigma} + Q_{\sigma\alpha\lambda} - Q_{\alpha\lambda\sigma}) k^{\mu\kappa} J^\alpha{}_\kappa g_{\alpha\beta} = K^\mu{}_\beta S^\mu{}_{\nu\rho} = K^\mu{}_\tau J^\lambda{}_\nu J^\sigma{}_\rho G^\tau{}_{\lambda\sigma} + K^\mu{}_\alpha H^\alpha{}_{\nu\rho} prove(eq9, eq4, eq5, eq7, eq8) proved from eq4, eq5, eq7, eq8

(a) in Carroll's form: Γ′ = K J J Γ + K ∂²x/∂x′∂x′.

\alpha, \beta, \gamma, \mu, \nu, \rho = indices J = tensor(1, 1) K = tensor(1, 1) H = tensor(1, 2) L = tensor(1, 2) M = tensor(1, 2) L^\mu{}_{\nu\alpha} J^\alpha{}_\rho + K^\mu{}_\alpha H^\alpha{}_{\nu\rho} = 0 L^\mu{}_{\nu\alpha} = J^\beta{}_\nu M^\mu{}_{\beta\alpha} K^\mu{}_\alpha H^\alpha{}_{\nu\rho} = - M^\mu{}_{\beta\gamma} J^\beta{}_\nu J^\gamma{}_\rho prove(eq3, eq1, eq2) proved from eq1, eq2

M = ∂²x′/∂x∂x, and L = ∂′K: eq1 is ∂′(KJ) = ∂′δ = 0, eq2 the chain rule. With this, the K H of the form above is the −M J J of the problem statement.

\alpha, \beta, \gamma, \kappa, \lambda, \mu, \nu, \rho, \sigma, \tau = indices \delta = kronecker J = tensor(1, 1) K = tensor(1, 1) H = tensor(1, 2) L = tensor(1, 2) S = tensor(1, 2) G = tensor(1, 2) A = tensor(1, 0) U = tensor(1, 1) T = tensor(1, 1) L^\mu{}_{\nu\alpha} J^\alpha{}_\rho + K^\mu{}_\alpha H^\alpha{}_{\nu\rho} = 0 S^\mu{}_{\nu\rho} = K^\mu{}_\tau J^\lambda{}_\nu J^\sigma{}_\rho G^\tau{}_{\lambda\sigma} + K^\mu{}_\alpha H^\alpha{}_{\nu\rho} T^\mu{}_\nu = L^\mu{}_{\nu\alpha} A^\alpha + K^\mu{}_\alpha J^\beta{}_\nu U^\alpha{}_\beta J^\alpha{}_\mu K^\mu{}_\beta = \delta^\alpha_\beta T^\mu{}_\nu + S^\mu{}_{\nu\rho} K^\rho{}_\gamma A^\gamma = K^\mu{}_\alpha J^\beta{}_\nu (U^\alpha{}_\beta + G^\alpha{}_{\beta\gamma} A^\gamma) prove(eq5, eq1, eq2, eq3, eq4) proved from eq1, eq2, eq3, eq4

(b): U = ∂A, T = ∂′A′ with A′ = K A (Leibniz and chain rule, eq3); ∇′A′ = K J ∇A — the non-tensorial term in Γ′ cancels the one in ∂′A′.

4. The derivatives of the metric solves

Problem 4 · Connection identities: derivatives of the metric

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 3, question 7, parts (a)–(c), p. 3 — ocw.mit.edu.

Prove the following connection identities:

(a) ∂λgμν = Γμνλ + Γνμλ.

(b) gμκ∂λgκν = −gκν∂λgμκ.

(c) ∂λgμν = −Γμλκgκν − Γνλκgκμ.

In Sucuri: all of it: (a), (b) and (c).

MIT 8.962, problem set 3, #7: with the Levi-Civita connection, (a) ∂λgμν = Γμνλ + Γνμλ; (b) gμκ∂λgκν = −gκν∂λgμκ; (c) ∂λgμν = −Γμλκgκν − Γνλκgκμ. expand(eq, g) also writes each Γ in terms of the metric, ½gαδ(∂g + ∂g − ∂g) — which holds only for Levi-Civita, and is only done with it declared.

\mu, \nu, \lambda, \sigma = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \partial_\lambda g_{\mu\nu} = g_{\mu\sigma} \Gamma^\sigma{}_{\nu\lambda} + g_{\nu\sigma} \Gamma^\sigma{}_{\mu\lambda} expand(eq2, g) True
\mu, \nu, \lambda, \kappa = indices g = metric g_{\mu\kappa} \partial_\lambda g^{\kappa\nu} = -g^{\kappa\nu} \partial_\lambda g_{\mu\kappa} simplify(eq1) True

∂λgμν = −gμαgνβ∂λgαβ is not a hypothesis: it is what “inverse” means, and simplify applies it whenever a metric is declared.

\mu, \nu, \lambda, \kappa = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \partial_\lambda g^{\mu\nu} = -\Gamma^\mu{}_{\lambda\kappa} g^{\kappa\nu} - \Gamma^\nu{}_{\lambda\kappa} g^{\kappa\mu} expand(eq2, g) True

17. The Lie derivative in components solves

Problem 17 · Lie derivative of a covector and of the metric

H. S. Reall, Part III General Relativity (2022), §5.3 “Lie derivative”, “Exercises (examples sheet 2)” 1 and 2, eqs. (5.18)–(5.21), pp. 55–56 (the same content is in D. Tong, Example Sheet 1, Q5, and Example Sheet 2, Q7) — damtp.cam.ac.uk.

(1) Derive the formula for the Lie derivative of a covector field ωa in a coordinate basis:

(LXω)μ = Xν∂νωμ + ων∂μXν   (5.18)

Show that this can be written in the basis-independent form

(LXω)a = Xb∇bωa + ωb∇aXb   (5.19)

where ∇ is any torsion-free connection.

(2) Show that the Lie derivative of a (0, 2) tensor field gab in a coordinate basis is

(LXg)μν = Xρ∂ρgμν + gμρ∂νXρ + gρν∂μXρ   (5.20)

and that this can be written in the basis-independent form

(LXg)ab = ∇aXb + ∇bXa   (5.21)

where ∇ is the Levi-Civita connection.

Note that we cannot use abstract indices in (5.18) and (5.20) because they hold only in a coordinate basis.

In Sucuri: all of it: (1) and (2): the forms with ∇ in place of ∂.

Reall §5.3: (ℒXω)μ = Xν∂νωμ + ων∂μXν can be written with ∇ in place of ∂ if the connection is torsion-free; and (ℒXg)μν = ∇μXν + ∇νXμ.

\mu, \nu, \lambda = indices V = tensor(1, 0) X = tensor(1, 0) W = tensor(0, 1) \nabla = levi-civita \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) X^\nu \partial_\nu W_\mu + W_\nu \partial_\mu X^\nu = X^\nu \nabla_\nu W_\mu + W_\nu \nabla_\mu X^\nu expand(eq2) True

without \nabla = levi-civita, Γ is not symmetric, the Γ terms do not cancel, and the answer is not True — which is the exercise's hypothesis, stated.

\mu, \nu, \rho, \lambda = indices V = tensor(1, 0) X = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) X^\rho \partial_\rho g_{\mu\nu} + g_{\mu\rho} \partial_\nu X^\rho + g_{\rho\nu} \partial_\mu X^\rho = g_{\nu\rho} \nabla_\mu X^\rho + g_{\mu\rho} \nabla_\nu X^\rho expand(eq2, g) True

65. Christoffel symbols of conformal metrics partly

Problem 65 · Conformal metrics and Nordstrøm's theory

J. M. Evans, Part II General Relativity, Example Sheet 1 (Cambridge, 2026), question 8, p. 2 — damtp.cam.ac.uk.

Two metrics gαβ and ĝαβ are conformally related if ĝαβ = Ω2gαβ for some scalar function Ω. Show that their Christoffel symbols Γαβγ and Γ̂αβγ are related by

Γ̂αβγ = Γαβγ + Ω−1(δαβΩ,γ + δαγΩ,β − gαδgβγΩ,δ).

In Nordstrøm's theory of gravity the metric is given by gαβ = e2φηαβ, where φ is a scalar function of position and ηαβ is the Minkowski metric. Compute the equation of a geodesic in Nordstrøm's theory and use your result to show that TαTβgαβ is constant on a geodesic, where Tα is the tangent vector corresponding to an affine parameter.

Using the results of question 5, show that a null geodesic is also a null geodesic in Minkowski spacetime, and hence deduce that light rays are not subject to gravitational deflection.

Show that for any time-like geodesic, with suitably chosen parameter μ,

d2xα/dμ2 = −ηαγψ,γ

for some function ψ.

[A comma denotes differentiation with respect to a coordinate, e.g. φ,α = ∂αφ, and the Christoffel symbols are defined by Γαβγ = ½ gαμ(gβμ,γ + gμγ,β − gβγ,μ).]

In Sucuri, partly: the relation between the Christoffel symbols; Nordstrøm's theory, no.

Evans, Part II GR, sheet 1, Q8: if ĝ = Ω²g, then Γ̂αβγ = Γαβγ + Ω−1(δαβ∂γΩ + δαγ∂βΩ − gασgβγ∂σΩ). The left-hand side is the Levi-Civita formula written for ĝ; the right, expanded through the metric g.

\alpha, \beta, \gamma, \sigma, \lambda = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\beta V^\alpha = \partial_\beta V^\alpha + \Gamma^\alpha{}_{\beta\lambda} V^\lambda \Gamma = christoffel(eq1) \frac{1}{2} \Omega^{-2} g^{\alpha\sigma} (\partial_\beta (\Omega^2 g_{\sigma\gamma}) + \partial_\gamma (\Omega^2 g_{\sigma\beta}) - \partial_\sigma (\Omega^2 g_{\beta\gamma})) = \Gamma^\alpha{}_{\beta\gamma} + \Omega^{-1} (g^\alpha{}_\beta \partial_\gamma \Omega + g^\alpha{}_\gamma \partial_\beta \Omega - g^{\alpha\sigma} g_{\beta\gamma} \partial_\sigma \Omega) expand(eq2, g) True

gαβ is the delta: with a declared metric, Sucuri reads it that way.

5, 64 and 66. The determinant of the metric solves (5, 66) partly (64)

Problem 5 · Covariant divergences and the d'Alembertian in a coordinate basis

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 3, question 7, parts (d)–(f), p. 3 — ocw.mit.edu.

(Continuation of question 7: “Prove the following identities of the connection”.) The next three parts rely on an identity I will prove on either Thursday March 2nd or on Tuesday March 7th. The quantity g is the determinant of the metric gμν.

(d) ∇νAμν = |g|−1/2∂ν(|g|1/2Aμν) − ΓλνμAλν in a coordinate basis.

(e) ∇νFμν = |g|−1/2∂ν(|g|1/2Fμν) in a coordinate basis, if Fμν is antisymmetric.

(f) □S ≡ gμν∇μ∇νS = |g|−1/2∂μ(|g|1/2gμν∂νS) in a coordinate basis. (S is a scalar function.)

In Sucuri: all of it: (d), (e) and (f); and the divergence of a vector, which is the special case.

Problem 64 · Variation of the determinant and contraction of the Christoffel symbols

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, 2026), question 10, parts (i)–(iv), p. 3 — damtp.cam.ac.uk.

(i) Let M be an invertible matrix. Show that under a small change δM, the corresponding change in the determinant is, to first order, δ(det M) = (det M) tr(M−1δM). [Hint: if the entries of a matrix A are small then, to first order, det(I + A) = 1 + tr A, where I is the identity matrix.]

(ii) Let gαβ be a metric with Lorentzian signature and let g = det(gαβ). Use the result in (i) to show that

Γαββ = (1/2g) ∂g/∂xα = (1/√−g) ∂(√−g)/∂xα,

where Γαγβ is the Levi-Civita connection. (Note that g < 0 for a metric with Lorentzian signature.)

(iii) A tensor density of weight q is defined to be a quantity that transforms as a tensor under a change in coordinates from {xμ} to {x̃α} but with an additional factor of Δq, where Δ = det(∂xμ/∂x̃α), the Jacobian. Show that g transforms as a scalar density of weight 2.

(iv) For ψ a scalar density of weight q, the covariant derivative is defined by

∇αψ = ∂ψ/∂xα − q Γαββ ψ.

Show that ∇αψ is a covector density of weight q.

In Sucuri, partly: part (ii); (i) is Jacobi's formula, which in Sucuri is the declaration of det g; (iii) and (iv), densities, no.

Problem 66 · Maxwell's equations with √−g

G. 't Hooft, Introduction to General Relativity (Utrecht, 2013), ch. 10 “Electromagnetism”, exercise after eq. (10.12), eqs. (10.13)–(10.14), p. 45 — webspace.science.uu.nl.

[Context: varying Aμ, we find that the inhomogeneous equation becomes DμFμν = gαβDαFβν = −Jν (10.12), and hence receives a contribution from the gravitational field Γλμν and the potential gαβ.]

Exercise: show that Eq. (10.12) can also be written as

∂μ(√−g Fμν) = −√−g Jν,   (10.13)

and that

∂μ(√−g Jμ) = 0.   (10.14)

Thus √−g Jμ is the real conserved current, and Eq. (10.13) implies that √−g acts as the dielectric constant of the vacuum.

In Sucuri: all of it: (10.13), via the identity ∇μFμν = (−g)−1/2∂μ(√−g Fμν), and (10.14): ∂ν of the left-hand side of (10.13) is identically zero.

MIT 8.962, problem set 3, #7(d)–(f): ∇νAμν = |g|−1/2∂ν(|g|1/2Aμν) − ΓλνμAλν, the version with antisymmetric F, and the d'Alembertian; the simplest case is the divergence of a vector. Evans, sheet 3, Q10: Γβαβ = ∂α ln √−g. 't Hooft, (10.13)–(10.14): with antisymmetric F, ∇μFμν = (−g)−1/2∂μ(√−g Fμν). Books write g, with no index, for det gμν; Sucuri reads it that way only once declared — g = det(g) — and without the declaration refuses. From it, ∂g = g gμν∂gμν (Jacobi's formula), and the sign of g comes from the signature.

\mu, \nu, \lambda = indices V = tensor(1, 0) \nabla = levi-civita g = metric g = det(g) \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \nabla_\mu V^\mu = \frac{1}{\sqrt{|g|}} \partial_\mu (\sqrt{|g|} V^\mu) expand(eq2, g) True

without the signature, g is only real and nonzero, and |g| stays |g| — which is enough. With −g in place of √|g|, it does not come out.

\mu, \nu, \lambda = indices V = tensor(1, 0) A = tensor(1, 1) \nabla = levi-civita g = metric g = det(g) \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \nabla_\nu A^\nu{}_\mu = \frac{1}{\sqrt{|g|}} \partial_\nu (\sqrt{|g|} A^\nu{}_\mu) - \Gamma^\lambda{}_{\nu\mu} A^\nu{}_\lambda expand(eq2, g) True

part (d) of the problem set; without the Γ term, it does not come out. Part (f), gμν∇μ∇νS, also comes out.

\alpha, \beta, \lambda = indices V = tensor(1, 0) \nabla = levi-civita g = metric(-,+,+,+) g = det(g) \nabla_\alpha V^\beta = \partial_\alpha V^\beta + \Gamma^\beta{}_{\alpha\lambda} V^\lambda \Gamma = christoffel(eq1) \Gamma^\beta{}_{\alpha\beta} = \partial_\alpha (\ln \sqrt{-g}) expand(eq2, g) True
\mu, \nu, \lambda = indices V = tensor(1, 0) F = tensor(2, 0, antisymmetric) \nabla = levi-civita g = metric(-,+,+,+) g = det(g) \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \nabla_\mu F^{\mu\nu} = (-g)^{-1/2} \partial_\mu (\sqrt{-g} F^{\mu\nu}) expand(eq2, g) True

the term ΓνμλFμλ drops out because Γ is symmetric and F is not: without antissimétrico, it does not come out.

\mu, \nu = indices F = tensor(2, 0, antisymmetric) g = metric g = det(g) \partial_\nu \partial_\mu (\sqrt{-g} F^{\mu\nu}) simplify(eq1) 0

66, (10.14): ∂ν applied to both sides of (10.13) kills the left one — symmetric ∂∂ against antisymmetric F — leaving ∂ν(√−g Jν) = 0.

10. The contracted Bianchi identity solves

Problem 10 · Contracted Bianchi identity

D. Tong, General Relativity: Example Sheet 2 (2019), question 4, p. 1 — davidtong.org.

The Riemann tensor constructed from the Levi-Civita connection obeys the Bianchi identity Rμν[ρσ;λ] = 0. Use this fact to derive the contracted Bianchi identity Gμν;μ = 0, where Gμν = Rμν − ½ R gμν is the Einstein tensor.

In Sucuri: all of it: the contracted identity, from the second Bianchi identity.

Tong, GR, sheet 2, Q4: from the second Bianchi identity, ∇[λRρσ]μν = 0, deduce ∇μGμν = 0. The search contracts the hypothesis twice with the metric — the traces the book takes by hand.

\mu, \nu, \rho, \sigma, \lambda = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) \nabla_\lambda R^\rho{}_{\sigma\mu\nu} + \nabla_\mu R^\rho{}_{\sigma\nu\lambda} + \nabla_\nu R^\rho{}_{\sigma\lambda\mu} = 0 \nabla^\mu (R_{\mu\nu} - \frac{1}{2} g_{\mu\nu} R) = 0 prove(eq4, eq3) contracted

with ∇μRμν = ∇νR — without the ½ — it does not come out.

14 and 68. Second-order Killing solves (14, 68)

Problem 14 · Second derivative of a Killing vector and the Killing vectors of Minkowski

D. Tong, General Relativity: Example Sheet 3 (2019), question 6, p. 2 — davidtong.org.

Show that a Killing vector field Kμ satisfies the equation

∇μ∇νKρ = Rρνμσ Kσ

[Hint: use the identity Rρ[μνσ] = 0.]

Deduce that in Minkowski spacetime the components of Killing covectors are linear functions of the coordinates.

In Sucuri: all of it: the identity, and, in Minkowski, ∂∂K = 0: K linear in the coordinates.

Problem 68 · Second-order equation for Killing covectors; the 10 Killing vectors of Minkowski

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, 2026), question 6, p. 2 — damtp.cam.ac.uk.

Let ξα be a Killing covector field, satisfying ξα;β + ξβ;α = 0 (see question 5 on Example Sheet 2). Use the Ricci identity and Rα[βγδ] = 0 to show that

ξα;βγ = −Rδγαβ ξδ .

In the case of Minkowski space, integrate this equation twice and deduce that there are 10 independent Killing vectors.

In Sucuri: all of it: the identity, in Evans's notation, and the 10 Killing fields of Minkowski.

Tong, sheet 3, Q6: a Killing field satisfies ∇μ∇νKρ = RρνμσKσ; the hint is Rρ[μνσ] = 0. Evans, sheet 3, Q6: the same, in semicolon notation, ξα;βγ = −Rδγαβξδ. The first Bianchi identity is a theorem of vanishing torsion: with \nabla = levi-civita it enters without being a hypothesis, and the proof says when it used it.

\mu, \nu, \rho, \sigma = indices V = tensor(1, 0) K = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) \nabla_\mu K_\nu + \nabla_\nu K_\mu = 0 \nabla_\mu \nabla_\nu K^\rho = R^\rho{}_{\nu\mu\sigma} K^\sigma prove(eq3, eq2) Bianchi (a theorem, from zero torsion)
\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) K = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) \nabla_\mu K_\nu + \nabla_\nu K_\mu = 0 \nabla_\gamma \nabla_\beta K_\alpha = -R^\delta{}_{\gamma\alpha\beta} K_\delta prove(eq3, eq2) Bianchi (a theorem, from zero torsion)

with the sign flipped, it does not come out. The rest of the two exercises — integrating in Minkowski and counting the 10 Killing fields — no: that is integrating and counting.

In Minkowski, R = 0, and the same search proves ∇μ∇νKρ = 0 from Killing's equation. In an inertial chart ∇ = ∂, and ∂μ∂νKρ = 0 means K is linear in the coordinates — what Tong asks for. Integrating twice, as Evans asks, is solving Killing's equation among the linear fields: killing(g, 1) finds them all.

\mu, \nu, \rho, \sigma = indices V = tensor(1, 0) K = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) \nabla_\mu K_\nu + \nabla_\nu K_\mu = 0 R^\rho{}_{\sigma\mu\nu} = 0 \nabla_\mu \nabla_\nu K^\rho = 0 prove(eq4, eq2, eq3) proved from eq2, eq3
x = coordinates(t, x, y, z) g = metric(-1, 1, 1, 1) killing(g, 1) 10

four translations, three rotations and three boosts: the 10 Killing fields of Minkowski. With killing(g, 2), still 10 — the quadratic terms are killed by Killing's equation.

69. |∇φ|² + R constant solves

Problem 69 · Scalar field with ∇∇φ = Ricci

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, 2026), question 4, p. 1 — damtp.cam.ac.uk.

Let φ be a scalar field in curved spacetime such that

∇α∇βφ = Rαβ,

where Rαβ is the Ricci tensor. Show that

∇α(∇β∇βφ) = −2Rαβ∇βφ

and hence deduce that ∇αφ∇αφ + R is constant.

[Hint: use the Ricci identity and the contracted Bianchi identity, ∇βRαβ = ½∇αR.]

In Sucuri: all of it: the conclusion, with the contracted Bianchi identity as a hypothesis, as in the hint.

Evans, sheet 3, Q4: if ∇α∇βφ = Rαβ, then ∇αφ∇αφ + R is constant. The book's proof uses the contracted Bianchi identity — which is exercise 10, and enters here as a hypothesis: a proof above becomes a lemma for the one below.

\alpha, \beta, \gamma, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) \nabla^\mu R_{\mu\nu} = \frac{1}{2} \nabla_\nu R \nabla_\alpha \nabla_\beta \phi = R_{\alpha\beta} \nabla_\gamma (\nabla^\alpha \phi \nabla_\alpha \phi + R) = 0 prove(eq5, eq3, eq4) proved from eq3, eq4

the decisive step is ∇μ∇μ∇νφ = ∇ν∇μ∇μφ + Rνμ∇μφ — the commutator, which the canonical form already knows — and the trace of the hypothesis. With |∇φ|² − R, it does not come out.

6. Second derivatives of a function solves

Problem 6 · Torsion and the commutator of covariant derivatives of a function

D. Tong, General Relativity: Example Sheet 2 (2019), question 2, p. 1 — davidtong.org.

Let ∇ be a connection that is not torsion-free. Let T(X, Y) = ∇XY − ∇YX − [X, Y], where X and Y are vector fields. Show that this defines a (1, 2) tensor field T. This is called the torsion tensor. Show that, for any function f,

2∇[μ∇ν] f = −Tρμν ∇ρ f.

In Sucuri: all of it: the identity, with the torsion written in a chart, Tρμν = Γρμν − Γρνμ; that T is a tensor comes out index-free (problem 21).

Tong, sheet 2, Q2: with torsion, 2∇[μ∇ν]f = −Tρμν∇ρf; without torsion, the derivatives commute. In a chart, the torsion is Tρμν = Γρμν − Γρνμ — and with Γ declared by its definition, without Levi-Civita, expand shows the identity.

\mu, \nu, \rho, \lambda = indices V = tensor(1, 0) \nabla_\mu V^\nu = \partial_\mu V^\nu + \Gamma^\nu{}_{\mu\lambda} V^\lambda \Gamma = christoffel(eq1) \nabla_\mu \nabla_\nu f - \nabla_\nu \nabla_\mu f = -(\Gamma^\rho{}_{\mu\nu} - \Gamma^\rho{}_{\nu\mu}) \nabla_\rho f expand(eq2) True
\mu, \nu = indices \nabla = levi-civita \nabla_\mu \nabla_\nu f - \nabla_\nu \nabla_\mu f simplify(eq1) 0

without \nabla = levi-civita, the same difference does not vanish: Sucuri does not assume zero torsion.

7 and 8. The Ricci identity solves (7, 8)

Problem 7 · Ricci identity

H. S. Reall, Part 3 General Relativity (2022), §4.2 “The Riemann tensor”, exercise with eq. (4.8), p. 42 — damtp.cam.ac.uk.

Context: we saw earlier that, with vanishing torsion, the second covariant derivatives of a function commute. The same is not true of covariant derivatives of tensor fields. The failure to commute arises from the Riemann tensor.

Exercise. Let ∇ be a torsion-free connection. Prove the Ricci identity:

∇c∇dZa − ∇d∇cZa = Rabcd Zb (4.8)

Hint. Show that the equation is true when multiplied by arbitrary vector fields Xc and Yd.

In Sucuri: all of it: the identity (4.8), from the components (4.6) of the Riemann tensor, with the torsion-free connection.

Problem 8 · Ricci identity for a 1-form

D. Tong, General Relativity: Example Sheet 2 (2019), question 3, p. 1 — davidtong.org.

Let ∇ be a torsion-free connection. Derive the analogue of the Ricci identity for a 1-form ω,

2∇[μ∇ν] ωρ = −Rσρμν ωσ.

In Sucuri: all of it: the identity on the 1-form, from the definition of the Riemann tensor on a vector.

Reall §4.2, eq. (4.8), asks to prove the Ricci identity on a vector, ∇c∇dZa − ∇d∇cZa = RabcdZb; Tong, sheet 2, Q3, the analogue on a 1-form. In Sucuri, the identity on the vector is the definition of the Riemann tensor — it is from it that riemann(eq) reads the sign and the slot order, because the books do not agree. From it, the 1-form one comes out:

\mu, \nu, \rho, \sigma = indices V = tensor(1, 0) W = tensor(0, 1) \nabla = levi-civita \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) \nabla_\mu \nabla_\nu W_\rho - \nabla_\nu \nabla_\mu W_\rho simplify(eq2) -R(L_0, -rho, -mu, -nu)*W(-L_0)

2∇[μ∇ν]ωρ = −Rσρμνωσ, Tong's answer. The link between identity (7) and the index-free definition R(X,Y)Z appears in indices, in the manual.

Reall defines the Riemann tensor by R(X, Y)Z and derives its components in a coordinate basis, (4.6). With Γ declared in his convention — the derivative index last — and the connection torsion-free, expand opens up both sides:

\mu, \nu, \rho, \sigma, \tau = indices V = tensor(1, 0) Z = tensor(1, 0) \nabla = levi-civita \nabla_\rho V^\mu = \partial_\rho V^\mu + \Gamma^\mu{}_{\nu\rho} V^\nu \Gamma = christoffel(eq1) \nabla_\rho \nabla_\sigma Z^\mu - \nabla_\sigma \nabla_\rho Z^\mu = (\partial_\rho \Gamma^\mu{}_{\nu\sigma} - \partial_\sigma \Gamma^\mu{}_{\nu\rho} + \Gamma^\tau{}_{\nu\sigma} \Gamma^\mu{}_{\tau\rho} - \Gamma^\tau{}_{\nu\rho} \Gamma^\mu{}_{\tau\sigma}) Z^\nu expand(eq2) True

without \nabla = levi-civita the torsion terms remain; with the sign of the Riemann tensor flipped, it does not come out.

9. The pair-exchange symmetry of the Riemann tensor solves

Problem 9 · Pair-exchange symmetry of the Riemann tensor

S. M. Carroll, Lecture Notes on General Relativity (1997), ch. 3 “Curvature”, paragraph right after eq. (3.81), p. 79 — arxiv.org.

This last property is equivalent to the vanishing of the antisymmetric part of the last three indices:

Rρ[σμν] = 0.   (3.81)

All of these properties have been derived in a special coordinate system, but they are all tensor equations; therefore they will be true in any coordinates. Not all of them are independent; with some effort, you can show that (3.64), (3.78) and (3.81) together imply (3.79). The logical interdependence of the equations is usually less important than the simple fact that they are true.

[References in the text: (3.64) Rρσμν = −Rρσνμ; (3.78) Rρσμν = −Rσρμν; (3.79) Rρσμν = Rμνρσ; (3.80) Rρσμν + Rρμνσ + Rρνσμ = 0.]

In Sucuri: all of it: the pair exchange, rewritten index-free, with the three properties as hypotheses.

Carroll, ch. 3, after eq. (3.81): the antisymmetries in the two pairs and the cyclic identity imply Rρσμν = Rμνρσ. Stated in indices; solved index-free, with the (0,4) tensor evaluated on the vectors and the three properties as general hypotheses:

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) W = tensor(1, 0) Q = tensor(0, 4) \forall A, B, C, D: Q(A,B,C,D) = -Q(B,A,C,D) \forall A, B, C, D: Q(A,B,C,D) = -Q(A,B,D,C) \forall A, B, C, D: Q(A,B,C,D) + Q(A,C,D,B) + Q(A,D,B,C) = 0 Q(X,Y,Z,W) = Q(Z,W,X,Y) prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

Q was declared with no symmetry at all — otherwise the pair exchange would be built in, and there would be nothing to prove. The answer's table lists the instances used: four of the cyclic identity, as in the textbook proof, and the antisymmetries that line them up.

11. The Riemann tensor from its symmetric part solves

Problem 11 · Riemann from its symmetrized part

H. S. Reall, Part III General Relativity (2022), §4.5 “Geodesic deviation”, exercise with eq. (4.25), p. 47 (the symmetries used are eqs. (4.16)–(4.17) of §4.4, p. 45) — damtp.cam.ac.uk.

[…] Hence by measuring the LHS above [of the geodesic deviation equation (4.24)] we can determine Ra(bc)d. From this we can determine Rabcd:

Exercise. Show that, for a torsion-free connection,

Rabcd = ⅔ (Ra(bc)d − Ra(bd)c).   (4.25)

[Symmetries available in the text: (4.16) Rab(cd) = 0 and, for a torsion-free connection, (4.17) Ra[bcd] = 0.]

In Sucuri: all of it: (4.25), rewritten index-free.

Reall §4.5, eq. (4.25): without torsion, Rabcd = ⅔(Ra(bc)d − Ra(bd)c), using Rab(cd) = 0 and Ra[bcd] = 0. Also rewritten index-free:

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) W = tensor(1, 0) Q = tensor(0, 4) \forall A, B, C, D: Q(A,B,C,D) = -Q(A,B,D,C) \forall A, B, C, D: Q(A,B,C,D) + Q(A,C,D,B) + Q(A,D,B,C) = 0 Q(W,X,Y,Z) = \frac{1}{3} \cdot (Q(W,X,Y,Z) + Q(W,Y,X,Z) - Q(W,X,Z,Y) - Q(W,Z,X,Y)) prove(eq3, eq1, eq2) proved from eq1, eq2

⅔ times the difference of the symmetrized parts is ⅓ times the sum of four terms. The \cdot before the parenthesis is not decoration: without it, \frac{1}{3} (…) is the same as writing f(…), and Sucuri would ask.

59. Vector identities in index notation solves

Problem 59 · Vector identities in index notation

R. E. Hunt, Vector Calculus: Example Sheet 1 (Part IA, Cambridge, Lent 2026), question 3, p. 1 — damtp.cam.ac.uk.

Use suffix notation to show that for vector fields u(x) and v(x),

∇ × (u × v) = (∇ · v)u − (∇ · u)v + (v · ∇)u − (u · ∇)v

and

(u · ∇)u = ∇(½|u|²) − u × (∇ × u).

In Sucuri: all of it: both identities.

Hunt, Part IA Vector Calculus, sheet 1, Q3: with ε and δ, show that ∇×(u×v) = (∇·v)u − (∇·u)v + (v·∇)u − (u·∇)v, and that (u·∇)u = ∇(½|u|²) − u×(∇×u). In index notation, with the Euclidean metric and the Levi-Civita tensor — which in Cartesian coordinates is the symbol — and with ∇, which there is ∂:

i, j, k, l, m = indices(3) \delta = kronecker g = metric(euclidiana) \epsilon = levi-civita(tensor) \nabla = levi-civita u = tensor(1, 0) v = tensor(1, 0) \epsilon^{ijk} \nabla_j (\epsilon_{klm} u^l v^m) = u^i \nabla_m v^m - v^i \nabla_l u^l + v^j \nabla_j u^i - u^j \nabla_j v^i simplify(eq1) True

the derivative opens by Leibniz, ∇ε = 0 by the declared connection, and εε becomes δδ − δδ. Drop two terms from the right-hand side and the same computation does not close. This exercise found three bugs (section VI).

i, j, k, l, m = indices(3) \delta = kronecker g = metric(cartesiana) \epsilon = levi-civita(tensor) u = tensor(1, 0) u^j \partial_j u^i = \frac{1}{2} \partial^i (u_j u^j) - \epsilon^{ijk} u_j \epsilon_{klm} \partial^l u^m simplify(eq1) True

with ∂, not ∇, the chart matters: uj = gjkuk sits inside the derivative, and ∂g is zero only in a Cartesian chart. With metric(euclidiana), which states only the signature, a ∂g term is left over — and the answer is not True. Sucuri used to assume ∂g = 0 silently (section VI).

60. Decomposition of a rank-2 tensor partly

Problem 60 · Decomposition of a second-rank tensor

R. E. Hunt, Vector Calculus: Example Sheet 4 (Part IA, Cambridge, Lent 2026), section “Symmetries of Tensors”, question 6, p. 2 — damtp.cam.ac.uk.

(i) By first decomposing into symmetric and antisymmetric parts, show that an arbitrary second rank tensor Tij can be written in the form

Tij = α δij + εijk ωk + Dij

for some scalar α, vector ω and symmetric traceless tensor Dij.

(ii) Let u(x) be a vector field and let Tij = ∂ui/∂xj. Show that α = ⅓ ∇ · u and ω = −½ ∇ × u. Find Dij (the “deviatoric strain tensor”) for the case u(x) = (xy², yz², zx²).

* Verify that (0, 0, 1) is a principal axis of Dij at the point x = (2, 3, 0), and find the others.

In Sucuri, partly: part (i), through the antisymmetric part; (ii), no.

Hunt, Part IA Vector Calculus, sheet 4, Q6: every Tij in 3D is αδij + εijkωk + Dij, with ωk = ½εabkTab and D symmetric and traceless. The part that needs ε is the antisymmetric one: T[ij] = ½εijkεabkTab.

i, j, k, a, b = indices(3) \delta = kronecker \epsilon = levi-civita(symbol) T = tensor(0, 2) T_{ij} - \frac{1}{2} \epsilon_{ijk} \epsilon^{abk} T_{ab} - T_{(ij)} simplify(eq1) 0

51. Curl of the gradient, divergence of the curl solves

Problem 51 · d² = 0 in ℝ³: curl grad = 0 and div curl = 0

C. Wendl, Differential Geometry I, Problem Set 6 (Humboldt-Universität zu Berlin, 2016–17), problem 6(c) — mathematik.hu-berlin.de.

Context (preamble of problem 6 and part (b)): given a volume form μ ∈ Ωn(M) on an n-manifold M, the divergence of X ∈ Vec(M) is the unique real-valued function div(X) : M → ℝ such that LXμ = div(X) μ. On M = ℝ³, with μ = dx ∧ dy ∧ dz and X = Xx∂x + Xy∂y + Xz∂z, div(X) = ∂xXx + ∂yXy + ∂zXz, an expression sometimes also denoted by ∇ · X.

(c) Recall that on ℝ³, the gradient of a function f : ℝ³ → ℝ is the vector field

grad(f) = ∇f := (∂xf)∂x + (∂yf)∂y + (∂zf)∂z,

and the curl of a vector field X = Xx∂x + Xy∂y + Xz∂z is the vector field

curl(X) = ∇ × X := (∂yXz − ∂zXy)∂x + (∂zXx − ∂xXz)∂y + (∂xXy − ∂yXx)∂z.

Using the relations of these operations to differential forms and the exterior derivative, deduce from d² = 0 the formulas

∇ × (∇f) = 0   and   ∇ · (∇ × X) = 0

for all f ∈ C∞(ℝ³) and X ∈ Vec(ℝ³).

In Sucuri: all of it: curl grad = 0 and div curl = 0, with forms and in index notation.

Wendl, Differential Geometry I (HU Berlin), problem set 6, 6(c): both identities come from d² = 0 (section IV). In index notation, they come from ∂ commuting against the antisymmetric ε:

i, j, k = indices(3) \epsilon = levi-civita(symbol) \epsilon^{ijk} \partial_j \partial_k f simplify(eq1) 0
i, j, k = indices(3) \epsilon = levi-civita(symbol) X = tensor(0, 1) \epsilon^{ijk} \partial_i \partial_j X_k simplify(eq1) 0

with the symbol, ∂ε = 0: it is ±1 in every chart.

63. The commutator on a rank-2 tensor solves

Problem 63 · Ricci identity for rank-2 tensors

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, 2026), question 1, parts (i)–(ii), p. 1 — damtp.cam.ac.uk.

The Ricci identity for a vector field Vμ is ∇α∇βVμ − ∇β∇αVμ = RμναβVν.

(i) Deduce the corresponding result for a covector field Wμ by considering ∇α∇β(VμWμ) − ∇β∇α(VμWμ). Is there an easier way to find this result?

(ii) [part corresponding to exercise 63] Given two vector fields Uμ and Vμ, evaluate ∇α∇β(UμVν) − ∇β∇α(UμVν). Deduce that, for an arbitrary tensor field Tμν,

∇α∇βTμν − ∇β∇αTμν = RμσαβTσν + RνσαβTμσ.

Hence show that ∇α∇βTαβ = ∇β∇αTαβ, for any tensor field Tαβ.

In Sucuri: all of it: (i) (problem 8) and (ii).

Evans, Part II GR, sheet 3, Q1(ii): from the Ricci identity on a vector, deduce [∇α,∇β]Tμν, and show that ∇α∇βTαβ = ∇β∇αTαβ. The second part needs the Ricci to be symmetric — a theorem, from the symmetries of the Riemann, which Sucuri uses once Levi-Civita and a metric are declared.

\mu, \nu, \rho, \sigma, \alpha, \beta = indices V = tensor(1, 0) \nabla = levi-civita \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) T = tensor(2, 0) \nabla_\alpha \nabla_\beta T^{\mu\nu} - \nabla_\beta \nabla_\alpha T^{\mu\nu} simplify(eq2) R(mu, L_0, -alpha, -beta)*T(-L_0, nu) + R(nu, L_0, -alpha, -beta)*T(mu, -L_0)
\mu, \nu, \rho, \sigma, \alpha, \beta = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) T = tensor(2, 0) \nabla_\alpha \nabla_\beta T^{\alpha\beta} - \nabla_\beta \nabla_\alpha T^{\alpha\beta} simplify(eq2) 0

without g = metric, two contracted Riemanns are left over: the pair exchange Rρσμν = Rμνρσ and the antisymmetry in ρσ come from ∇g = 0, and without a metric they are not assumed.

13. The Weyl tensor solves

Problem 13 · The Weyl tensor

D. Tong, General Relativity: Example Sheet 2 (2019), question 9, p. 3 — davidtong.org.

In a d-dimensional spacetime, define a tensor

Cμνρσ = Rμνρσ + α(Rμρgνσ + Rνσgμρ − Rμσgνρ − Rνρgμσ) + βR(gμρgνσ − gμσgνρ),

where α and β are constants. Show that Cμνρσ has the same symmetries as Rμνρσ.

What values of α and β give Cμνμσ = 0? Determine them. With this extra condition Cμνρσ is called the Weyl tensor. Show that it vanishes if d = 2, 3.

Setting d = 4, how many independent components do Rμν and Cμνρσ have? Show that in vacuum

∇μCμνρσ = 0.

What does the Weyl tensor represent physically?

In Sucuri: all of it: the symmetries, α and β, the vanishing Weyl in d = 2, 3, the counts in d = 4 and ∇μCμνρσ = 0 in vacuum; the physics question is left to the reader.

Tong, GR, sheet 2, Q9: in d dimensions, Cμνρσ = Rμνρσ + α(Rμρgνσ + Rνσgμρ − Rμσgνρ − Rνρgμσ) + βR(gμρgνσ − gμσgνρ); which α and β give Cμνμσ = 0? Textbooks write the Riemann, the Ricci and the scalar with the same letter, and which pair the Ricci contracts varies: R = ricci(eq2) reads it from the definition. The dimension can be a letter, indices(d).

\mu, \nu, \rho, \sigma = indices(d) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) a, b = constant g^{\mu\rho} (R_{\mu\nu\rho\sigma} + a (R_{\mu\rho} g_{\nu\sigma} + R_{\nu\sigma} g_{\mu\rho} - R_{\mu\sigma} g_{\nu\rho} - R_{\nu\rho} g_{\mu\sigma}) + b R (g_{\mu\rho} g_{\nu\sigma} - g_{\mu\sigma} g_{\nu\rho})) simplify(eq3) (R*a + R*b*d - R*b)*g(-nu, -sigma) + (a*d - 2*a + 1)*Ric(-nu, -sigma)

setting both coefficients to zero gives α = −1/(d−2) and β = 1/((d−1)(d−2)). Ric is R with two indices: internally it gets a different head, because a tensor has only one rank.

\mu, \nu, \rho, \sigma = indices(d) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) g^{\mu\rho} (R_{\mu\nu\rho\sigma} - \frac{1}{d-2} (R_{\mu\rho} g_{\nu\sigma} + R_{\nu\sigma} g_{\mu\rho} - R_{\mu\sigma} g_{\nu\rho} - R_{\nu\rho} g_{\mu\sigma}) + \frac{1}{(d-1)(d-2)} R (g_{\mu\rho} g_{\nu\sigma} - g_{\mu\sigma} g_{\nu\rho})) simplify(eq3) 0

that C has the symmetries of the Riemann also comes out.

That the Weyl vanishes in d = 2, 3, and how many components it has in d = 4, is counting: independent(C, eq…) counts the components of a tensor with the declared symmetries and the given equations. Zero means only the zero tensor has them.

\mu, \nu, \rho, \sigma = indices(3) g = metric C = tensor(0, 4, riemann) C_{\mu\nu\rho\sigma} + C_{\mu\rho\sigma\nu} + C_{\mu\sigma\nu\rho} = 0 C^\mu{}_{\nu\mu\sigma} = 0 independent(C, eq1, eq2) 0

with indices(4), 10; and the Ricci, symmetric, 10 as well.

\mu, \nu, \rho, \sigma, \lambda = indices(4) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) \nabla_\lambda R^\rho{}_{\sigma\mu\nu} + \nabla_\mu R^\rho{}_{\sigma\nu\lambda} + \nabla_\nu R^\rho{}_{\sigma\lambda\mu} = 0 R_{\mu\nu} = 0 \nabla^\mu (R_{\mu\nu\rho\sigma} - \frac{1}{2} (R_{\mu\rho} g_{\nu\sigma} + R_{\nu\sigma} g_{\mu\rho} - R_{\mu\sigma} g_{\nu\rho} - R_{\nu\rho} g_{\mu\sigma}) + \frac{1}{6} R (g_{\mu\rho} g_{\nu\sigma} - g_{\mu\sigma} g_{\nu\rho})) = 0 prove(eq5, eq3, eq4) proved from eq3, eq4

in vacuum, ∇μCμνρσ = 0. Without Rμν = 0, it does not come out.

12. The Riemann in two dimensions solves

Problem 12 · Independent components of the Riemann; the two-dimensional case

D. Tong, General Relativity: Example Sheet 2 (2019), question 8, p. 2 — davidtong.org.

How many independent components does the Riemann tensor (of the Levi-Civita connection) have in two, three and four dimensions? Show that in two dimensions

Rμνρσ = ½ R (gμρgνσ − gμσgνρ).

Discuss the implications for general relativity in two spacetime dimensions.

In Sucuri: all of it: the counts (1, 6 and 20), the 2D formula and, for the discussion, that Gμν ≡ 0 in 2D.

Tong, sheet 2, Q8: how many independent components does the Riemann have in 2, 3 and 4 dimensions? In 2, Rμνρσ = ½R(gμρgνσ − gμσgνρ); what does that say about relativity in 2 dimensions? With the Riemann declared and Levi-Civita, the count uses the symmetries and the first Bianchi identity — 1, 6 and 20.

\mu, \nu, \rho, \sigma = indices(4) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) independent(R) 20

The 2D formula: in_components(eq) checks an identity with the most general tensor the declarations allow — here, the Riemann with its symmetries and Bianchi — and an arbitrary metric, component by component. If it holds for the most general one, it holds for all.

\mu, \nu, \rho, \sigma = indices(2) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) R_{\mu\nu\rho\sigma} = \frac{1}{2} R (g_{\mu\rho} g_{\nu\sigma} - g_{\mu\sigma} g_{\nu\rho}) in_components(eq3) True

without the ½, or in indices(3), it gives False. And the answer to Tong's question: in 2D, Rμν − ½Rgμν = 0 with in_components too — the Einstein tensor vanishes identically, and Einstein's equation says nothing.

61. Components of an antisymmetric tensor partly

Problem 61 · Counting components of antisymmetric tensors

A. Guth, MIT 8.962 General Relativity (Spring 2018), Problem Set 4, Problem 1 (“Counting Tensors”), p. 1 — web.mit.edu.

In lecture we derived the fact that if Sμ1···μR is fully symmetric in its indices, and each index can take on D different values, then the number of independent components is

N = (R + D − 1)! / [R! (D − 1)!].

Suppose that Aμ1···μR is fully antisymmetric in its indices, with each index again taking on D different values. How many independent components does Aμ1···μR have?

In Sucuri, partly: the count dimension by dimension; the general formula in R and D, no.

MIT 8.962 (2018), problem set 4, P1: how many independent components does an antisymmetric tensor have? n(n−1)/2 for two indices, and, for p indices, the binomial of n and p. independent counts, dimension by dimension.

\mu, \nu = indices(4) F = tensor(0, 2, antisymmetric) independent(F) 6

with tensor(0, 3, antisymmetric), 4; with simétrico, 10.

62. Killing fields of Euclidean space solves

Problem 62 · Rotational Killing fields in 3D Euclidean space

Y. Ali-Haïmoud, General Relativity (NYU, Fall 2019), Homework 6, Exercise 1, parts (v) and (vi) (with the preamble that precedes them), p. 1 — cosmo.nyu.edu.

Consider flat 3-D space, with cartesian coordinates x1, x2, x3 and associated coordinate basis {∂(i)}, i = 1, 2, 3, and line element dℓ² = δijdxidxj. Define the three vector fields V(i) ≡ εijkxj∂(k), where εijk is the Levi-Civita symbol (fully antisymmetric, such that ε123 = 1), and we sum over repeated indices, even if they are not up and down.

(v) Show that these three vector fields are Killing vector fields.

(vi) Show that they satisfy the commutation relations [V(i), V(j)] = −εijkV(k).

In Sucuri: all of it: (v) and (vi).

NYU, homework 6, 1(v, vi): in ℝ³ with dℓ² = δijdxidxj, the fields V(i) = εijkxj∂k are Killing, and [V(i), V(j)] = −εijkV(k). The coordinates come in with an index: x = coordinates says that ∂jxi = δij. The chart is Cartesian, and that is a declaration: metric(cartesiana).

i, j, k, l, m, n = indices(3) \delta = kronecker g = metric(cartesiana) \epsilon = levi-civita(tensor) x = coordinates \partial_k (\epsilon_{ijl} x^j) + \partial_l (\epsilon_{ijk} x^j) simplify(eq1) 0
i, j, k, l, m, n = indices(3) \delta = kronecker g = metric(cartesiana) \epsilon = levi-civita(tensor) x = coordinates \epsilon_{i m l} x^m \partial^l (\epsilon_{j n k} x^n) - \epsilon_{j m l} x^m \partial^l (\epsilon_{i n k} x^n) = -\epsilon_{i j m} \epsilon^m{}_{n k} x^n simplify(eq1) True

the commutator [V,W]k = Vl∂lWk − Wl∂lVk. With the sign flipped, it does not come out; with metric(euclidiana), which does not state the chart, ∂ε is left over.

71. The linearized Riemann solves

Problem 71 · Linearized curvature

J. M. Evans, Part II General Relativity, Example Sheet 4 (Cambridge, 2026), question 10, part (i), p. 4 — damtp.cam.ac.uk.

(i) A weak gravitational field has the spacetime metric gαβ = ηαβ + εhαβ + O(ε2), where ηαβ is the Minkowski metric and ε is small. Show that

Rαβγδ = ½ε(hαδ,βγ + hβγ,αδ − hαγ,βδ − hβδ,αγ) + O(ε2).

Let h = hγγ and define h̄αβ = hαβ − ½h ηαβ. Check that hαβ = h̄αβ − ½h̄ ηαβ where h̄ = h̄γγ, and show that

Rαβ = ½ε(−□h̄αβ + h̄αγ,βγ + h̄βγ,αγ + ½ηαβ□h̄) + O(ε2),

where □ = ηαβ∂α∂β. What is the linearised vacuum Einstein equation for h̄αβ?

In Sucuri: all of it: the linearized Riemann and Ricci; the vacuum equation is read off the result.

Evans, sheet 4, Q10(i): with gαβ = ηαβ + εhαβ + O(ε²), Rαβγδ = ½ε(hαδ,βγ + hβγ,αδ − hαγ,βδ − hβδ,αγ) + O(ε²), and the Ricci via the trace-reversed h. linearize(eq, h) opens ∇ into Γ and Γ into ∂g, replaces g by η + εh — the inverse by η − εh — and truncates at order ε. The η keeps the metric's name, and it is what raises and lowers the indices of h.

\alpha, \beta, \gamma, \delta, \rho, \lambda = indices V = tensor(1, 0) h = tensor(0, 2, symmetric) \nabla = levi-civita g = metric \nabla_\alpha V^\beta = \partial_\alpha V^\beta + \Gamma^\beta{}_{\alpha\lambda} V^\lambda \Gamma = christoffel(eq1) g_{\alpha\rho} (\nabla_\gamma \nabla_\delta V^\rho - \nabla_\delta \nabla_\gamma V^\rho) = \frac{\epsilon}{2} (\partial_\beta \partial_\gamma h_{\alpha\delta} + \partial_\alpha \partial_\delta h_{\beta\gamma} - \partial_\beta \partial_\delta h_{\alpha\gamma} - \partial_\alpha \partial_\gamma h_{\beta\delta}) V^\beta linearize(eq2, h) True

the left-hand side is RαβγδVβ, by the definition of the Riemann in Evans's convention. With a sign flipped, it does not come out. The Ricci with h̄ comes out too, and is in the tests.

67 and 70. The exterior derivative in components solves (67) partly (70)

Problem 67 · Bianchi identity for F = dA and conservation of the Maxwell stress tensor

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, 2026), question 5, p. 2 — damtp.cam.ac.uk.

The Maxwell tensor Fαβ for the electromagnetic field in curved spacetime is given in terms of a vector potential Aα by Fαβ = ∇αAβ − ∇βAα. Show that this implies ∇[γFαβ] = 0. Show further that if ∇βFαβ = 0, then the energy momentum tensor

Tαβ = FαγFβγ − ¼ gαβFγδFγδ,

is conserved, i.e. ∇βTαβ = 0.

In Sucuri: all of it: ∇[γFαβ] = 0 for F = dA written with ∇, via the first Bianchi identity, and the conservation of Tαβ.

Problem 70 · Exterior derivative in coordinates: the case of 1-forms and 2-forms

C. Wendl, Differential Geometry I (Humboldt-Universität zu Berlin, 2016–17), Problem Set 6, Problem 4(b) (pp. 2–3) — mathematik.hu-berlin.de.

4. Recall that if (x1, …, xn) : U → ℝn is a chart defined on an open subset in some n-manifold M, any k-form ω ∈ Ωk(M) can be written on U as

ω = ωi1…ik dxi1 ⊗ … ⊗ dxik = (1/k!) ωi1…ik dxi1 ∧ … ∧ dxik = Σi1<…<ik ωi1…ik dxi1 ∧ … ∧ dxik,

where the first two expressions use the Einstein summation convention and the third one does not. Here the component functions ωi1…ik : U → ℝ can be written in terms of the coordinate vector fields ∂1, …, ∂n as ωi1…ik = ω(∂i1, …, ∂ik). In order to write down a coordinate formula for the exterior derivative, we introduce the following notation: given any collection of functions Ti1…ik on U labeled by the indices i1, …, ik, define

T[i1…ik] := (1/k!) Σσ∈Sk (−1)|σ| Tiσ(1)…iσ(k),

so for instance if Ti1…ik are the components of a tensor field T, then Alt(T)i1,…,ik = T[i1…ik], and the wedge product of α ∈ Ωk(M) and β ∈ Ωℓ(M) can now be written in coordinates as

(α ∧ β)i1…ikj1…jℓ = α[i1…ikβj1…jℓ].

(a) Prove that the exterior derivative d : Ωk(M) → Ωk+1(M) satisfies (dω)i1…ik+1 = (k + 1) ∂[i1ωi2…ik+1].

(b) Show that for λ ∈ Ω1(M) and ω ∈ Ω2(M), the above formula reduces to

(dλ)ij = ∂iλj − ∂jλi,  and  (dω)ijk = ∂iωjk + ∂jωki + ∂kωij.

In Sucuri, partly: the cyclic sum, by explicit antisymmetrization; ∂[iωjk], with the antisymmetrization reaching across the derivative's index, cannot be written yet.

Evans, sheet 3, Q5: with F = dA, the Bianchi identity ∂[γFαβ] = 0 (the conservation of the Maxwell stress tensor, the second part, needs hypotheses in index notation). Wendl, problem set 6, 4(b): for a 2-form, (dω)ijk is the cyclic sum — the explicit antisymmetrization over the six permutations, with ω antisymmetric, gives the three terms.

\alpha, \beta, \gamma = indices A = tensor(0, 1) \partial_\gamma (\partial_\alpha A_\beta - \partial_\beta A_\alpha) + \partial_\alpha (\partial_\beta A_\gamma - \partial_\gamma A_\beta) + \partial_\beta (\partial_\gamma A_\alpha - \partial_\alpha A_\gamma) simplify(eq1) 0
i, j, k = indices \omega = tensor(0, 2, antisymmetric) \frac{1}{2} (\partial_i \omega_{jk} - \partial_i \omega_{kj} + \partial_j \omega_{ki} - \partial_j \omega_{ik} + \partial_k \omega_{ij} - \partial_k \omega_{ji}) - (\partial_i \omega_{jk} + \partial_j \omega_{ki} + \partial_k \omega_{ij}) simplify(eq1) 0

antisymmetrization that reaches across the derivative's index, ∂[iωjk], cannot be written yet: brackets work only within a single factor.

With ∇, as Evans writes Fαβ = ∇αAβ − ∇βAα, the cyclic sum goes through commutators — and what is left is the cyclic sum of the Riemann, which only the first Bianchi identity closes. prove uses it as a theorem; and the conservation of the energy-momentum tensor follows from the two hypotheses.

\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) A = tensor(0, 1) \nabla_\gamma (\nabla_\alpha A_\beta - \nabla_\beta A_\alpha) + \nabla_\alpha (\nabla_\beta A_\gamma - \nabla_\gamma A_\beta) + \nabla_\beta (\nabla_\gamma A_\alpha - \nabla_\alpha A_\gamma) = 0 prove(eq2) Bianchi (a theorem, from zero torsion)
\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) F = tensor(0, 2, antisymmetric) \nabla_\beta F^{\alpha\beta} = 0 \nabla_\gamma F_{\alpha\beta} + \nabla_\alpha F_{\beta\gamma} + \nabla_\beta F_{\gamma\alpha} = 0 \nabla_\beta (F^\alpha{}_\gamma F^{\beta\gamma} - \frac{1}{4} g^{\alpha\beta} F_{\gamma\delta} F^{\gamma\delta}) = 0 prove(eq4, eq2, eq3) proved from eq2, eq3

without the Bianchi identity, eq3, it does not come out.

II. Components in a chart

Here Sucuri is at home: with the coordinates and the diagonal metric declared, christoffel, riemann, ricci and scalar do the whole calculation.

35 and 36. The sphere solves (35) partly (36)

Problem 35 · Unit sphere: Christoffel, Ricci and Ricci scalar

J. M. Cline, A Short Course in General Relativity (2026), ch. 7, §7.5 Problems, problem 7.2, p. 43 — arxiv.org.

(a) Compute the nonvanishing Christoffel symbols for the unit sphere metric, with ds² = dθ² + sin²θ dϕ².

(b) Compute the Ricci tensor.

(c) Compute the Ricci scalar.

In Sucuri: all of it: (a), (b) and (c).

Problem 36 · Curvature of a sphere

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 5, question 4 (“Curvature of a sphere”), parts (a)–(c), p. 3 — ocw.mit.edu.

(a) Compute all the nonvanishing components of the Riemann tensor Rijkl [(i, j, k, l) ∈ (θ, φ)] for the surface of a 2-sphere.

(b) Consider the parallel transport of a tangent vector A = Aθ eθ + Aφ eφ on the sphere around an infinitesimal parallelogram of sides eθ dθ and eφ dφ. Using the results of part (a), show that to first order in dΩ ≡ sin θ dθ dφ, the length of A is unchanged, but its direction rotates through an angle equal to dΩ.

(c) Show that, if A is parallel transported around the boundary of any simply connected solid angle Ω, its direction rotates through an angle Ω. (“Simply connected” is a topological term meaning that the boundary of the region could be shrunk to a point; it tells us that there are no holes in the manifold or other pathologies.) Compare with the result of Problem 3 (parallel transport around the parallel θ = θ0 on the 2-sphere of radius a, ds² = a²(dθ² + sin²θ dφ²)).

In Sucuri, partly: (a) the Riemann tensor, and (b): δA = R(·)A dθ dφ is orthogonal to A and has |δA| = sin θ |A| dθ dφ — a rotation by dΩ; (c), for an arbitrary region, is summing the infinitesimal rotations, an argument, and does not count.

Cline, problem 7.2 (unit sphere: Christoffel, Ricci, R = 2); MIT 8.962, PS5 #4, with Carroll's worked example, eq. (3.105) (radius a: R = 2/a²).

x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) ricci(g) R_{{\phi}{\phi}}=sin(theta)**2

and Rθθ = 1, Γθφφ = −sin 2θ/2 = −sin θ cos θ, Γφθφ = cot θ.

x = coordinates(\theta, \phi) g = metric(a^2, a^2 \sin^2\theta) scalar(g) 2/a**2
a, b, c, d, e, f, h, k = indices(2) V = tensor(1, 0) \nabla = levi-civita x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) \nabla_a \nabla_b V^c - \nabla_b \nabla_a V^c = R^c{}_{dab} V^d R = riemann(eq1) A = field() R^a{}_{bcd} A^b in_chart(eq2) ^{\theta}_{\theta}_{\phi}: A^phi(theta, phi)*sin(theta)**2

36(b): around the parallelogram, δAa = RabθφAbdθdφ: δAθ = sin²θ Aφ, δAφ = −Aθ.

a, b, c, d, e, f, h, k = indices(2) V = tensor(1, 0) \nabla = levi-civita x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) \nabla_a \nabla_b V^c - \nabla_b \nabla_a V^c = R^c{}_{dab} V^d R = riemann(eq1) A = field() g_{ab} A^a R^b{}_{cef} A^c in_chart(eq2) 0

δA ⟂ A: the length does not change to first order.

a, b, c, d, e, f, h, k = indices(2) V = tensor(1, 0) \nabla = levi-civita x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) \nabla_a \nabla_b V^c - \nabla_b \nabla_a V^c = R^c{}_{dab} V^d R = riemann(eq1) A = field() g_{ab} R^a{}_{cef} A^c R^b{}_{dhk} A^d - (g_{eh} g_{fk} - g_{ek} g_{fh}) g_{cd} A^c A^d in_chart(eq2) 0

in the θφθφ component, |δA|² = sin²θ|A|² dθ²dφ²: the rotation is sin θ dθ dφ = dΩ.

37. The plane in polar coordinates solves

Problem 37 · Covariant derivative and vector Laplacian in polar coordinates

J. M. Cline, A Short Course in General Relativity (2026), Problem 7.6, §7.5 “Problems”, pp. 43–44 — arxiv.org.

7.6. (a) Compute the components Ar;r, Ar;θ, Aθ;r, Aθ;θ of a vector in 2D polar coordinates. Recall that Γrθθ = −r and Γθrθ = 1/r.

(b) Using the results of (a), find Ar;rr, Ar;θθ, Aθ;rr and Aθ;θθ.

(c) Put these together to obtain the covariant Laplacian gij∇i∇jAk for the covariant components of A⃗. It is known as the vector Laplacian. Hint: you can consider Aθ to be dimensionless, and Ar to have dimensions of 1/r, as a way of checking dimensions. I get

∇2Ar = Ar,rr + r−2(Ar,θθ − 2r−1Aθ,θ + rAr,r − Ar),
∇2Aθ = Aθ,rr + r−2Aθ,θθ − r−1Aθ,r + 2r−1Ar,θ.   (7.44)

(This does not agree with Wolfram MathWorld, because the latter assumes a basis where the θ unit vector is normalized to 1, whereas x̂θ · x̂θ = gθθ = 1/r2 here. In other words, Aθ is divided by r in Wolfram to give it the same dimensions as Ar. See the box on coordinate versus orthonormal bases below.)

In Sucuri: all of it: (a) and (b), the first and second covariant derivatives, and (c) the vector Laplacian, which agrees with (7.44).

Fields by components in a chart: A = field() is the generic vector, Ar(r, θ)∂r + Aθ(r, θ)∂θ, and covector() the covector. nabla(g, A) gives Aj;i and Aj;ik; laplacian(g, A), gik∇i∇kA.

x = coordinates(r, \theta) g = metric(1, r^2) A = field() nabla(g, A) -r*A^theta(r, theta) + Derivative(A^r(r, theta), theta)

(a): Ar;θ = ∂θAr − rAθ, from Γrθθ = −r; the table lists the four first derivatives and the eight second ones, (b).

x = coordinates(r, \theta) g = metric(1, r^2) A = covector() laplacian(g, A) Derivative(A_theta(r, theta), (r, 2)) + 2*Derivative(A_r(r, theta), theta)/r - Derivative(A_theta(r, theta), r)/r + Derivative(A_theta(r, theta), (theta, 2))/r**2

(c): the covariant Laplacian of the covariant components, gij∇i∇jAk — Cline's (7.44) term by term, both components.

38. Rindler solves

Problem 38 · Connection in Rindler spacetime

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 4, question 1, p. 1 — ocw.mit.edu.

The spacetime for an accelerated observer that we derived on Pset 2,

ds² = −(1 + g x̄)² dt̄² + dx̄² + dȳ² + dz̄²   (1)

is known as “Rindler spacetime”. Compute all non-zero Christoffel symbols for this spacetime. (Carroll problem 3.3 will help you quite a bit here.)

In Sucuri: all of it: the Christoffel symbols.

MIT 8.962, PS4 #1: the Christoffel symbols of ds² = −(1 + g x)²dt² + dx² + dy² + dz². The acceleration is called k below, because g is the metric.

x = coordinates(t, x, y, z) g = metric(-(1 + k x)^2, 1, 1, 1) christoffel(g) \Gamma^{x}_{{t}{t}}=k*(k*x + 1)

and Γttx = k/(1 + kx). And riemann(g) answers all components vanish: it is Minkowski in another chart.

39. Static spacetime in 1+1 partly

Problem 39 · Riemann tensor for 1+1 static spacetimes

MIT OpenCourseWare, 8.962 General Relativity (Spring 2020), Problem Set 5, question 6, parts (a)–(b), p. 3 — ocw.mit.edu.

(a) Compute all the nonvanishing components of the Riemann tensor for the spacetime with line element

ds² = −e2φ(x) dt² + e−2ψ(x) dx².

(b) For the case φ = ψ = ½ ln|g(x − x0)| where g and x0 are constants, show that the spacetime is flat and find a coordinate transformation to globally flat coordinates (t̄, x̄) such that ds² = −dt̄² + dx̄².

In Sucuri, partly: part (a); in (b), the change to flat coordinates, no.

MIT 8.962, PS5 #6: for ds² = −e2φ(x)dt² + e−2ψ(x)dx², the Riemann tensor and the curvature scalar.

x = coordinates(t, x) e = euler \phi = \phi(x) \psi = \psi(x) g = metric(-e^{2\phi}, e^{-2\psi}) scalar(g) -2*(Derivative(phi(x), x)**2 + Derivative(phi(x), x)*Derivative(psi(x), x) + Derivative(phi(x), (x, 2)))*exp(2*psi(x))

R = −2e2ψ(φ″ + φ′² + φ′ψ′). With φ = ψ = ½ ln|g(x−x₀)| the bracket vanishes and the space is flat, as part (b) asks. This exercise found two bugs (section VI): the declaration \phi = \phi(x) was not accepted, and the metric came out flat.

40 and 41. Schwarzschild and the f(r) family solves (40, 41)

Problem 40 · Christoffel symbols and curvature of a static spherically symmetric metric

D. Tong, General Relativity: Example Sheet 2 (2019), question 11*, p. 3 — davidtong.org.

Consider metrics of the form

ds2 = −f(r)2 dt2 + f(r)−2 dr2 + r2(dθ2 + sin2θ dφ2).

Use the action for a test particle to write down the geodesic equations in this metric, and hence extract the Christoffel symbols in coordinates (t, r, θ, φ).

Use a basis of vierbeins to determine the curvature 2-form, and hence the components of the Riemann tensor in coordinates (t, r, θ, φ).

In Sucuri: all of it: the geodesic equations — from the Christoffel symbols and checked against the test-particle action —, the Christoffel symbols, the tetrad, the curvature 2-forms and the Riemann tensor, in the frame and in coordinates.

Problem 41 · Schwarzschild curvature 2-forms in an orthonormal basis

H. S. Reall, Part 3 General Relativity (Cambridge, 2022), §8.4 “Curvature 2-forms”, Exercise (examples sheet 4) after the Example, eqs. (8.45)–(8.51), pp. 105–106 — damtp.cam.ac.uk.

Context (§8.2, eq. (8.27), p. 102): the Schwarzschild spacetime admits the obvious tetrad

e0 = f dt,   e1 = f−1 dr,   e2 = r dθ,   e3 = r sin θ dφ,   where f = √(1 − 2M/r),

with connection 1-forms ω01 = f′e0, ω21 = (f/r)e2, ω31 = (f/r)e3, ω32 = (1/r) cot θ e3 (and those related by ωμν = −ωνμ).

Context (§8.4): in an orthonormal basis the curvature 2-forms are Θμν = ½ Rμνρσ eρ ∧ eσ (8.45); the antisymmetry of the Riemann tensor implies Θμν = −Θνμ (8.46); and Θμν = dωμν + ωμρ ∧ ωρν (8.47). In the Example, from ω01 = f′e0 one obtains

Θ01 = −Θ01 = (ff″ + f′2) e0 ∧ e1  (8.50),   R0101 = −R0110 = ff″ + f′2 = ½(f2)″ = −2M/r3  (8.51).

Exercise (examples sheet 4). Determine the remaining curvature 2-forms Θ02, Θ03, Θ12, Θ13, Θ23 (all others are related to these by (8.46)). Hence determine the Riemann tensor components. Check that the Ricci tensor vanishes.

In Sucuri: all of it: the five requested 2-forms, the frame components of the Riemann tensor, and the vanishing Ricci tensor.

The orthonormal basis: cartan(g) takes the tetrad ea = √|gaa| dxa of a diagonal metric, finds the connection forms from dea = −ωab∧eb with ωab = −ωba — both checked —, and gives Θab = dωab + ωac∧ωcb, in the frame and in coordinates, the Riemann tensor Rabcd and the Ricci tensor in the frame. geodesics(g) checks its equations against the Euler–Lagrange equations of the action ∫ g(ẋ, ẋ) dλ, which is Tong's route.

x = coordinates(t, r, \theta, \varphi) f = f(r) g = metric(-f^2, \frac{1}{f^2}, r^2, r^2 \sin^2\theta) christoffel(g) \Gamma^{r}_{{t}{t}}=f(r)**3*Derivative(f(r), r)

and Γttr = f′/f, Γrrr = −f′/f, Γrθθ = −rf², all as on the sheet; geodesics(g) gives the equations the sheet extracts them from.

x = coordinates(t, r, \theta, \phi) f = f(r) g = metric(-f^2, f^{-2}, r^2, r^2 \sin^2\theta) cartan(g) (1 - f(r)**2)/r**2*e^2∧e^3

Tong's tetrad, e0 = f dt, e1 = dr/f, …, and the six 2-forms; Θ01 = −(ff″ + f′²) e0∧e1, which is Reall's (8.50).

x = coordinates(t, r, \theta, \phi) g = metric(-(1 - \frac{2m}{r}), \frac{1}{1 - \frac{2m}{r}}, r^2, r^2 \sin^2\theta) cartan(g) Ricci na base = todas as componentes nulas

41: with f = √(1 − 2m/r), Θ02 = Θ03 = −(m/r³)e0∧e2,3, Θ12 = Θ13 = −(m/r³)e1∧e2,3, Θ23 = (2m/r³)e2∧e3; R0101 = −2m/r³, R2323 = 2m/r³; and the Ricci tensor vanishes.

42. The hyperbolic plane solves

Problem 42 · Christoffel symbols and geodesics of the hyperbolic plane

J. Ross, Part III 2016: Differential Geometry — Example Sheets (Cambridge, 2016), Example Sheet 4, Question 9 (p. 13 of the compiled PDF) — dec41.user.srcf.net.

9. Define hyperbolic space H to be the upper half plane {(x, y) : y > 0} and set

g = (dx ⊗ dx + dy ⊗ dy) / y2.

Verify that g defines a Riemannian metric on H, and compute the Christoffel symbols of the Levi-Civita connecion [sic]. Show that γ = (u, v) is a geodesic if and only if

u″v = 2u′v′  and  v″v = v′2 − u′2

and use this to find the geodesics.

In Sucuri: all of it: g is Riemannian (the diagonal is positive for y > 0), the Christoffel symbols, the geodesic equations and the geodesics: semicircles (x − x₀)² + y² = κ/L² and, for L = 0, vertical lines.

The Christoffel symbols and the curvature come from the metric; the geodesics, from geodesics(g) — ẍa + Γabcẋbẋc = 0, with the conserved quantities — and from orbits(g), which integrates them by quadrature when a coordinate is cyclic.

x = coordinates(x, y) g = metric(\frac{1}{y^2}, \frac{1}{y^2}) scalar(g) -2

R = 2K, so Gaussian curvature K = −1. And Γxxy = −1/y, Γyxx = 1/y, Γyyy = −1/y.

x = coordinates(x, y) g = metric(\frac{1}{y^2}, \frac{1}{y^2}) geodesics(g) Derivative(x(lambda), (lambda, 2)) - 2*Derivative(x(lambda), lambda)*Derivative(y(lambda), lambda)/y(lambda)

times y, the two equations are u″v = 2u′v′ and v″v = v′² − u′², with γ = (u, v): those of the problem statement.

x = coordinates(x, y) g = metric(\frac{1}{y^2}, \frac{1}{y^2}) orbits(g) Eq(x - x_0, -sqrt(-L**2*y**2 + kappa)/L)

with L = gxxẋ and κ = g(ẋ, ẋ): (x − x₀)² + y² = κ/L², semicircles centred on the x axis. The case L = 0, which the quadrature excludes, is ẋ = 0: the vertical lines.

43. The 3-sphere solves

Problem 43 · The unit 3-sphere

J. M. Cline, A Short Course in General Relativity (2026), ch. 7, §7.5 Problems, problem 7.10, pp. 44–45 — arxiv.org.

The unit 3-sphere can be embedded into 4D Euclidean space with coordinates (w, x, y, z) by the constraint x² + y² + z² + w² = 1.

(a) Show that the constraint is satisfied by parametrizing the Euclidean coordinates by

w = cψ,   x = sψsθcϕ,   y = sψsθsϕ,   z = sψcθ   (7.46)

where c and s stand for cosine and sine of the subscripted angles 1 ≤ θ, ψ ≤ π and 0 ≤ ϕ ≤ 2π. [sic: read 0 ≤ θ, ψ ≤ π]

(b) By computing the line element in Euclidean space, show that the induced metric is given by diag(1, sψ2, sψ2sθ2) for the coordinates (ψ, θ, ϕ). Hint: start with dx² + dy² to isolate the dϕ² term, then add dz² to get the dθ² term. Add dw² last. Find the volume element and compute the volume of the 3-sphere.

(c) Plug the metric into a symbolic manipulation package and show that the Ricci curvature is R = 6. Also find Rij.

(d) Compute the volume V of a small region ψ < ϵ to O(ϵ5) and compare it to the corresponding volume VE of a sphere of radius ϵ in flat (Euclidean) space. Thereby check Eq. (7.20), V/VE = 1 − r²R/[6(d + 2)] + O(r4), for the 3-sphere.

In Sucuri: all of it: (a) the constraint, (b) the induced metric, the volume element and the volume 2π², (c) R = 6 and Rij = 2gij, (d) V to O(ϵ⁵) and the ratio V/VE = 1 − ϵ²/5, which is (7.20) with R = 6 and d = 3.

The induced metric is a pull-back: h = induced(E, X^1, …, X^n) pulls the metric E back through the parametrization, with the new coordinates declared last. volume(h, …) integrates √|det h| over the given limits, and series expands.

\cos^2\psi + \sin^2\psi \sin^2\theta \cos^2\phi + \sin^2\psi \sin^2\theta \sin^2\phi + \sin^2\psi \cos^2\theta simplify(eq1) 1

(a): the parametrization (7.46) satisfies the constraint.

X = coordinates(w, x, y, z) E = metric(1, 1, 1, 1) u = coordinates(\psi, \theta, \phi) h = induced(E, \cos\psi, \sin\psi \sin\theta \cos\phi, \sin\psi \sin\theta \sin\phi, \sin\psi \cos\theta) volume(h, \psi = 0 .. \pi, \theta = 0 .. \pi, \phi = 0 .. 2\pi) 2*pi**2

(b): the induced metric is diag(1, sin²ψ, sin²ψ sin²θ) — element(h) writes it out —, the volume element sin²ψ sin θ, and the volume 2π².

x = coordinates(\psi, \theta, \phi) g = metric(1, \sin^2\psi, \sin^2\psi \sin^2\theta) ricci(g) R_{{\phi}{\phi}}=2*sin(psi)**2*sin(theta)**2

(c): which is 2gφφ; and scalar(g) gives 6.

\frac{2 \pi \cdot (\epsilon - \sin\epsilon \cos\epsilon)}{\frac{4}{3} \pi \epsilon^3} series(eq1, \epsilon, 4) 1 - epsilon**2/5 + O(epsilon**4)

(d): the numerator is the volume of ψ < ϵ, which volume(h, \psi = 0 .. \epsilon, …) gives; the denominator, VE = 4πϵ³/3. (7.20), 1 − ϵ²R/(6(d + 2)) with R = 6 and d = 3, is 1 − ϵ²/5.

44. Killing fields of the sphere solves

Problem 44 · Isometries and the so(3) algebra

D. Tong, General Relativity: Example Sheet 3 (2019), question 4, p. 1 — davidtong.org.

(i) Let X and Y be two vector fields. Show that ℒX(ℒYQ) − ℒY(ℒXQ) = ℒ[X,Y]Q, when Q is either a function or a vector field. Use the Leibniz property of the Lie derivative to show that this also holds when Q is a one-form.

(ii) Demonstrate that if a Riemannian or Lorentzian manifold has two “independent” isometries then it has a third, and define what is meant by independent here.

(iii) Consider the unit sphere with metric ds² = dθ² + sin²θ dφ². Show that X = ∂/∂φ and Y = sin φ ∂/∂θ + cot θ cos φ ∂/∂φ are Killing vectors. Find a third, and show that they obey the Lie algebra of so(3).

In Sucuri: all of it: (i) on functions (the definition of the bracket), on vector fields (Jacobi) and on 1-forms (Leibniz, problem 46); (ii) the bracket of two Killing fields is Killing — “independent” means it is not a combination of the two; (iii) all of it.

With fields by components, killing(g, X) computes ℒXg and bracket(X, Y), [X, Y]. The third field is the bracket of the two — which is what (ii) asks to see in general.

x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) X = field(0, 1) Y = field(\sin\phi, \cot\theta \cos\phi) killing(g, Y) True
x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) X = field(0, 1) Y = field(\sin\phi, \cot\theta \cos\phi) bracket(X, Y) (cos(phi))*∂_\theta + (-sin(phi)*cot(theta))*∂_\phi

Z = [X, Y] = cos φ ∂θ − cot θ sin φ ∂φ.

x = coordinates(\theta, \phi) g = metric(1, \sin^2\theta) Y = field(\sin\phi, \cot\theta \cos\phi) Z = field(\cos\phi, -\cot\theta \sin\phi) bracket(Y, Z) (1)*∂_\phi

[Y, Z] = X; and bracket(Z, X) gives Y, and killing(g, Z), True: so(3).

X = tensor(1, 0) Y = tensor(1, 0) A = tensor(1, 0) B = tensor(1, 0) g = metric \forall C, D: \nabla_X g(C,D) = g([X,C],D) + g(C,[X,D]) \forall C, D: \nabla_Y g(C,D) = g([Y,C],D) + g(C,[Y,D]) \forall C, D, F: [C,[D,F]] + [D,[F,C]] + [F,[C,D]] = 0 \nabla_{[X,Y]} g(A,B) = g([[X,Y],A],B) + g(A,[[X,Y],B]) prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

(ii): ℒXg = 0 and ℒYg = 0, written out by the Leibniz rule, give ℒ[X,Y]g = 0.

45. Robertson–Walker solves

Problem 45 · Curvature of the Robertson–Walker metric

S. M. Carroll, Lecture Notes on General Relativity (1997), ch. 8 “Cosmology”, metric (8.7) on p. 219 and eqs. (8.12)–(8.14) on p. 220 — arxiv.org.

[Metric (8.7): ds2 = −dt2 + a2(t)[dr2/(1 − kr2) + r2(dθ2 + sin2θ dφ2)].]

With the metric in hand, we can set about computing the connection coefficients and curvature tensor. Setting ȧ ≡ da/dt, the Christoffel symbols are given by

Γ011 = aȧ/(1 − kr2),   Γ022 = aȧr2,   Γ033 = aȧr2 sin2θ,
Γ101 = Γ110 = Γ202 = Γ220 = Γ303 = Γ330 = ȧ/a,
Γ122 = −r(1 − kr2),   Γ133 = −r(1 − kr2) sin2θ,
Γ212 = Γ221 = Γ313 = Γ331 = 1/r,
Γ233 = −sin θ cos θ,   Γ323 = Γ332 = cot θ.   (8.12)

The nonzero components of the Ricci tensor are

R00 = −3ä/a,   R11 = (aä + 2ȧ2 + 2k)/(1 − kr2),   R22 = r2(aä + 2ȧ2 + 2k),   R33 = r2(aä + 2ȧ2 + 2k) sin2θ,   (8.13)

and the Ricci scalar is then

R = (6/a2)(aä + ȧ2 + k).   (8.14)

In Sucuri: all of it: the Christoffel symbols, the Ricci tensor and the scalar.

Carroll, ch. 8, eqs. (8.12)–(8.14) (worked example in the text): the curvature scalar of Robertson–Walker, with spatial curvature k.

x = coordinates(t, r, \theta, \phi) a = a(t) g = metric(-1, \frac{a^2}{1 - k r^2}, a^2 r^2, a^2 r^2 \sin^2\theta) scalar(g) 6*(k + a(t)*Derivative(a(t), (t, 2)) + Derivative(a(t), t)**2)/a(t)**2

R = 6(ä/a + ȧ²/a² + k/a²), eq. (8.14).

72. The hyperbolic half-space in 3D solves

Problem 72 · Negative curvature: 3D hyperbolic space

A. Guth, MIT 8.962 General Relativity (Spring 2018), Problem Set 6, Problem 4 (“Negative Curvature”), p. 4 — web.mit.edu.

Consider a 3-dimensional space given by the set of points {(x, y, z), x ∈ ℝ, y ∈ ℝ, z > 0} with the metric

ds² = (a/z²)(dx² + dy² + dz²).   (4.1)

Compute the metric, connection, Riemann curvature tensor, Ricci tensor, and Ricci scalar.

In Sucuri: all of it: the connection, the Riemann tensor, the Ricci tensor and the scalar.

MIT 8.962 (2018), problem set 6, P4: ds² = (a/z²)(dx² + dy² + dz²), with Ricci Rij = −(2/a)gij and R = −6/a.

x = coordinates(x, y, z) g = metric(\frac{a}{z^2}, \frac{a}{z^2}, \frac{a}{z^2}) scalar(g) -6/a

73. Surface of revolution solves

Problem 73 · Surfaces of revolution

C. Hirata, Ph 236 – Homework 4 (Caltech, 2011), problem 1 “Surfaces of revolution”, p. 1 — tapir.caltech.edu.

Consider a 2-dimensional surface M embedded in 3-dimensional Euclidean space ℝ3. Suppose further that it is a surface of revolution, i.e. it can be described in cylindrical coordinates (z, ϖ, φ) by an equation of the form

z = f(ϖ).

The two-dimensional surface can be written locally with coordinates (ϖ, φ).

(a) Show that the line element for the surface is ds2 = F(ϖ) dϖ2 + ϖ2 dφ2, where F(ϖ) = 1 + [f′(ϖ)]2.

(b) Compute the Christoffel symbols for this surface.

(c) Compute the Riemann tensor Rijkl for this surface (hint: use symmetries of the Riemann tensor to avoid having to do 16 tedious calculations).

(d) Show that the Ricci scalar is R = F′(ϖ) / (ϖ[F(ϖ)]2).

(e) Prove that a 2-dimensional surface has zero Riemann tensor if and only if its Ricci scalar vanishes.

(f) Use (e) to completely classify the surfaces of revolution with zero Riemann tensor.

In Sucuri: all of it: (a) the induced metric, (b)–(c) Christoffel and Riemann, (d) the scalar, (e) via the 2-dimensional identity of problem 12, (f) by solving F′ = 0: z = C₁ + C₂ϖ, planes and cones.

X = coordinates(z, \varpi, \phi) E = metric(1, 1, \varpi^2) f = f(\varpi) u = coordinates(\varpi, \phi) h = induced(E, f(\varpi), \varpi, \phi) element(h) dphi**2*varpi**2 + dvarpi**2*(Derivative(f(varpi), varpi)**2 + 1)

(a): ds² = F dϖ² + ϖ²dφ² with F = 1 + f′². christoffel(h) and riemann(h) give (b) and (c).

X = coordinates(z, \varpi, \phi) E = metric(1, 1, \varpi^2) f = f(\varpi) u = coordinates(\varpi, \phi) h = induced(E, f(\varpi), \varpi, \phi) scalar(h) 2*Derivative(f(varpi), varpi)*Derivative(f(varpi), (varpi, 2))/(varpi*(Derivative(f(varpi), varpi)**2 + 1)**2)

(d): F′/(ϖF²), with F′ = 2f′f″. This exercise found a silent error: the variable \varpi kept its backslash, and F came out constant (section VI). (e) is the identity Rμνρσ = ½R(gg − gg) in 2D, from problem 12: the Riemann tensor vanishes if and only if R vanishes.

f = f(\varpi) f' f'' = 0 solve(eq1) C1 + C2*varpi

(f): since F ≥ 1, R = 0 means F′ = 2f′f″ = 0, and z = f(ϖ) is linear: planes and cones.

74. The cone solves

Problem 74 · Geodesics, connection and curvature on the sphere and the cone

P. van Baal, Problem Set – Theory of General Relativity (Leiden), problem 5 “Geodesics, connections, curvature”, pp. 1–2 — lorentz.leidenuniv.nl.

To get more familiar to curved spaces and geodesics we now consider the two dimensional sphere and cone embedded in ℝ3 with ds2 = dx2 + dy2 + dz2:

sphere: x = a sin θ cos φ, y = a sin θ sin φ, z = a cos θ;
cone: x = r cos φ, y = r sin φ, z = r.

The following questions are to be answered for both the sphere and the cone. (Hint: see example FN.2.2.1 and problem FN.2.3.)

a) What is the metric tensor?

b) Calculate the Christoffel symbols Γabc.

c) Calculate the curvature tensor Rabcd ≡ ∂cΓabd − ∂dΓabc + ΓebdΓaec − ΓebcΓaed. Furthermore, calculate Rbd ≡ Rabad and R ≡ Raa. Give an interpretation of the results.

d) What are the geodesic equations? Determine the most general solutions. Hint for the sphere: Show that the angular momentum is conserved, and use this to determine the geodesics. Hint for the cone: If you cut open the cone, what are the geodesics? Does this agree with what you found for the curvature?

In Sucuri: all of it: for the sphere and the cone: the induced metric, Christoffel, Riemann, Ricci and R, the geodesic equations and their solutions — great circles, cot θ = M sin(φ − φ₀); and, on the cone, straight lines of the unrolled cone. The interpretation is left to the reader.

X = coordinates(x, y, z) E = metric(1, 1, 1) u = coordinates(\theta, \phi) h = induced(E, a \sin\theta \cos\phi, a \sin\theta \sin\phi, a \cos\theta) scalar(h) 2/a**2

the sphere of radius a: element(h) gives a²dθ² + a²sin²θ dφ², and R = 2/a².

X = coordinates(x, y, z) E = metric(1, 1, 1) u = coordinates(\theta, \phi) h = induced(E, a \sin\theta \cos\phi, a \sin\theta \sin\phi, a \cos\theta) orbits(h) sin(phi - phi_0)

cot θ proportional to sin(φ − φ₀): the great circles, with L the angular momentum the hint tells you to use.

X = coordinates(x, y, z) E = metric(1, 1, 1) u = coordinates(r, \phi) h = induced(E, r \cos\phi, r \sin\phi, r) scalar(h) 0

the cone: ds² = 2dr² + r²dφ², and the whole Riemann tensor vanishes.

X = coordinates(x, y, z) E = metric(1, 1, 1) u = coordinates(r, \phi) h = induced(E, r \cos\phi, r \sin\phi, r) orbits(h) Eq(-sqrt(2)/r, sqrt(2)*sqrt(kappa)*sin(sqrt(2)*(phi - phi_0)/2)/L)

1/r = c · sin((φ − φ₀)/√2): with the angle of the unrolled cone, ϑ = φ/√2, and ρ = √2 r, this is ρ sin(ϑ − ϑ₀) = constant — straight lines, as the hint says.

75. A uniform gravitational field solves

Problem 75 · Christoffel symbols in the Newtonian limit

A. Guth, MIT 8.962 General Relativity (Spring 2018), Problem Set 3, Problem 2 (“An Inertial Coordinate System in the Newtonian Limit”), part (a), p. 2 — web.mit.edu.

Consider the spacetime metric with components

g00 = −(1 + az),   g0i = 0,   gij = δij

in a region near the origin O = (0, 0, 0, 0) where az ≪ 1, and a is a constant.

(a) Compute the nonzero Christoffel coefficients of this metric.

In Sucuri: all of it: the Christoffel symbols.

MIT 8.962 (2018), problem set 3, P2(a): g00 = −(1 + az). The Christoffel symbols asked for and, beyond them, the curvature scalar.

x = coordinates(t, x, y, z) g = metric(-(1 + a z), 1, 1, 1) scalar(g) a**2/(2*(a*z + 1)**2)

and Γttz = a/(2(1 + az)), Γztt = a/2.

76. de Sitter in flat slicing solves

Problem 76 · de Sitter spacetime in flat coordinates

C. Hirata, Ph 236 – Homework 5 (Caltech, 2011), problem 1 “Cosmological constant”, p. 1 — tapir.caltech.edu.

Consider the highly symmetrical spacetime with line element given by

ds2 = (1/(H2η2)) (−dη2 + dx2 + dy2 + dz2),

where H > 0 is a constant, and in the domain η < 0. (This is known as de Sitter spacetime, and H is the Hubble constant.)

(a) Find the Einstein tensor for this spacetime.

(b) Prove that de Sitter spacetime is a vacuum solution of Einstein's equations in the presence of a positive cosmological constant for a particular choice of H. What is the relation between H and Λ?

(c) Consider an observer whose world line is given by fixed spatial coordinates (x, y, z). Explain why this trajectory is a geodesic, and show that an infinite amount of proper time elapses before the observer reaches η = 0.

In Sucuri: all of it: (a) and (b), via the Ricci tensor and the scalar; (c) the geodesic equations, which constant x, y, z solve, and the proper time, which diverges.

x = coordinates(\eta, x, y, z) g = metric(-\frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}) scalar(g) 12*H**2

and ricci(g) gives Rμν = 3H²gμν: Gμν = −3H²gμν, and the vacuum equation with Λ requires Λ = 3H².

x = coordinates(\eta, x, y, z) g = metric(-\frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}, \frac{1}{H^2 \eta^2}) geodesics(g) Derivative(x(lambda), (lambda, 2)) - 2*Derivative(eta(lambda), lambda)*Derivative(x(lambda), lambda)/eta(lambda)

(c): the equations for x, y and z only have terms with ẋ, ẏ, ż — constant x, y, z solve them, and the η equation is left with a single function.

\int_{0}^{1} \frac{1}{\sigma} d\sigma evaluate(eq1) oo

the proper time up to η = 0 is (1/H)∫dη/|η|; with σ = −η/η₀, this is what diverges.

77. de Sitter in 2D partly

Problem 77 · Two-dimensional de Sitter spacetime

J. M. Evans, Part II General Relativity, Example Sheet 1 (Cambridge, 2026), question 9, p. 3 — damtp.cam.ac.uk.

2-dimensional de Sitter space-time has the line element

ds2 = −du2 + cosh2u dφ2,

where −∞ < u < ∞ and 0 ≤ φ < 2π. Compute the Christoffel symbols and hence the geodesic equations. Verify that the equations of an affinely parametrized geodesic xα = xα(λ) can be derived from the variational principle

δ ∫ (u̇2 − cosh2u φ̇2) dλ = 0,

where u̇ = du/dλ etc. Verify from the variational principle that there are two first integrals

cosh2u φ̇ = K,   cosh2u u̇2 = K2 + L cosh2u,

along the geodesic, where K and L are constants.

Show that if K = 0 the geodesics are φ = const. If K ≠ 0, show that λ may be eliminated in favour of φ as a parameter along the geodesic and obtain the equation

v′2 = M2 − v2,

where v = tanh u, v′ = dv/dφ and M is a constant depending on L and K. Hence show that the K ≠ 0 geodesics are given by

tanh u = M sin(φ − φ0),

where φ0 is a constant. Show also that M2 > 1 for timelike geodesics, M2 = 1 for null geodesics and M2 < 1 for spacelike geodesics. Regarding u and φ as cartesian coordinates, sketch the set of geodesics starting from (0, 0).

Show from your diagram that no two such timelike geodesics will meet again, but that spacelike geodesics may recross each other. Demonstrate also that there are pairs of points which cannot be joined by a geodesic. Which, if any, of these statements would be valid in Minkowski spacetime?

In Sucuri, partly: the Christoffel symbols, the geodesic equations and the variational principle (checked), the two first integrals, and the geodesics tanh u = M sin(φ − φ₀) with M² = 1 − κ/L² — M² > 1, = 1, < 1 for timelike, null, spacelike; and R = 2. The sketch, and the arguments drawn from it — which geodesics meet again, which pairs of points cannot be joined — no.

x = coordinates(u, \phi) g = metric(-1, \cosh^2 u) orbits(g) Eq(tanh(u), sqrt(L**2 - kappa)*sin(phi - phi_0)/L)

with the substitution v = tanh u, which is Evans's: v = M sin(φ − φ₀), M² = (L² − κ)/L². geodesics(g) gives the equations and the first integrals cosh²u φ̇ = K and g(ẋ, ẋ) = κ.

x = coordinates(u, \phi) g = metric(-1, \cosh^2 u) scalar(g) 2

78. The family −f dt² + dr²/f solves

Problem 78 · Spherically symmetric Einstein–Maxwell: Reissner–Nordström and Λ

J. M. Evans, Part II General Relativity, Example Sheet 4 (Cambridge, 2026), question 6, p. 2 — damtp.cam.ac.uk.

Show that the Einstein–Maxwell equations (i.e. the Einstein equations with energy momentum tensor for an electromagnetic field Tαβ = FαγFβγ − ¼ FγδFγδgαβ) can be written

Rαβ = κ(FαγFβγ − ¼ gαβFγδFγδ).

For a line element of the form

ds2 = −f(r)c2dt2 + (1/f(r)) dr2 + r2(dθ2 + sin2θ dφ2) ,

the only non-zero components of the Ricci tensor are given by

Rtt/(c2f) = −f Rrr = ½ f″ + f′/r , Rθθ = Rφφ/sin2θ = 1 − r f′ − f .

In the case

Ftr = −Frt = Q/r2 and Fαβ = 0 otherwise,

show that a solution can be found that reduces to the Schwarzschild solution when Q = 0.

Find an analogous solution in the case Rαβ = Λgαβ.

In Sucuri: all of it: the form Rαβ = κ(…), via the trace in 4D; the Ricci components the problem statement gives, checked; f for Einstein–Maxwell, which with Q = 0 is Schwarzschild; and f for Rαβ = Λgαβ. (With c = 1.)

in_chart(eq) evaluates an indexed expression component by component, in the chart of the declared metric: the metric, forms in a chart, fields, and Riemann, Ricci and the scalar through the Christoffel symbols — in the convention that riemann(eq) and ricci(eq) read.

\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices(4) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) F = tensor(0, 2, antisymmetric) \kappa = constant R_{\alpha\beta} - \frac{1}{2} g_{\alpha\beta} R = \kappa \cdot (F_{\alpha\gamma} F_\beta{}^\gamma - \frac{1}{4} g_{\alpha\beta} F_{\gamma\delta} F^{\gamma\delta}) R_{\alpha\beta} = \kappa \cdot (F_{\alpha\gamma} F_\beta{}^\gamma - \frac{1}{4} g_{\alpha\beta} F_{\gamma\delta} F^{\gamma\delta}) prove(eq4, eq3) proved from eq3

the trace, in 4D, gives −R = 0 — T is traceless —, and hence Rαβ = κTαβ.

\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita x = coordinates(t, r, \theta, \phi) f = f(r) g = metric(-f, \frac{1}{f}, r^2, r^2 \sin^2\theta) \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) R_{\mu\nu} in_chart(eq3) -r*Derivative(f(r), r) - f(r) + 1

the components in the problem statement: Rtt/f = −fRrr = ½f″ + f′/r, Rθθ = Rφφ/sin²θ = 1 − rf′ − f.

f = f(r) 1 - r f'(r) - f(r) = \frac{\kappa Q^2}{2 r^2} solve(eq1) Eq(f(r), C1/r + Q**2*kappa/(2*r**2) + 1)

the θθ component, with Tθθ = Q²/(2r²) (from in_chart, with F = forma(Q/r² dt∧dr)). With C1 = −2m, Q = 0 gives Schwarzschild.

\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita x = coordinates(t, r, \theta, \phi) g = metric(-(1 - \frac{2m}{r} + \frac{\kappa Q^2}{2 r^2}), \frac{1}{1 - \frac{2m}{r} + \frac{\kappa Q^2}{2 r^2}}, r^2, r^2 \sin^2\theta) \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) F = form(\frac{Q}{r^2} dt \wedge dr) \kappa = constant R_{\alpha\beta} = \kappa \cdot (F_{\alpha\gamma} F_\beta{}^\gamma - \frac{1}{4} g_{\alpha\beta} F_{\gamma\delta} F^{\gamma\delta}) in_chart(eq3) True

the solution checks out in all four components, not just θθ.

f = f(r) 1 - r f'(r) - f(r) = \Lambda r^2 solve(eq1) Eq(f(r), C1/r - Lambda*r**2/3 + 1)
\alpha, \beta, \gamma, \delta, \mu, \nu, \rho, \sigma = indices V = tensor(1, 0) \nabla = levi-civita x = coordinates(t, r, \theta, \phi) g = metric(-(1 - \frac{2m}{r} - \frac{\Lambda r^2}{3}), \frac{1}{1 - \frac{2m}{r} - \frac{\Lambda r^2}{3}}, r^2, r^2 \sin^2\theta) \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) R_{\alpha\beta} = \Lambda g_{\alpha\beta} in_chart(eq3) True

Schwarzschild–de Sitter.

79. Reissner–Nordström solves

Problem 79 · Ricci curvature of Reissner–Nordström

T. Dray, MTH 437/537 – HW #6 (Oregon State, Spring 2024), problem 1, p. 1 — sites.science.oregonstate.edu.

Consider the Reissner–Nordström geometry, with line element

−f(r) dt2 + dr2/f(r) + r2(dθ2 + sin2θ dφ2),

where f(r) = 1 − 2m/r + q2/r2. Determine the components Rij of the Ricci curvature.

You may use any formalism you wish, and any coordinate system you wish. Tedious computation by hand is discouraged, as is trying to fully typeset your results. Other options include using computer algebra and/or using or adapting a computation published elsewhere. Such strategies must of course be documented. However, make sure you understand — and clearly describe — the conventions being used! One reasonable starting point would be the curvature 2-forms, and the relations

Ωij = ½ Rijkl σk ∧ σl,   Rij = Rmimj,

where repeated indices are summed over.

In Sucuri: all of it: the components Rij.

Dray, MTH 437, HW6 #1: with f = 1 − 2m/r + q²/r², Rθθ = q²/r², and R = 0.

x = coordinates(t, r, \theta, \phi) g = metric(-(1 - \frac{2m}{r} + \frac{q^2}{r^2}), \frac{1}{1 - \frac{2m}{r} + \frac{q^2}{r^2}}, r^2, r^2 \sin^2\theta) ricci(g) R_{{\theta}{\theta}}=q**2/r**2

80. Global AdS₃ solves

Problem 80 · Metric of AdS3

T. Hartman, Lectures on Quantum Gravity and Black Holes (Cornell, 2015), §9.1 “Exercise: Metric of AdS3”, eqs. (9.1)–(9.4), pp. 97–98 — hartmanhep.net.

Anti-de Sitter space is a constant-negative-curvature spacetime. It is the maximally symmetric solution of Einstein's equation with a negative cosmological constant. AdSD can be realized as a hyperboloid embedded in a D + 1-dimensional geometry. In this section we will talk about AdS3, which is the hyperboloid

XAXA = −ℓ2   (9.1)

where A = 0, 1, 2, 3 is an index in the space Minkowski2 × Minkowski2, with metric

HAB dXAdXB = −dX02 + dX12 + dX22 − dX32.   (9.2)

To find intrinsic coordinates on AdS3, we just need to solve (9.1). One way to solve this equation is by

X0 = ℓ cosh ρ cos t,  X1 = ℓ sinh ρ sin φ,  X2 = ℓ sinh ρ cos φ,  X3 = ℓ cosh ρ sin t.   (9.3)

(a) Check that this solves (9.1), and use (9.2) to find the induced metric on the hyperboloid.

Answer: ds2 = ℓ2(−cosh2ρ dt2 + dρ2 + sinh2ρ dφ2).   (9.4)

These are global coordinates on AdS3. Although on the hyperboloid (9.1) we can see from (9.3) that t is a periodic coordinate, when we say “AdS3” we will always mean the space in which t is “unwrapped”, t ∈ (−∞, ∞) (the universal covering space of the hyperboloid).

(b) Find the cosmological constant in terms of the AdS radius ℓ.

In Sucuri: all of it: (a) the hyperboloid and the induced metric (9.4); (b) Λ = −1/ℓ², via the Ricci tensor.

-(\ell \cosh\rho \cos t)^2 + (\ell \sinh\rho \sin\phi)^2 + (\ell \sinh\rho \cos\phi)^2 - (\ell \cosh\rho \sin t)^2 simplify(eq1) -ell**2

(a): the parametrization (9.3) lies on the hyperboloid (9.1).

X = coordinates(X_0, X_1, X_2, X_3) H = metric(-1, 1, 1, -1) u = coordinates(t, \rho, \phi) g = induced(H, \ell \cosh\rho \cos t, \ell \sinh\rho \sin\phi, \ell \sinh\rho \cos\phi, \ell \cosh\rho \sin t) element(g) dphi**2*ell**2*sinh(rho)**2 + drho**2*ell**2 - dt**2*ell**2*cosh(rho)**2

the induced metric is (9.4).

x = coordinates(t, \rho, \phi) g = metric(-\ell^2 \cosh^2\rho, \ell^2, \ell^2 \sinh^2\rho) scalar(g) -6/ell**2

(b): with Rμν = −(2/ℓ²)gμν, Gμν = gμν/ℓ², and Gμν + Λgμν = 0 gives Λ = −1/ℓ².

81. Hyperbolic space H³ partly

Problem 81 · Geometry of hyperbolic space

C. Hirata, Ph 236 – Homework 16 (Caltech, 2012), problem 1 “Geometry of hyperbolic space”, p. 1 — tapir.caltech.edu.

Consider the hyperbolic space H3 with unit negative curvature,

ds2 = dχ2 + sinh2χ (dθ2 + sin2θ dφ2).

(a) Show that this is the geometry that one obtains by taking 4-dimensional Minkowski space M4 and restricting to the spacelike 3-surface Σ ⊂ M4 given by −t2 + x2 + y2 + z2 = −1. What is the mapping Φ : H3 → Σ of the coordinates (χ, θ, φ) to (t, x, y, z)?

(b) Using your knowledge of infinitesimal Lorentz transformations, find the 6 Killing fields corresponding to Lorentz rotations and boosts, and express them as contravariant vectors in the χθφ coordinate system.

(c) Show that given two points P, Q ∈ H3 that the separation s between them (i.e. length of the shortest¹ geodesic that connects them) is given by cosh s = −Φ(P)·Φ(Q), where the dot product is that defined in M4.

(d) Use the result of (c) to prove that in the hyperbolic space, the law of cosines is given by cosh c = cosh a cosh b − sinh a sinh b cos γ, where we denote the side lengths of a “triangle” (i.e. set of 3 geodesics connecting the 3 vertices) as a, b, and c, and the angles opposite them as α, β, and γ.

(e) Taylor expand the result from (d) and show that if a, b, and c are ≪ 1 then the usual law of cosines is recovered.

¹ In hyperbolic space, the statement about the shortest geodesic is redundant: there is exactly one geodesic connecting any two points. This is actually true for any simply connected space of Euclidean signature (+ + +) and negative-semidefinite curvature, i.e. where Rαβγδξαηβξγηδ ≤ 0 for any vectors ξ and η.

In Sucuri, partly: (a) the induced metric and the map Φ, and (b) the six Killing fields — rotations and boosts of M⁴ restricted to the hyperboloid, and checked; (c) the distance, which needs invariance under the isometries, and with it (d) and (e), no.

X = coordinates(t, x, y, z) M = metric(-1, 1, 1, 1) u = coordinates(\chi, \theta, \phi) h = induced(M, \cosh\chi, \sinh\chi \sin\theta \cos\phi, \sinh\chi \sin\theta \sin\phi, \sinh\chi \cos\theta) element(h) dchi**2 + dphi**2*sin(theta)**2*sinh(chi)**2 + dtheta**2*sinh(chi)**2

(a): Φ(χ, θ, φ) = (cosh χ, sinh χ sin θ cos φ, sinh χ sin θ sin φ, sinh χ cos θ), which lies on −t² + x² + y² + z² = −1.

x = coordinates(\chi, \theta, \phi) g = metric(1, \sinh^2\chi, \sinh^2\chi \sin^2\theta) scalar(g) -6

curvature −1: Rij = −2gij, R = −6.

X = coordinates(t, x, y, z) M = metric(-1, 1, 1, 1) B = field(x, t, 0, 0) u = coordinates(\chi, \theta, \phi) h = induced(M, \cosh\chi, \sinh\chi \sin\theta \cos\phi, \sinh\chi \sin\theta \sin\phi, \sinh\chi \cos\theta) restrict(B, h) killing(h, B) True

(b): the boost x∂t + t∂x of M⁴, restricted to Σ, is sin θ cos φ ∂χ + cos φ cos θ coth χ ∂θ − sin φ/(sin θ tanh χ) ∂φ, and it is Killing on H³. The other two boosts and the three rotations likewise; a translation, restrict refuses: it is not tangent.

82. The 4-sphere solves

Problem 82 · Ricci curvature of the N-sphere

G. 't Hooft, Introduction to General Relativity (Utrecht, 2013), ch. 14 “The Robertson-Walker metric”, exercise after eq. (14.8), with eq. (14.9), p. 60 — webspace.science.uu.nl.

[Context: for the three-dimensional isotropic space dω2 = B(ϱ)dϱ2 + ϱ2(dθ2 + sin2θ dϕ2) (14.2), one requires Rij = λgij (14.5), which leads to B = 1/(1 − ½λϱ2) (14.8).]

Exercise: show that with ϱ = √(2/λ) sin ψ, this gives the metric of the 3-sphere, in terms of its three angular coordinates ψ, θ, ϕ. Indeed, the metric for an N-sphere can be written as

dωN2 = dψN2 + sin2ψN dωN−12,   Rij(N) = λNgij;  λN = N − 1.   (14.9)

In Sucuri: all of it: the metric of the 3-sphere via ϱ = √(2/λ) sin ψ, and Rij = (N − 1)gij.

X = coordinates(\varrho, \theta, \phi) g = metric(\frac{1}{1 - \frac{1}{2} \lambda \varrho^2}, \varrho^2, \varrho^2 \sin^2\theta) u = coordinates(\psi, \theta, \phi) h = induced(g, \sqrt{\frac{2}{\lambda}} \sin\psi, \theta, \phi) element(h) 2*dphi**2*sin(psi)**2*sin(theta)**2/lambda + 2*dpsi**2/lambda + 2*dtheta**2*sin(psi)**2/lambda

the exercise: with ϱ = √(2/λ) sin ψ, (14.2) with (14.8) is (2/λ)(dψ² + sin²ψ dΩ²), the 3-sphere. A change of coordinates is the same pull-back as an embedding.

x = coordinates(\chi, \psi, \theta, \phi) g = metric(1, \sin^2\chi, \sin^2\chi \sin^2\psi, \sin^2\chi \sin^2\psi \sin^2\theta) scalar(g) 12

and (14.9) for N = 4: Rij = 3gij, R = 12.

85. Flamm's paraboloid solves

Problem 85 · The shape of the Schwarzschild geometry (Flamm's paraboloid)

T. Dray, MTH 437/537 – HW #5 (Oregon State, Spring 2024), problem 1 “Shape of Schwarzschild geometry”, p. 1 — sites.science.oregonstate.edu.

(a) Plot the graph of 8mr = h2 + 16m2 with r as the horizontal coordinate. Equivalently, solve this equation for h, and plot h as a function of r.

(b) Express arclength along this graph in terms of r (and dr). Assume that (r, h) are rectangular, Euclidean coordinates.

(c) Consider the surface of revolution obtained by rotating your curve about the h-axis. What is the line element for this surface? Compare your answer with the Schwarzschild geometry.

(d) Find the Gaussian curvature of this surface. Is it positive or negative?

In Sucuri: all of it: (a) h(r), (b) the arclength, (c) the surface of revolution, which is the equatorial slice of Schwarzschild, and (d) K = −m/r³ < 0. The plot is left to the reader.

8 m r = h^2 + 16 m^2 solve(eq1) 2*sqrt(2)*sqrt(m*(-2*m + r))

(a): h = ±√(8m(r − 2m)).

X = coordinates(r, h) E = metric(1, 1) u = coordinates(r) c = induced(E, r, \sqrt{8 m r - 16 m^2}) element(c) dr**2*r/(-2*m + r)

(b): ds² = dr²/(1 − 2m/r) along the graph.

X = coordinates(r, \phi, h) E = metric(1, r^2, 1) u = coordinates(r, \phi) g = induced(E, r, \phi, \sqrt{8 m r - 16 m^2}) scalar(g) -2*m/r**3

(c): the surface has ds² = dr²/(1 − 2m/r) + r²dφ², the equatorial slice of Schwarzschild; (d): K = R/2 = −m/r³, negative.

83. The Ricci tensor of a weak field solves

Problem 83 · Ricci scalar of a special metric

Y. Ali-Haïmoud, General Relativity (NYU, Fall 2019), Homework 6, Exercise 3, p. 1 — cosmo.nyu.edu.

At linear order in Φ, compute the 10 components of the Ricci tensor, as well as the Ricci scalar of the metric

ds² = −[1 + 2Φ(x⃗)] dt² + [1 − 2Φ(x⃗)] δij dxidxj.   (3)

Note that Φ(x⃗) is assumed to depend only on the spatial coordinates (i.e. not on t).

In Sucuri: all of it: the ten components and the scalar, to first order.

ricci(g, \Phi) computes to first order in Φ: it replaces Φ by εΦ, does the exact computation, truncates the series in ε after order 1 and sets ε = 1.

x = coordinates(t, x, y, z) \Phi = \Phi(x, y, z) g = metric(-(1 + 2\Phi), 1 - 2\Phi, 1 - 2\Phi, 1 - 2\Phi) ricci(g, \Phi) R_{{t}{t}}=Derivative(Phi(x, y, z), (x, 2)) + Derivative(Phi(x, y, z), (y, 2)) + Derivative(Phi(x, y, z), (z, 2))

Rtt = ∇²Φ and Rij = δij∇²Φ; the six off-diagonal ones vanish at this order.

x = coordinates(t, x, y, z) \Phi = \Phi(x, y, z) g = metric(-(1 + 2\Phi), 1 - 2\Phi, 1 - 2\Phi, 1 - 2\Phi) scalar(g, \Phi) 2*Derivative(Phi(x, y, z), (x, 2)) + 2*Derivative(Phi(x, y, z), (y, 2)) + 2*Derivative(Phi(x, y, z), (z, 2))

R = 2∇²Φ.

84. Connection forms on the sphere solves

Problem 84 · Spherical coordinates, III

T. Dray, MTH 434/534 – HW #6 (Oregon State, Winter 2024), problem 1, p. 1 — sites.science.oregonstate.edu.

Consider the sphere of radius r, in spherical coordinates (θ, φ), with line element ds² = r²(dθ² + sin²θ dφ²). [The original prints r²dθ² + sin²θ dφ², without the parentheses.]

(a) Find the connection 1-forms ωij in this basis.

(b) Compute Ωij = dωij + ωik ∧ ωkj for i, j = 1, 2 (and where there is an implicit sum over k).

(c) (Optional for MTH 434/Required for MTH 534.) Compare your answers (and your computations) with those from the previous homework assignment.

In Sucuri: all of it: (a) and (b); the comparison in (c) is left to the reader.

x = coordinates(\theta, \phi) g = metric(r^2, r^2 \sin^2\theta) cartan(g) sin(theta)*dtheta∧dphi

in Dray's ordering, e1 = r dθ and e2 = r sin θ dφ — here e0, e1 —: ω12 = −cos θ dφ, and Ω12 = (1/r²) e1∧e2 = sin θ dθ∧dφ.

III. Index-free notation

18 and 19. The Lie bracket solves

Problem 18 · The Lie bracket makes 𝔛(M) a Lie algebra

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 2.2.3, §2.2, p. 19 — math.nagoya-u.ac.jp.

Context: The set 𝔛(M) has also the structure of a Lie algebra, meaning that there exists a Lie bracket [ , ] : 𝔛(M) × 𝔛(M) → 𝔛(M) satisfying:

(i) It is linear in each argument;
(ii) It is antisymmetric: [Y, X] = −[X, Y];
(iii) The Jacobi identity holds: [X, [Y, Z]] + [Y, [Z, X]] + [Z, [X, Y]] = 0.

[…] by defining [X, Y]p(f) := Xp(Yf) − Yp(Xf), this verifies the linearity property and the Leibniz rule. As a consequence, [X, Y]p is an element of Tp(M) for each p ∈ M.

Exercise 2.2.3. Check that 𝔛(M) endowed with the product [X, Y] is a Lie algebra, see [2, p. 152–153].

In Sucuri: all of it: the properties of the bracket, index-free.

Problem 19 · Bracket of fields multiplied by functions

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 2.2.4, §2.2, p. 19 — math.nagoya-u.ac.jp.

Exercise 2.2.4. For any X, Y ∈ 𝔛(M) and f, g ∈ C∞(M), check the following equality:

[fX, gY] = fg[X, Y] + f(Xg)Y − g(Yf)X.

In Sucuri: all of it: the formula.

Reall §1.10 and Richard, ex. 2.2.3–2.2.4: antisymmetry, additivity, Leibniz, Jacobi; and [fX, gY] = fg[X,Y] + f X(g) Y − g Y(f) X.

X = tensor(1, 0) Y = tensor(1, 0) [f X, g Y] = f g [X, Y] + f \nabla_X g \, Y - g \nabla_Y f \, X prove(eq1) proved from no hypothesis

\nabla_X g, with g a function, is the directional derivative X(g). The Leibniz rule of the bracket is something the engine knows on its own.

The Jacobi identity between vector fields is not built in — what is built in is its form on functions, [A,B](f) = A(B(f)) − B(A(f)), which is the definition of the bracket. And that is how Jacobi comes out, applied to a function h:

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \nabla_{[X,[Y,Z]]} h + \nabla_{[Y,[Z,X]]} h + \nabla_{[Z,[X,Y]]} h = 0 prove(eq1) proved from no hypothesis

20. The Lie derivative of the bracket solves

Problem 20 · Properties of the Lie derivative of vector fields

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 2.3.12, §2.3, p. 23 — math.nagoya-u.ac.jp.

Context (Lemma 2.3.11): For any X, Y ∈ 𝔛(M), one has LX(Y) = [X, Y].

Exercise 2.3.12. Prove (and understand) the following equality, for any X, Y, Z ∈ 𝔛(M) and f ∈ C∞(M):

(i) LX[Y, Z] = [LXY, Z] + [Y, LXZ],
(ii) LX ∘ LY − LY ∘ LX = L[X,Y],
(iii) LX(fY) = (LXf)Y + f(LXY).

In Sucuri: all of it: (i) and (ii), with Jacobi as a hypothesis — ℒX is the bracket, and both are Jacobi rewritten —, and (iii), the Leibniz rule of the bracket, with no hypothesis.

Richard, ex. 2.3.12: with ℒXY = [X,Y], ℒX is a derivation of the bracket — which is Jacobi, rewritten. With Jacobi as a hypothesis:

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \forall A, B, C: [A,[B,C]] + [B,[C,A]] + [C,[A,B]] = 0 [X, [Y, Z]] = [[X, Y], Z] + [Y, [X, Z]] prove(eq2, eq1) eq1[A→X, B→Y, C→Z]

without the hypothesis, the same call does not prove it, and says that no hypothesis talks about those nested brackets.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \forall A, B, C: [A,[B,C]] + [B,[C,A]] + [C,[A,B]] = 0 [X, [Y, Z]] - [Y, [X, Z]] = [[X, Y], Z] prove(eq2, eq1) eq1[A→X, B→Y, C→Z]

(ii): (ℒXℒY − ℒYℒX)Z = ℒ[X,Y]Z.

X = tensor(1, 0) Y = tensor(1, 0) [X, f Y] = \nabla_X f \, Y + f [X, Y] prove(eq1) proved from no hypothesis

(iii): ℒX(fY) = (ℒXf)Y + fℒXY. Without the X(f)Y term, it does not come out.

21. Torsion and curvature are tensors solves

Problem 21 · Torsion and curvature are C∞(M)-linear

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 5.4.3 (after Definition 5.4.2), §5.4, p. 49 — math.nagoya-u.ac.jp.

Context (Definition 5.4.2): For a Riemannian manifold M and an affine connection ∇ one sets for any X, Y ∈ 𝔛(M)

T(X, Y) := ∇XY − ∇YX − [X, Y] ∈ 𝔛(M)  and  R(X, Y) := ∇X∇Y − ∇Y∇X − ∇[X,Y] ∈ End(𝔛(M)),

the former being called the torsion while the latter is called the curvature of the connection. It is then interesting to observe that the map (X, Y) ↦ T(X, Y) is C∞(M)-linear in both arguments, and that the map (X, Y, Z) ↦ R(X, Y)Z is C∞(M)-linear in the three arguments.

Exercise 5.4.3. Show the C∞(M)-linearity mentioned above, see also [8, Prop. 6.3].

In Sucuri: all of it: the C∞(M)-linearity of T and of R.

Richard, ex. 5.4.3 (also Tong, sheet 2, Q2): T(fX, gY) = fg T(X,Y) and R(fX, gY)(hZ) = fgh R(X,Y)Z. Written out from the definitions, assuming nothing:

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \nabla_{f X} \nabla_{g Y} (h Z) - \nabla_{g Y} \nabla_{f X} (h Z) - \nabla_{[f X, g Y]} (h Z) = f g h \cdot (\nabla_X \nabla_Y Z - \nabla_Y \nabla_X Z - \nabla_{[X,Y]} Z) prove(eq1) proved from no hypothesis

all the derivatives of f, g and h that Leibniz produces cancel — including the second ones, via the bracket acting on h. The torsion one comes out the same way.

22. The difference of two connections solves

Problem 22 · The difference of two connections is a tensor

H. S. Reall, Part III General Relativity (2022), §3.1 “Introduction”, exercise right after eq. (3.10), p. 33 (p. 34 of the PDF) — damtp.cam.ac.uk.

Let ∇ and ∇̃ be two different connections on M. Show that ∇ − ∇̃ is a (1, 2) tensor field. You can do this either from the definition of a connection, or from the transformation law for the connection components.

In Sucuri: all of it: from the definition of a connection: (∇ − ∇̃)XY is linear over functions in both slots — which makes it a (1, 2) tensor —, and, if both are torsion-free, symmetric.

An accent on ∇ — \tilde\nabla, \hat\nabla, \bar\nabla — is another connection. About it the engine only knows what holds for every connection: linear over functions in the direction, Leibniz in the operand, and ∇̃Xf = X(f) on a function.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \nabla_{f X + Z} (h Y) - \tilde\nabla_{f X + Z} (h Y) = f h \cdot (\nabla_X Y - \tilde\nabla_X Y) + h \cdot (\nabla_Z Y - \tilde\nabla_Z Y) prove(eq1) proved from no hypothesis

the two Leibniz X(h)Y terms cancel: the difference is linear over functions in both slots, and is a tensor — unlike each ∇.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \nabla_X (h Y) - \tilde\nabla_X (h Y) = h \nabla_X Y prove(eq1) I found no combination of the hypotheses that gives this

and the engine does not confuse the two: without ∇̃ = 0 declared, this does not come out.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \tilde\nabla_X Y - \tilde\nabla_Y X = [X, Y] \nabla_X Y - \nabla_Y X = [X, Y] (\nabla_X Y - \tilde\nabla_X Y) - (\nabla_Y X - \tilde\nabla_Y X) = 0 prove(eq3, eq1, eq2) proved from eq1, eq2

with both torsion-free, the difference is symmetric — like the δΓμνρ Reall uses when varying the Einstein–Hilbert action.

23. The Koszul formula solves

Problem 23 · Koszul formula

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 5.4.5, §5.4, p. 49 — math.nagoya-u.ac.jp.

Context: a connection is called torsion free if T(X, Y) = 0 for all X, Y, and is called Riemannian (or Levi-Civita) if, in addition,

X⟨Y1, Y2⟩ = ⟨∇XY1, Y2⟩ + ⟨Y1, ∇XY2⟩.   (5.4.1)

Exercise 5.4.5 (Koszul formula). For X, Y, Z ∈ 𝔛(M) set

⟨∇XY, Z⟩ := ½ { X⟨Y, Z⟩ + Y⟨Z, X⟩ − Z⟨X, Y⟩ − ⟨X, [Y, Z]⟩ + ⟨Y, [Z, X]⟩ + ⟨Z, [X, Y]⟩ }.

Check that the connection defined by this relation is torsion free and satisfies (5.4.1). This equality is called Koszul formula and shows the existence of a Riemannian connection.

In Sucuri: all of it: both directions: the one asked for — with ∇ defined by the formula, it is torsion-free and metric-compatible — and the converse, uniqueness. From g(W, Z) = 0 for all Z follows W = 0: that is the metric being non-degenerate.

Richard, ex. 5.4.5 (and Reall §3.2): the Koszul formula, from the two Levi-Civita conditions — metric compatibility and zero torsion. Richard asks for the opposite direction (define ∇ by the formula, check the two conditions); Sucuri proves this one, which is uniqueness.

g = metric X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \forall A, B, C: \nabla_A g(B,C) = g(\nabla_A B, C) + g(B, \nabla_A C) \forall A, B: \nabla_A B - \nabla_B A = [A,B] 2 g(\nabla_X Y, Z) = \nabla_X g(Y,Z) + \nabla_Y g(Z,X) - \nabla_Z g(X,Y) - g(X,[Y,Z]) + g(Y,[Z,X]) + g(Z,[X,Y]) prove(eq3, eq1, eq2) g(eq2[A→X, B→Y], Z)

three instances of compatibility and three of zero torsion, the latter placed inside g(□, ·). Without zero torsion, it does not come out.

The direction Richard asks for: with the Koszul formula as the definition of ∇ — a hypothesis for all fields —, the two Levi-Civita conditions come out.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) g = metric \forall A, B, C: 2 g(\nabla_A B, C) = \nabla_A g(B,C) + \nabla_B g(C,A) - \nabla_C g(A,B) - g(A,[B,C]) + g(B,[C,A]) + g(C,[A,B]) g(\nabla_X Y, Z) - g(\nabla_Y X, Z) - g([X,Y], Z) = 0 prove(eq2, eq1) proved from eq1

torsion-free: g(∇XY − ∇YX − [X, Y], Z) = 0 for all Z.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) g = metric \forall A, B, C: 2 g(\nabla_A B, C) = \nabla_A g(B,C) + \nabla_B g(C,A) - \nabla_C g(A,B) - g(A,[B,C]) + g(B,[C,A]) + g(C,[A,B]) \nabla_X g(Y,Z) = g(\nabla_X Y, Z) + g(Y, \nabla_X Z) prove(eq2, eq1) proved from eq1

metric-compatible, (5.4.1). With the sign flipped, it does not come out.

24. The first Bianchi identity solves

Problem 24 · First Bianchi identity

S. Richard, Differential geometry (Nagoya University, 2024), Exercise 6.2.2 (on Lemma 6.2.1), §6.2, p. 60 — math.nagoya-u.ac.jp.

Context (Lemma 6.2.1): Assume that the connection on smooth manifold is torsion free, then the following identity holds for the curvature: for any X, Y, Z ∈ 𝔛(M),

R(X, Y)Z + R(Y, Z)X + R(Z, X)Y = 0,

where 0 means the 0-vector field.

Exercise 6.2.2. Prove Bianchi identity getting inspiration from [4, Prop. 4.1.3].

In Sucuri: all of it: the identity, from the definitions and zero torsion.

Richard, ex. 6.2.2: torsion-free, R(X,Y)Z + R(Y,Z)X + R(Z,X)Y = 0.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) R = curvature \forall A, B, W: R(A,B)W = \nabla_A \nabla_B W - \nabla_B \nabla_A W - \nabla_{[A,B]} W \forall A, B: \nabla_A B - \nabla_B A = [A,B] \forall A, B, C: [A,[B,C]] + [B,[C,A]] + [C,[A,B]] = 0 R(X,Y)Z + R(Y,Z)X + R(Z,X)Y = 0 prove(eq4, eq1, eq2, eq3) eq3[A→X, B→Y, C→Z]

ten steps: the three definitions of R, six instances of zero torsion (three of them inside ∇) and Jacobi once.

25. The second Bianchi identity, and Schur's lemma solves

Problem 25 · Second Bianchi identity and Einstein manifolds

A. G. Kovalev, Part III: Riemannian Geometry — Example Sheets (Cambridge, Lent 2017), Example Sheet 1, Question 2, p. 1 — dec41.user.srcf.net.

2. (i) Let M be a Riemannian manifold. Show that the Levi–Civita covariant derivative of R(X, Y) ∈ Γ(End TM) is given by

∇ZR(X, Y) = [∇Z, R(X, Y)] − R(∇ZX, Y) − R(X, ∇ZY).

Deduce from this a version of the second Bianchi identity for the Levi–Civita connection

∇XR(Y, Z) + ∇YR(Z, X) + ∇ZR(X, Y) = 0.   (∗)

(ii) When dim M ≥ 3, show, using (∗), that if Ric = f g for some smooth function f, then f is constant (M then is said to be an Einstein manifold).
(You might like to consider a map δ : Γ(Sym2 T*M) → Γ(T*M) = Ω1(M) defined by (δh)(X) = −∑i=1n (∇eih)(ei, X), where {ei} is any local orthonormal frame field on M, and put h = Ric.)

In Sucuri: all of it: (i) the second Bianchi identity, from the Ricci identity — the formula for ∇ZR(X, Y) is the Leibniz rule for ∇ of a tensor, and the deduction comes out in indices; (ii) Ric = fg ⇒ ∇f = 0, with the condition d ≠ 2 stated by the proof itself.

(i) in indices: the Ricci identity, ∇μ∇νV − ∇ν∇μV = RV, is something the canonical form already knows — on its own, it simplifies to zero —, but its ∇ it does not: prove differentiates it, sums the three cyclic permutations, and closes with the first Bianchi identity.

\mu, \nu, \rho, \sigma, \lambda = indices V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) (\nabla_\lambda R^\rho{}_{\sigma\mu\nu} + \nabla_\mu R^\rho{}_{\sigma\nu\lambda} + \nabla_\nu R^\rho{}_{\sigma\lambda\mu}) V^\sigma = 0 prove(eq2, eq1) proved from eq1, and from Bianchi

for all V: ∇[λRρ|σ|μν] = 0. With one sign flipped in the sum, it does not come out.

(ii) Schur's lemma: with the contracted Bianchi identity (exercise 10) and Rμν = fgμν, ∇f = 0 — and the combination divides by d − 2, which the proof states. In dimension 2, every Ricci is fg (problem 12), and f need not be constant.

\mu, \nu, \rho, \sigma, \lambda = indices(d) V = tensor(1, 0) \nabla = levi-civita g = metric \nabla_\mu \nabla_\nu V^\rho - \nabla_\nu \nabla_\mu V^\rho = R^\rho{}_{\sigma\mu\nu} V^\sigma R = riemann(eq1) R_{\mu\nu} = R^\rho{}_{\mu\rho\nu} R = ricci(eq2) \nabla^\mu R_{\mu\nu} = \frac{1}{2} \nabla_\nu R R_{\mu\nu} = f g_{\mu\nu} \nabla_\nu f = 0 prove(eq5, eq3, eq4) holds if d - 2 ≠ 0

26. Lie group with a bi-invariant metric partly

Problem 26 · Lie group with a bi-invariant metric

A. G. Kovalev, Part III: Riemannian Geometry — Example Sheets (Cambridge, Lent 2017), Example Sheet 1, Question 5 (pp. 1–2 of the PDF) — dec41.user.srcf.net.

5. Let G be a Lie group endowed with a Riemannian metric g which is left and right invariant and let X, Y, Z be left invariant vector fields of G.

(i) Show that g([X, Y], Z) + g(Y, [X, Z]) = 0. (Consider the flow of X.)
(ii) Show that ∇XX = 0. (Hint: consider g(Y, ∇XX).)
(iii) Show that ∇XY = ½[X, Y].
(iv) Prove that R(X, Y)Z = ¼[[X, Y], Z].
(v) Suppose that X and Y are orthonormal, and let K(σ) be the sectional curvature of the 2-plane σ spanned by X and Y. Prove that

K(σ) = ¼ |[X, Y]|g2.

In Sucuri, partly: (ii)–(v), with part (i) — g([X,Y], Z) + g(Y, [X,Z]) = 0 — as a hypothesis, and g(Y, Z) constant for left-invariant fields; (iv) and (v) in Kovalev's sign convention, stated as a definition. (i) itself, via the flow of X, no.

Left-invariant fields of a left-invariant metric have constant inner products, ∇Ag(B, C) = 0; right invariance, differentiated, is part (i). With these two hypotheses and the Koszul formula, (ii) and (iii) come out directly — without the trick of applying ∇XX = 0 to X + Y.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) g = metric \forall A, B, C: 2 g(\nabla_A B, C) = \nabla_A g(B,C) + \nabla_B g(C,A) - \nabla_C g(A,B) - g(A,[B,C]) + g(B,[C,A]) + g(C,[A,B]) \forall A, B, C: \nabla_A g(B,C) = 0 \forall A, B, C: g([A,B], C) + g(B, [A,C]) = 0 g(\nabla_X X, Z) = 0 prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

(ii): g(∇XX, Z) = 0 for all Z, hence ∇XX = 0.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) g = metric \forall A, B, C: 2 g(\nabla_A B, C) = \nabla_A g(B,C) + \nabla_B g(C,A) - \nabla_C g(A,B) - g(A,[B,C]) + g(B,[C,A]) + g(C,[A,B]) \forall A, B, C: \nabla_A g(B,C) = 0 \forall A, B, C: g([A,B], C) + g(B, [A,C]) = 0 g(\nabla_X Y, Z) = \frac{1}{2} g([X,Y], Z) prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

(iii): ∇XY = ½[X, Y].

(iv) depends on the sign of R. Kovalev uses R(X, Y) = ∇[X,Y] − [∇X, ∇Y], the opposite of Carroll and Reall; with that definition declared, ¼[[X, Y], Z] comes out — with Carroll's it would be −¼, and Sucuri says so.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) R = curvature \forall A, B, W: R(A,B)W = \nabla_B \nabla_A W - \nabla_A \nabla_B W + \nabla_{[A,B]} W \forall A, B: \nabla_A B = \frac{1}{2} [A,B] \forall A, B, W: [A,[B,W]] + [B,[W,A]] + [W,[A,B]] = 0 R(X,Y)Z = \frac{1}{4} [[X,Y],Z] prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3
X = tensor(1, 0) Y = tensor(1, 0) g = metric R = curvature \forall A, B, W: R(A,B)W = \frac{1}{4} [[A,B],W] \forall A, B, C: g([A,B], C) + g(B, [A,C]) = 0 g(R(X,Y)X, Y) = \frac{1}{4} g([X,Y], [X,Y]) prove(eq3, eq1, eq2) proved from eq1, eq2

(v): for orthonormal X, Y, K(σ) = g(R(X, Y)X, Y) — in Kovalev's convention — is ¼|[X, Y]|².

27. dω on vectors partly

Problem 27 · dω(X, Y) for a 1-form

J. Ross, Part III: Differential Geometry, Example Sheet 2 (Cambridge, 2016, version 3 of 16/11/2016), question 9, p. 6 of the PDF — dec41.user.srcf.net.

9. (Differential Form Identity) Prove the identity

dω(X, Y) = Xω(Y) − Yω(X) − ω([X, Y]),

for a 1-form ω and vector fields X, Y. *Can you generalize this result to the case when ω is a p-form?

In Sucuri, partly: the formula for 1-forms, with dωμν = ∂μων − ∂νωμ (the normalization is the formula's); the generalization to p-forms, the starred part, no.

In components, the formula is Leibniz and the bracket. The normalization of d is the one the formula asks for: with another book's factor, another factor comes out — and that is why, index-free, Sucuri takes it as a hypothesis (section V).

\mu, \nu, \lambda = indices X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \omega = tensor(0, 1) X^\mu Y^\nu (\partial_\mu \omega_\nu - \partial_\nu \omega_\mu) = X^\mu \partial_\mu (\omega_\nu Y^\nu) - Y^\nu \partial_\nu (\omega_\mu X^\mu) - \omega_\mu (X^\nu \partial_\nu Y^\mu - Y^\nu \partial_\nu X^\mu) simplify(eq1) True

53. The bracket via the connection solves

Problem 53 · Commutator of vector fields via the covariant derivative

Y. Ali-Haïmoud, General Relativity (NYU, Fall 2019), Homework 6, Exercise 1, part (iii), p. 1 — cosmo.nyu.edu.

Given a vector field X, for any smooth function f : 𝓜 → ℝ, we can define the smooth function X(f) which, at each point p ∈ 𝓜, associates X|p(f). Given two vector fields X and Y, we define their commutator [X, Y] such that

[X, Y](f) ≡ X(Y(f)) − Y(X(f)).   (1)

(iii) Recall the notation convention ∇XYα ≡ Xβ∇βYα. Prove that [X, Y]α = ∇XYα − ∇YXα.

In Sucuri: all of it: the identity.

Ali-Haïmoud (NYU), HW6, 1(iii): with Levi-Civita, [X,Y] = ∇XY − ∇YX — zero torsion, read backwards.

X = tensor(1, 0) Y = tensor(1, 0) \forall A, B: \nabla_A B - \nabla_B A = [A,B] [X, Y] = \nabla_X Y - \nabla_Y X prove(eq2, eq1) eq1[A→X, B→Y]

58. The Lie derivative of a vector and of a (0, 2) solves

Problem 58 · Lie derivative of vectors, covectors and (0,2) tensors; Killing condition

J. M. Evans, Part II General Relativity, Example Sheet 3 (Cambridge, Lent 2026), question 9, p. 2 — damtp.cam.ac.uk.

9. The Lie derivative (ℒξV)α of a vector field Vα with respect to a vector field ξα (assumed to be timelike) is defined by the following conditions: (i) If {xα} is a coordinate system in which ξα = (1, 0, 0, 0), then

(ℒξV)α = ∂Vα/∂x0 = ξβ ∂Vα/∂xβ,

and (ii) (ℒξV)α transforms as a vector. Show that, in a general coordinate system,

(ℒξV)α = ξβ∇βVα − Vβ∇βξα.

Suppose, in addition, that the Lie derivative ℒξφ of a scalar field φ with respect to a vector field ξα is defined in a general coordinate system {xα} by

ℒξφ = ξα ∂φ/∂xα

and that the Lie derivative obeys the usual Leibniz rule when applied to a tensor product. Find the Lie derivative (ℒξU)α of a covector field Uα.

Write down an expression for the Lie derivative with respect to ξα of a (0, 2) tensor Tαβ and show that the condition for ξα to be a Killing vector field (as in question 6 above, i.e. ξα;β + ξβ;α = 0) is (ℒξg)αβ = 0, where gαβ is the metric tensor.

In Sucuri: all of it: ℒξV from adapted coordinates, the covariant form (torsion-free), the covector and the (0, 2) — both as in problem 17 — and the Killing equation as ℒξg = 0.

x = coordinates(u, v, w, s) xi = field(1, 0, 0, 0) V = field() bracket(xi, V) (Derivative(V^u(u, v, w, s), u))*∂_u

in adapted coordinates, ξ = ∂u: [ξ, V] = ∂uV, which is the definition; and in any chart [ξ, V]α = ξβ∂βVα − Vβ∂βξα.

\alpha, \beta, \lambda = indices V = tensor(1, 0) \xi = tensor(1, 0) \nabla = levi-civita \nabla_\alpha V^\beta = \partial_\alpha V^\beta + \Gamma^\beta{}_{\alpha\lambda} V^\lambda \Gamma = christoffel(eq1) \xi^\beta \partial_\beta V^\alpha - V^\beta \partial_\beta \xi^\alpha = \xi^\beta \nabla_\beta V^\alpha - V^\beta \nabla_\beta \xi^\alpha expand(eq2) True

the covariant form; without \nabla = levi-civita, torsion is left over.

The covector and the (0, 2) are problem 17: (ℒξω)α = ξβ∇βωα + ωβ∇αξβ, and (ℒξg)αβ = ∇αξβ + ∇βξα — the Killing equation of question 6 is ℒξg = 0.

IV. Differential forms

28. dα from the connection partly

Problem 28 · Exterior derivative via a torsion-free connection

A. G. Kovalev, Part III: Riemannian Geometry, Example Sheet 3 (Cambridge, Lent 2017), question 3 — dec41.user.srcf.net.

Γ(T*M ⊗ Λ^p T*M) --alt--> Γ(Λ^{p+1} T*M), with d : Γ(Λ^p T*M) → Γ(Λ^{p+1} T*M), where alt(ξ ⊗ α) = ξ ∧ α denotes projection to the subspace of anti-symmetric tensors (p > 0). Deduce the formula for the exterior derivative of one-forms dα(X, Y) = (∇_X α)(Y) − (∇_Y α)(X) as stated in the Lectures. ∗ Show that these results hold for any torsion-free connection ∇ on M." -->

3. Show that for the Levi-Civita connection, the following diagram commutes:

Γ(ΛpT*M) —∇→ Γ(T*M ⊗ ΛpT*M) —alt→ Γ(Λp+1T*M),  with  d : Γ(ΛpT*M) → Γ(Λp+1T*M),

where alt(ξ ⊗ α) = ξ ∧ α denotes projection to the subspace of anti-symmetric tensors (p > 0). Deduce the formula for the exterior derivative of one-forms

dα(X, Y) = (∇Xα)(Y) − (∇Yα)(X),

as stated in the Lectures. (*) Show that these results hold for any torsion-free connection ∇ on M.

In Sucuri, partly: the formula for dα for any torsion-free connection, with the formula for dα on vectors as a hypothesis; the diagram for p-forms, no.

Kovalev, sheet 3, Q3: for a torsion-free connection, dα(X,Y) = (∇Xα)(Y) − (∇Yα)(X). With (∇Xα)(Y) = X(α(Y)) − α(∇XY), and the formula for dα on vectors as a hypothesis (its factor changes with the normalization — see 27):

\alpha = form(1) X = tensor(1, 0) Y = tensor(1, 0) \forall A, B: \iota_B \iota_A \mathrm{d}\alpha = \nabla_A (\alpha(B)) - \nabla_B (\alpha(A)) - \alpha([A,B]) \forall A, B: \nabla_A B - \nabla_B A = [A,B] \iota_Y \iota_X \mathrm{d}\alpha = \nabla_X (\alpha(Y)) - \alpha(\nabla_X Y) - \nabla_Y (\alpha(X)) + \alpha(\nabla_Y X) prove(eq3, eq1, eq2) alpha(eq2[A→X, B→Y])

zero torsion goes inside α: a relation between vectors carried to one between scalars.

29–33. The algebra of forms solves (29, 30, 32) partly (31, 33)

Problem 29 · Contraction (interior product) as an antiderivation

J. Ross, Part III: Differential Geometry, Example Sheet 2 (Cambridge, 2016, version 2 of 14/10/2016), question 8 — dec41.user.srcf.net.

8. (Contractions) Let V be a vector space, X ∈ V and ω ∈ ΛkV*. Define iX(ω) by

iX(ω)(Y1, …, Yk−1) = ω(X, Y1, …, Yk−1).

Prove that iX(ω) ∈ Λk−1V*. Prove also that if ω ∈ ΛkV* and η ∈ ΛlV* then

iX(ω ∧ η) = iXω ∧ η + (−1)k ω ∧ (iXη).

Now on a smooth manifold M suppose X ∈ Vect(M) and ω ∈ Ωp(M). Show how the formula

iXω|p := iXp(ωp)

gives an element of Ωp−1(M).

In Sucuri: all of it: the antiderivation rule for ιX; the rest of the problem statement is definition.

Problem 30 · Lie derivative of a 1-form and Cartan's formula

D. Tong, General Relativity: Example Sheet 1 (2019), question 5*, p. 2 — davidtong.org.

Use the Leibniz rule to derive the formula for the Lie derivative of a 1-form ω, valid in any coordinate basis:

(LXω)μ = Xν∂νωμ + ων∂μXν

[Hint: consider (LXω)(Y) for a vector field Y.] Show that the Lie derivative of a (0, 2) tensor g is

(LXg)μν = Xρ∂ρgμν + gμρ∂νXρ + gρν∂μXρ

For a p-form η, define ιXη to be the (p − 1)-form that results from contracting a vector field X with the first index of η. Show that for a 1-form ω,

LXω = ιX(dω) + d(ιXω).

In Sucuri: all of it: both formulas in coordinates, via the Leibniz rule, and Cartan's formula on a 1-form, from the Leibniz definition of ℒX and that of dω.

Problem 31 · Properties of the exterior derivative

H. S. Reall, Part 3 General Relativity (2022), §8.1 “Introduction” (ch. 8, differential forms), “Exercises (examples sheet 4)”, eqs. (8.7)–(8.9), p. 100 — damtp.cam.ac.uk.

Context: the exterior derivative of a p-form X is the (p + 1)-form dX defined in a coordinate basis by (dX)μ1…μp+1 = (p + 1) ∂[μ1Xμ2…μp+1] (8.5).

Exercises (examples sheet 4). Show that the exterior derivative enjoys the following properties:

d(dX) = 0 (8.7)

d(X ∧ Y) = (dX) ∧ Y + (−1)p X ∧ dY (8.8)

(where Y is a q-form) and

d(φ*X) = φ* dX (8.9)

(where φ : N → M), i.e. the exterior derivative commutes with pull-back.

In Sucuri, partly: d(dX) = 0 and the graded Leibniz rule; the pull-back, no.

Problem 32 · Graded commutativity and associativity of the wedge product

H. S. Reall, Part III General Relativity (2022), §8.1 “Introduction”, exercise with eqs. (8.2)–(8.3), p. 99 — damtp.cam.ac.uk.

Definition. The wedge product of a p-form X and a q-form Y is the (p + q)-form X ∧ Y defined by

(X ∧ Y)a₁…apb₁…bq = [(p + q)!/(p! q!)] X[a₁…apYb₁…bq].   (8.1)

Exercise. Show that

X ∧ Y = (−1)pq Y ∧ X   (8.2)

(so X ∧ X = 0 if p is odd); and

(X ∧ Y) ∧ Z = X ∧ (Y ∧ Z),   (8.3)

i.e. the wedge product is associative so we don't need to include the brackets.

In Sucuri: all of it: (8.2) and (8.3).

Problem 33 · Hodge dual: ⋆⋆ and ⋆d⋆

H. S. Reall, Part 3 General Relativity (2022), §8.5 “Volume form”, Lemma with eqs. (8.56)–(8.57) (“Proof. Exercise (use (8.54))”), pp. 107–108 — damtp.cam.ac.uk.

Context: ε is the volume form of an oriented n-dimensional manifold with metric, and

εa1…apcp+1…cn εb1…bpcp+1…cn = ±p!(n − p)! δa1[b1 … δapbp] (8.54),

where the upper (lower) sign holds for Riemannian (Lorentzian) signature.

Definition. On an oriented manifold with metric, the Hodge dual of a p-form X is the (n − p)-form ⋆X defined by

(⋆X)a1…an−p = (1/p!) εa1…an−pb1…bp Xb1…bp (8.55)

Lemma. For a p-form X,

⋆(⋆X) = ±(−1)p(n−p) X (8.56)

(⋆d⋆X)a1…ap−1 = ±(−1)p(n−p) ∇bXa1…ap−1b (8.57)

where the upper (lower) sign holds for Riemannian (Lorentzian) signature.

Proof. Exercise (use (8.54)).

In Sucuri, partly: ⋆⋆ = ±(−1)p(n−p), with the declared signature; ⋆d⋆ in terms of ∇, no.

Cambridge, Part III Differential Geometry, sheet 2, Q8 (ιX is an antiderivation); Tong, sheet 1, Q5–6 (Cartan on a 1-form; d² = 0 and graded Leibniz); Reall §8.1 (α∧β = (−1)pqβ∧α, associativity) and §8.5 (⋆⋆ = ±(−1)p(n−p)).

\omega = form(1) \eta = form(2) X = tensor(1, 0) \iota_X (\omega \wedge \eta) = \iota_X \omega \, \eta - \omega \wedge \iota_X \eta prove(eq1) proved from no hypothesis
\omega = form(1) \eta = form(2) \mathrm{d}(\omega \wedge \eta) = \mathrm{d}\omega \wedge \eta - \omega \wedge \mathrm{d}\eta prove(eq1) proved from no hypothesis

likewise d²η = 0, ω∧η = η∧ω (degree 1 with degree 2), associativity. Cartan's formula also goes through, but be clear: in Sucuri it is the definition of ℒX on forms, and the exercise becomes an identity.

\alpha = form(1) F = form(2) F = \mathrm{d}\alpha \mathrm{d} F = 0 prove(eq2, eq1) d(eq1)

Maxwell: F = dA implies dF = 0 — the step is applying d to both sides of the hypothesis.

g = metric(-,+,+,+) \star = hodge F = form(2) \star \star F = -F prove(eq1) proved from no hypothesis

in 4D Lorentzian, ⋆⋆ = −1 on 2-forms and +1 on 1- and 3-forms, as Reall says. The expression of ⋆d⋆ via ∇ (the second half of the exercise) does not come out.

30, by a route other than the forms engine — where ℒX is defined by Cartan's formula, and the formula becomes an identity. With ℒX from the Leibniz rule, (ℒXω)(Y) = X(ω(Y)) − ω([X, Y]), the components come out by expanding, and Cartan's formula is a consequence:

\mu, \nu, \rho = indices X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \omega = tensor(0, 1) h = tensor(0, 2, symmetric) X^\nu \partial_\nu (\omega_\mu Y^\mu) - \omega_\mu (X^\nu \partial_\nu Y^\mu - Y^\nu \partial_\nu X^\mu) = (X^\nu \partial_\nu \omega_\mu + \omega_\nu \partial_\mu X^\nu) Y^\mu simplify(eq1) True

(ℒXω)μ = Xν∂νωμ + ων∂μXν.

\mu, \nu, \rho = indices X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \omega = tensor(0, 1) h = tensor(0, 2, symmetric) X^\rho \partial_\rho (h_{\mu\nu} Y^\mu Z^\nu) - h_{\mu\nu} (X^\rho \partial_\rho Y^\mu - Y^\rho \partial_\rho X^\mu) Z^\nu - h_{\mu\nu} Y^\mu (X^\rho \partial_\rho Z^\nu - Z^\rho \partial_\rho X^\nu) = (X^\rho \partial_\rho h_{\mu\nu} + h_{\mu\rho} \partial_\nu X^\rho + h_{\rho\nu} \partial_\mu X^\rho) Y^\mu Z^\nu simplify(eq1) True

(ℒXg)μν, with g a symmetric (0, 2).

X = tensor(1, 0) Y = tensor(1, 0) \omega = tensor(0, 1) L = tensor(0, 1) D = tensor(0, 2) \forall B: L(B) = \nabla_X \omega(B) - \omega([X,B]) \forall A, B: D(A,B) = \nabla_A \omega(B) - \nabla_B \omega(A) - \omega([A,B]) L(Y) = D(X,Y) + \nabla_Y \omega(X) prove(eq3, eq1, eq2) proved from eq1, eq2

L = ℒXω and D = dω: (ℒXω)(Y) = (ιXdω)(Y) + (dιXω)(Y), Cartan's formula.

47, 48 and 51. d, ι and ℒ solves (48, 51) partly (47)

Problem 47 · d commutes with the Lie derivative; Leibniz rule for LX

C. Wendl, Differential Geometry I, Problem Set 7 (Humboldt-Universität zu Berlin, 2016–17), problem 1(a) — mathematik.hu-berlin.de.

1. The goal of this problem is to prove Cartan's formula for the Lie derivative of a differential form,

LXω = dιXω + ιXdω.   (1)

Recall that for a k-form ω ∈ Ωk(M) on a smooth n-manifold M, LXω ∈ Ωk(M) is defined as ∂t φt*ω|t=0 where φt : M → M denotes the time t flow of the vector field X ∈ Vec(M).

(a) Use the definition of the Lie derivative to prove the relation

d(LXω) = LX(dω)

and the Leibniz rule

LX(α ∧ β) = LXα ∧ β + α ∧ LXβ.   (2)

Observe that this formula determines the action of LX on arbitrary differential forms if we know how it acts on 0-forms (i.e. smooth functions) and exact 1-forms (i.e. differentials of smooth functions).

In Sucuri, partly: d ℒX = ℒX d and the Leibniz rule for ℒX (the example under 48). But the exercise asks for both from the flow definition of ℒX, in order to then prove Cartan's formula; in Sucuri, ℒX on forms is defined by Cartan's formula, and that route cannot be walked.

Problem 48 · dιX + ιXd satisfies the Leibniz rule (Cartan's formula)

C. Wendl, Differential Geometry I, Problem Set 7 (Humboldt-Universität zu Berlin, 2016–17), problem 1(d) — mathematik.hu-berlin.de.

Context: problem 1 aims to prove Cartan's formula LXω = dιXω + ιXdω (1). Part (a) proves the Leibniz rule LX(α ∧ β) = LXα ∧ β + α ∧ LXβ (2). Part (c) proves that ιv satisfies the graded Leibniz rule ιv(α ∧ β) = ιvα ∧ β + (−1)|α| α ∧ ιvβ.

(d) For any fixed X ∈ Vec(M), use the graded Leibniz rules satisfied by d and ιX to show that the operator (d ∘ ιX + ιX ∘ d) : Ωk(M) → Ωk(M) also satisfies the Leibniz rule (2). Deduce that this operator matches LX.

Author's note (footnote 2): in contrast to the exterior derivative and the interior product, the Leibniz rule satisfied by LX does not include any annoying signs. This is consistent with thinking of LX : Ωk(M) → Ωk(M) as an object of degree zero (hence even), of d as an object of degree one (hence odd), and of ιX : Ωk(M) → Ωk−1(M) as an object of degree −1 (hence odd).

In Sucuri: all of it: that dιX + ιXd satisfies the Leibniz rule.

Problem 51 · d² = 0 on ℝ³: curl grad = 0 and div curl = 0

C. Wendl, Differential Geometry I, Problem Set 6 (Humboldt-Universität zu Berlin, 2016–17), problem 6(c) — mathematik.hu-berlin.de.

Context (preamble of problem 6 and part (b)): given a volume form μ ∈ Ωn(M) on an n-manifold M, the divergence of X ∈ Vec(M) is the unique function div(X) : M → ℝ such that LXμ = div(X) μ. On M = ℝ³, with μ = dx ∧ dy ∧ dz and X = Xx∂x + Xy∂y + Xz∂z, one has div(X) = ∂xXx + ∂yXy + ∂zXz, an expression sometimes also denoted by ∇ · X.

(c) Recall that on ℝ³, the gradient of a function f : ℝ³ → ℝ is the vector field

grad(f) = ∇f := (∂xf)∂x + (∂yf)∂y + (∂zf)∂z,

and the curl of a vector field X = Xx∂x + Xy∂y + Xz∂z is the vector field

curl(X) = ∇ × X := (∂yXz − ∂zXy)∂x + (∂zXx − ∂xXz)∂y + (∂xXy − ∂yXx)∂z.

Using the relations of these operations to differential forms and the exterior derivative, deduce from d² = 0 the formulas

∇ × (∇f) = 0   and   ∇ · (∇ × X) = 0

for all f ∈ C∞(ℝ³) and X ∈ Vec(ℝ³).

In Sucuri: all of it: curl grad = 0 and div curl = 0, with forms and with indices.

Wendl, problem set 7, 1(a) and 1(d): [d, ℒX] = 0; and P = dιX + ιXd is a degree-0 derivation — the step that proves Cartan's formula in any degree. Wendl, problem set 6, 6(c): d² = 0, which on ℝ³ is curl grad = 0 and div curl = 0.

\omega = form(1) \eta = form(2) X = tensor(1, 0) \mathrm{d} \mathcal{L}_X \eta = \mathcal{L}_X \mathrm{d} \eta prove(eq1) proved from no hypothesis
\omega = form(1) \eta = form(2) X = tensor(1, 0) \mathrm{d} \iota_X (\omega \wedge \eta) + \iota_X \mathrm{d} (\omega \wedge \eta) = (\mathrm{d} \iota_X \omega + \iota_X \mathrm{d} \omega) \wedge \eta + \omega \wedge (\mathrm{d} \iota_X \eta + \iota_X \mathrm{d} \eta) prove(eq1) proved from no hypothesis

the signs of d and of ι cancel in pairs: what is left is a derivation with no sign, like that of ℒX.

46. ℒ[X,Y] = [ℒX, ℒY] partly

Problem 46 · ℒ[X,Y] = [ℒX, ℒY] on forms

C. Wendl, Differential Geometry I, Problem Set 7 (Humboldt-Universität zu Berlin, 2016–17), problem 4(a), p. 3 — mathematik.hu-berlin.de.

4. (a) Use Leibniz rules as in Problem 1 to show that for all X, Y ∈ Vec(M) and ω ∈ Ωk(M),

ℒ[X,Y]ω = ℒXℒYω − ℒYℒXω.

In Sucuri, partly: on functions — it is the definition of the bracket — and on 1-forms, by Leibniz and Jacobi; on higher-degree forms, the argument that two derivations agreeing on functions and on df agree everywhere, no.

On a 1-form, with ℒX from the Leibniz rule, (ℒXω)(B) = X(ω(B)) − ω([X, B]): P = ℒYω, Q = ℒXω, U = ℒXP, W = ℒYQ and E = ℒ[X,Y]ω.

X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \omega = tensor(0, 1) P = tensor(0, 1) Q = tensor(0, 1) U = tensor(0, 1) W = tensor(0, 1) E = tensor(0, 1) \forall B: P(B) = \nabla_Y \omega(B) - \omega([Y,B]) \forall B: Q(B) = \nabla_X \omega(B) - \omega([X,B]) \forall B: U(B) = \nabla_X P(B) - P([X,B]) \forall B: W(B) = \nabla_Y Q(B) - Q([Y,B]) \forall B: E(B) = \nabla_{[X,Y]} \omega(B) - \omega([[X,Y],B]) \forall A, B, C: [A,[B,C]] + [B,[C,A]] + [C,[A,B]] = 0 U(Z) - W(Z) = E(Z) prove(eq7, eq1, eq2, eq3, eq4, eq5, eq6) proved from eq1, eq2, eq3, eq4, eq5, eq6

(ℒXℒY − ℒYℒX)ω = ℒ[X,Y]ω — and it is also part (i) of problem 44.

49. dω(X, Y, Z) of a 2-form solves

Problem 49 · dω(X, Y, Z) for a 2-form

C. Wendl, Differential Geometry I, Problem Set 6 (Humboldt-Universität zu Berlin, 2016–17), problem 4 (preamble) and part (c), pp. 2–3 — mathematik.hu-berlin.de.

4. Recall that if (x1, …, xn) : U → ℝn is a chart defined on an open subset in some n-manifold M, any k-form ω ∈ Ωk(M) can be written on U as

ω = ωi1…ik dxi1 ⊗ … ⊗ dxik = (1/k!) ωi1…ik dxi1 ∧ … ∧ dxik = ∑i1<…<ik ωi1…ik dxi1 ∧ … ∧ dxik,

where the first two expressions use the Einstein summation convention and the third one does not. Here the component functions ωi1…ik : U → ℝ can be written in terms of the coordinate vector fields ∂1, …, ∂n as ωi1…ik = ω(∂i1, …, ∂ik). In order to write down a coordinate formula for the exterior derivative, we introduce the following notation: given any collection of functions Ti1…ik on U labeled by the indices i1, …, ik, define

T[i1…ik] := (1/k!) ∑σ∈Sk (−1)|σ| Tiσ(1)…iσ(k),

so for instance if Ti1…ik are the components of a tensor field T, then Alt(T)i1…ik = T[i1…ik], and the wedge product of α ∈ Ωk(M) and β ∈ Ωℓ(M) can now be written in coordinates as

(α ∧ β)i1…ikj1…jℓ = α[i1…ikβj1…jℓ].

[…]

(c) It now follows from Problem Set 4 #1(a) that the exterior derivative of a 1-form λ can also be written as

dλ(X, Y) = ℒX(λ(Y)) − ℒY(λ(X)) − λ([X, Y]).

Indeed, the right hand side is C∞-linear with respect to vector fields X, Y ∈ Vec(M) and thus defines a tensor field, whose component functions we've seen match the formula from part (b). Prove the corresponding formula for the exterior derivative of a 2-form,

dω(X, Y, Z) = ℒX(ω(Y, Z)) + ℒY(ω(Z, X)) + ℒZ(ω(X, Y)) − ω([X, Y], Z) − ω([Y, Z], X) − ω([Z, X], Y).

Remark: Similar formulas exist for the exterior derivatives of k-forms for all k > 2, though I cannot recall ever having needed to use them.

In Sucuri: all of it: the formula, in components, with (dω)λμν = ∂λωμν + ∂μωνλ + ∂νωλμ, its normalization.

\mu, \nu, \lambda = indices X = tensor(1, 0) Y = tensor(1, 0) Z = tensor(1, 0) \omega = tensor(0, 2, antisymmetric) X^\lambda Y^\mu Z^\nu (\partial_\lambda \omega_{\mu\nu} + \partial_\mu \omega_{\nu\lambda} + \partial_\nu \omega_{\lambda\mu}) = X^\lambda \partial_\lambda (\omega_{\mu\nu} Y^\mu Z^\nu) + Y^\lambda \partial_\lambda (\omega_{\mu\nu} Z^\mu X^\nu) + Z^\lambda \partial_\lambda (\omega_{\mu\nu} X^\mu Y^\nu) - \omega_{\mu\nu} (X^\lambda \partial_\lambda Y^\mu - Y^\lambda \partial_\lambda X^\mu) Z^\nu - \omega_{\mu\nu} (Y^\lambda \partial_\lambda Z^\mu - Z^\lambda \partial_\lambda Y^\mu) X^\nu - \omega_{\mu\nu} (Z^\lambda \partial_\lambda X^\mu - X^\lambda \partial_\lambda Z^\mu) Y^\nu simplify(eq1) True

the right-hand side is Wendl's formula: ℒX(ω(Y, Z)) + ℒY(ω(Z, X)) + ℒZ(ω(X, Y)) − ω([X, Y], Z) − ω([Y, Z], X) − ω([Z, X], Y).

34. Maurer–Cartan equations partly

Problem 34 · Maurer–Cartan equations

J. Ross, Part III: Differential Geometry, Example Sheet 2 (Cambridge, 2016, version 3 of 16/11/2016), question 11, p. 7 of the PDF — dec41.user.srcf.net.

11. Let G be Lie group that is a subgroup of GLn(ℝ) and Xi, i = 1, …, d = dim G, be linearly independent left-invariant vector fields on G induced by a basis of TIG. Show that the condition that ωi(Xj) = δij identically on G defines a system of pointwise linearly independent smooth 1-forms ωi on G. Show further that the 1-forms ωi are left-invariant in the sense that

Lg*(ωi) = ωi,   for every g ∈ G.

Let Ckij be a set of real constants determined by [Xi, Xj] = ∑k Ckij Xk. Deduce from the identity of the previous question the formula

dωk = −½ ∑i,j Ckij ωi ∧ ωj.

In Sucuri, partly: the formula dωk = −½ Ckij ωi∧ωj, deduced from the identity of the previous question — first without indices, dω(X, Y) = −ω([X, Y]) for constant ω(X), ω(Y); then with indices, summing over the basis. That the ωi exist and are smooth, independent and left-invariant is argument, and does not count.

X = tensor(1, 0) Y = tensor(1, 0) \omega = tensor(0, 1) D = tensor(0, 2, antisymmetric) \forall A, B: D(A, B) = \nabla_A \omega(B) - \nabla_B \omega(A) - \omega([A, B]) \forall A: \nabla_A \omega(X) = 0 \forall A: \nabla_A \omega(Y) = 0 D(X, Y) = -\omega([X, Y]) prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

D is dω, by the identity of the previous question; X, Y are two of the Xi, and ω one of the ωk: ωk(Xj) = δkj is constant.

Now the sum over the basis. With indices, the basis Xi is Xai and the dual is ωia: the index i numbers the fields, and runs from 1 to d = dim G like the manifold index a — which is why the two can be of the same type. Dkab are the components of dωk, Bcij those of [Xi, Xj], and δ is the declared delta. The two dualities, ω(X) = δ and X ω = δ, are hypotheses: the second says the Xi form a basis.

a, b, c, i, j, k, l = indices \delta = kronecker D = tensor(1, 2) B = tensor(1, 2) C = tensor(1, 2) \omega = tensor(1, 1) X = tensor(1, 1) D^k{}_{ab} X^a{}_i X^b{}_j = -\omega^k{}_c B^c{}_{ij} B^c{}_{ij} = C^l{}_{ij} X^c{}_l \omega^i{}_a X^a{}_j = \delta^i_j D^k{}_{ab} X^a{}_i X^b{}_j = -C^k{}_{ij} prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

eq1 is what was proved above, on the Xi; eq2, the definition of C: dωk(Xi, Xj) = −Ckij.

a, b, c, i, j, k, l = indices \delta = kronecker D = tensor(1, 2) B = tensor(1, 2) C = tensor(1, 2) \omega = tensor(1, 1) X = tensor(1, 1) D^k{}_{ab} X^a{}_i X^b{}_j = -C^k{}_{ij} X^a{}_i \omega^i{}_b = \delta^a_b C^k{}_{ij} + C^k{}_{ji} = 0 D^k{}_{ab} = -\frac{1}{2} C^k{}_{ij} (\omega^i{}_a \omega^j{}_b - \omega^i{}_b \omega^j{}_a) prove(eq4, eq1, eq2, eq3) proved from eq1, eq2, eq3

(ωi∧ωj)ab = ωiaωjb − ωibωja: this is dωk = −½ Ckij ωi∧ωj.

Forms in a chart

So far the forms were abstract. In a chart, a form is written with the d's of the coordinates — \alpha = form(a\,dr + b\,d\theta + c\,d\phi), H = form(H_x dy \wedge dz + …) — and the operations are computations: wedge, exterior (the d), star (the ⋆, with the metric), interior and lie (with a field), and equal. An operation with a name on the left stores the result: S = star(\beta, g). The orientation is that of the coordinate order: the volume form is √|g| dx¹∧…∧dxⁿ, and swapping two coordinates flips the sign of ⋆.

50. The divergence from the Lie derivative solves

Problem 50 · The divergence via ℒXμ = div(X) μ

C. Wendl, Differential Geometry I, Problem Set 6 (Humboldt-Universität zu Berlin, 2016–17), problem 6 (preamble) and part (b), p. 3 — mathematik.hu-berlin.de.

6. Given a volume form μ ∈ Ωn(M) on an n-manifold M, one can define volumes of compact regions U ⊂ M by

Vol(U) := ∫U μ.

The divergence of a vector field X ∈ Vec(M) can then be defined in terms of the Lie derivative of μ with respect to X: let div(X) : M → ℝ be the unique real-valued function such that

ℒXμ = div(X) μ.

Note that this is well defined since ℒXμ is an n-form and the space ΛnTp*M of n-forms at each point p ∈ M is 1-dimensional. Observe also that dμ = 0 since Ωn+1(M) = {0}, so Cartan's formula implies div(X) μ = dιXμ, which matches the formula we saw in lecture for the case M = ℝn.

[…]

(b) Show that in the case M = ℝ3 with μ = dx ∧ dy ∧ dz and X = Xx∂x + Xy∂y + Xz∂z ∈ Vec(ℝ3) in standard Cartesian coordinates (x, y, z),

div(X) = ∂xXx + ∂yXy + ∂zXz.

The latter expression is sometimes also denoted by ∇ · X.
Note: One can show more generally that if μ = dx1 ∧ … ∧ dxn on ℝn then div(X) = ∂iXi.

In Sucuri: all of it: div(X) = ∂xXx + ∂yXy + ∂zXz, from ℒXμ.

x = coordinates(x, y, z) X = field() \mu = form(dx \wedge dy \wedge dz) lie(X, \mu) (Derivative(X^x(x, y, z), x) + Derivative(X^y(x, y, z), y) + Derivative(X^z(x, y, z), z))*dx∧dy∧dz

ℒXμ = dιXμ + ιXdμ, with a generic X: the coefficient of μ is the divergence.

52. Contact forms partly

Problem 52 · Contact forms on 3-manifolds

C. Wendl, Differential Geometry I, Problem Set 12 (Humboldt-Universität zu Berlin, 2016–17), problem 4, parts (a)–(d), pp. 3–4 — mathematik.hu-berlin.de.

4. Suppose M is an oriented 3-manifold and λ ∈ Ω1(M) is nowhere zero, i.e. for all p ∈ M there exist vectors X ∈ TpM with λ(X) ≠ 0. Then at every p ∈ M, the kernel ker λp = {X ∈ TpM | λ(X) = 0} is a 2-dimensional subspace of TpM, and the union of these for all p defines a smooth 2-dimensional distribution

ξ := ker λ ⊂ TM.

(a) Show that the following conditions are equivalent:
i. λ ∧ dλ ≡ 0;
ii. For all p ∈ M and X, Y ∈ TpM, dλ(X, Y) = 0 [sic; from context, read X, Y ∈ ξp];
iii. ξ is integrable.
Hint: In Problem Set 6 #4(c), you will find a useful formula for dλ as a C∞-bilinear form on vector fields. Combine this with the Frobenius theorem.

The 1-form λ is called a contact form if λ ∧ dλ is a volume form; the distribution ξ (called the contact structure) is then “as non-integrable as possible.” An example on ℝ3 is shown in the figure in the original. Such examples can be constructed by the following trick. Let (ρ, φ, z) denote the standard cylindrical coordinates on ℝ3, so x = ρ cos φ and y = ρ sin φ. Choose smooth real-valued functions f(ρ), g(ρ) and define λ at (ρ, φ, z) by

λ = f(ρ) dz + g(ρ) dφ.   (2)

(b) Since the coordinates (ρ, φ, z) are not well defined at ρ = 0, there is of course some danger that the 1-form defined in Equation (2) might be singular at the z-axis. Show that λ is in fact smooth on all of ℝ3 and satisfies λ ∧ dλ ≠ 0 near the z-axis if we assume f(ρ) = 1 and g(ρ) = ρ2 for ρ sufficiently close to 0. Hint: Convert to Cartesian coordinates.

(c) Assuming f and g take the form described above for ρ near 0, show that λ is a contact form if and only if

f(ρ)g′(ρ) − f′(ρ)g(ρ) ≠ 0

for all ρ > 0. What does this mean geometrically about the curve ρ ↦ (f(ρ), g(ρ)) ∈ ℝ2? Interpret this in terms of “twisting” of the planes ξp as p ∈ ℝ3 moves along radial paths away from the z-axis.

(d) By a fundamental result in contact geometry known as Gray's theorem, contact structures have the following remarkable “stability” property: if M is a closed 3-manifold and {ξt}t∈[0,1] is any smooth family of contact structures on M, then they are all “equivalent” in the sense that there exists a smooth family of diffeomorphisms {ϕt : M → M}t∈[0,1] such that ϕ0 = Id and Tϕt(ξ0) = ξt for all t ∈ [0, 1]. Show that this is not true in general for arbitrary smooth families of distributions, e.g. it becomes false if we assume that ξ0 is integrable but ξt is a contact structure for each t > 0.

In Sucuri, partly: (c): λ∧dλ = (f g′ − f′ g) dρ∧dφ∧dz, and λ is contact if and only if that coefficient does not vanish; the equivalences in (a), the regularity on the axis in (b) and the stability in (d), no.

x = coordinates(\rho, \phi, z) f = f(\rho) h = h(\rho) \lambda = form(f dz + h d\phi) D = exterior(\lambda) wedge(\lambda, D) (f(rho)*Derivative(h(rho), rho) - h(rho)*Derivative(f(rho), rho))*drho∧dphi∧dz

with Wendl's g renamed h, because g is the metric.

54. ⋆ in spherical coordinates solves

Problem 54 · Hodge dual in spherical coordinates

T. Dray, MTH 434/534, HW #3 (Oregon State University, winter 2024), problem 1 — sites.science.oregonstate.edu.

1. Hodge dual in spherical coordinates. Consider spherical coordinates in 3-dimensional Euclidean space with the usual orientation, namely ω = r2 sin θ dr ∧ dθ ∧ dφ.
Warning: these are “physics” conventions: θ is the angle from the north pole (colatitude), and φ is the angle in the xy-plane (longitude).

(a) Determine the Hodge dual operator ∗ on all forms (expressed in spherical coordinates) by computing its action on basis forms at each rank.

(b) Compute the dot and cross products of two generic 1-forms in spherical coordinates using the expressions

α · β = ∗(α ∧ ∗β),
α × β = ∗(α ∧ β).

You may express your results either with respect to an orthonormal basis or with respect to a “coordinate” (non-orthonormal) spherical basis; make sure you know which you're doing. (“Generic” means for any two 1-forms, e.g. in terms of their components.)

In Sucuri: all of it: (a) ⋆ on all basis forms, and (b) α·β and α×β of two generic 1-forms — in the coordinate basis, which the problem allows, and says so.

x = coordinates(r, \theta, \phi) g = metric(1, r^2, r^2 \sin^2\theta) \alpha = form(dr) star(\alpha, g) r**2*sin(theta)*dtheta∧dphi

and ⋆dθ = sin θ dφ∧dr, ⋆dφ = dr∧dθ/sin θ, ⋆(dr∧dθ) = sin θ dφ, ⋆(dθ∧dφ) = dr/(r² sin θ), ⋆(dφ∧dr) = dθ/sin θ; ⋆1 = r² sin θ dr∧dθ∧dφ and ⋆(dr∧dθ∧dφ) = 1/(r² sin θ).

x = coordinates(r, \theta, \phi) g = metric(1, r^2, r^2 \sin^2\theta) \alpha = form(a_1 dr + a_2 d\theta + a_3 d\phi) \beta = form(b_1 dr + b_2 d\theta + b_3 d\phi) S = star(\beta, g) C = wedge(\alpha, S) star(C, g) (a_{1}*b_{1}*r**2 + a_{2}*b_{2} + a_{3}*b_{3}/sin(theta)**2)/r**2

α·β = ⋆(α∧⋆β) = gijaibj; with wedge(\alpha, \beta) instead of S, α×β = ⋆(α∧β).

55. ⋆ in Minkowski space solves

Problem 55 · Hodge dual in Minkowski space

T. Dray, MTH 434/534, HW #3 (Oregon State University, winter 2024), problem 2 — sites.science.oregonstate.edu.

2. Hodge dual in Minkowski space. 4-dimensional Minkowski space has an orthonormal, oriented basis of 1-forms given by

{dx, dy, dz, dt},

with g(dt, dt) = −1, g(dx, dx) = g(dy, dy) = g(dz, dz) = 1, and all others zero. The “volume element” (choice of orientation) is given by ω = dx ∧ dy ∧ dz ∧ dt.

(a) Determine the Hodge dual operator ∗ on all forms by computing its action on basis forms at each rank.

(b) How does your answer change if the opposite orientation is chosen, namely

ω = dt ∧ dx ∧ dy ∧ dz?

In Sucuri: all of it: (a) ⋆ on the basis forms, with ω = dx∧dy∧dz∧dt, and (b) the opposite orientation, which flips the sign of all of them.

x = coordinates(x, y, z, t) g = metric(1, 1, 1, -1) \alpha = form(dt) star(\alpha, g) dx∧dy∧dz

⋆dx = dy∧dz∧dt, ⋆dt = dx∧dy∧dz, ⋆(dx∧dt) = −dy∧dz, ⋆(dx∧dy) = dz∧dt, …

x = coordinates(t, x, y, z) g = metric(-1, 1, 1, 1) \alpha = form(dt) star(\alpha, g) -dx∧dy∧dz

(b): with the coordinates in the order (t, x, y, z), the orientation is dt∧dx∧dy∧dz = −ω, and every ⋆ flips sign.

56. Vector operators in parabolic coordinates solves

Problem 56 · Gradient, curl, divergence and Laplacian with d and ∗

T. Dray, MTH 434/534, HW #4 (Oregon State University, winter 2024), problem 1 — sites.science.oregonstate.edu.

1. Orthogonal coordinates. Choose any orthogonal coordinate system in 3-dimensional Euclidean space ℝ3 other than Cartesian, cylindrical, or spherical coordinates. (You may see me for suggestions.) Working in an orthonormal basis, compute the gradient and Laplacian of an arbitrary function, and the curl and divergence of an arbitrary “vector field” (again, really a 1-form), using the expressions

∇f = df,
∇ × α = ∗dα,
∇ · α = ∗d∗α,
△f = ∇ · ∇f = ∗d∗df.

You may check your answer in standard reference books, but you should use exterior differentiation and Hodge duality in your computation.

In Sucuri: all of it: gradient, curl, divergence and Laplacian from the four formulas, in parabolic coordinates — and in the orthonormal cobasis, as the problem asks.

Parabolic: ds² = (u² + v²)(du² + dv²) + u²v²dφ². ortonormal(α, g) writes α in the cobasis σi = √|gii| dxi.

x = coordinates(u, v, \phi) g = metric(u^2 + v^2, u^2 + v^2, u^2 v^2) f = f(u, v, \phi) F = form(f) D = exterior(F) ortonormal(D, g) Derivative(f(u, v, phi), u)/sqrt(u**2 + v**2)*σ^u

∇f = df, in the orthonormal cobasis.

x = coordinates(u, v, \phi) g = metric(u^2 + v^2, u^2 + v^2, u^2 v^2) f = f(u, v, \phi) F = form(f) D = exterior(F) S = star(D, g) T = exterior(S) star(T, g) (u**2*Derivative(f(u, v, phi), (phi, 2)) + u*v*(u*v*Derivative(f(u, v, phi), (u, 2)) + u*v*Derivative(f(u, v, phi), (v, 2)) + u*Derivative(f(u, v, phi), v) + v*Derivative(f(u, v, phi), u)) + v**2*Derivative(f(u, v, phi), (phi, 2)))/(u**2*v**2*(u**2 + v**2))

△f = ⋆d⋆df = [∂u(u∂uf)/u + ∂v(v∂vf)/v]/(u² + v²) + ∂²φf/(u²v²), the one in the tables. The curl ⋆dα and the divergence ⋆d⋆α of α = a du + b dv + c dφ come out the same way.

57. A 2-form that does not decompose solves

Problem 57 · Decomposable forms

T. Dray, MTH 434/534 – HW #2 (Oregon State, winter 2024), problem 1 “Decomposable forms”, p. 1 — sites.science.oregonstate.edu.

Denote the p-forms in ℝn by ⋀p(ℝn). A typical 1-form in ℝ2 would therefore take the form F = Fx dx + Fy dy ∈ ⋀1(ℝ2). A p-form β ∈ ⋀p(ℝn) is called decomposable if there exist 1-forms αi ∈ ⋀1(ℝn) with β = α1 ∧ … ∧ αp.

(a) Show that all elements of ⋀2(ℝ3), that is, all 2-forms in ℝ3, are decomposable. In other words, show that H = Hx dy∧dz + Hy dz∧dx + Hz dx∧dy is decomposable. Hint: consider the previous assignment! You may cite your solution to the previous assignment without proof, so long as an explicit reference is given (“see HW #1”). If you do this, it wouldn't hurt to include a copy of your previous assignment.

(b) Find an example of an indecomposable 2-form γ ∈ ⋀p(ℝn). Hint: don't work in ℝ3…

(c) Is γ ∧ γ = 0? Should it be? Can it be?

(d) MTH 534 students only (or extra credit): show that all 3-forms are decomposable in ℝ4.

Extra credit: can you argue that all elements of ⋀n−1(ℝn) are decomposable?

In Sucuri: all of it: (a) every 2-form on ℝ³ is α∧β — the decomposition, for Hx ≠ 0, checked; the other cases permute the coordinates —, (b) and (c), and (d): every 3-form on ℝ⁴ is α∧β∧γ, the same way. The extra credit, in ℝⁿ, is left to the reader.

Dray, MTH 434, HW2 #1: every 2-form on ℝ³ is α∧β, but on ℝ⁴ γ = dx∧dy + dz∧dw is not — γ∧γ = 2 dx∧dy∧dz∧dw ≠ 0, whereas (α∧β)∧(α∧β) = 0.

\alpha = form(1) \beta = form(1) \gamma = form(1) \kappa = form(1) (\alpha \wedge \beta + \gamma \wedge \kappa) \wedge (\alpha \wedge \beta + \gamma \wedge \kappa) = 2 \alpha \wedge \beta \wedge \gamma \wedge \kappa prove(eq1) proved from no hypothesis
x = coordinates(x, y, z) H = form(H_x dy \wedge dz + H_y dz \wedge dx + H_z dx \wedge dy) \alpha = form(H_x dy - H_y dx) \beta = form(\frac{1}{H_x} (H_x dz - H_z dx)) P = wedge(\alpha, \beta) equal(P, H) True

(a): H = (Hxdy − Hydx) ∧ (Hxdz − Hzdx)/Hx.

x = coordinates(x, y, z, w) H = form(v_1 dy \wedge dz \wedge dw - v_2 dx \wedge dz \wedge dw + v_3 dx \wedge dy \wedge dw - v_4 dx \wedge dy \wedge dz) \alpha = form(v_1 dy - v_2 dx) \beta = form(v_1 dz - v_3 dx) \gamma = form(\frac{1}{v_1^2} (v_1 dw - v_4 dx)) A = wedge(\alpha, \beta) P = wedge(A, \gamma) equal(P, H) True

(d): the general 3-form on ℝ⁴, with v1 ≠ 0, is a product of three 1-forms.

V. What does not come out

Saying this is part of the job. Each item below is an exercise from the problem sets, and the reason is what the program lacks — not the exercise.

Those that come out only partly because the problem asks for more than the central computation — geodesics, embeddings, volumes, parallel transport — are not listed here one by one: the note under each problem statement says what is left out.

No exercise in the list comes out with nothing. What the program lacks, in the partial ones, is this:

VI. What the exercises taught

A notation reader can only be tested with notation written by other people, for another purpose. Working through these problem sets found sixteen bugs — eight of them silent, the kind of error this program exists not to make. All fixed, each with a test.

Sources

The manual explains every declaration and every verb used here: manual. The code is at github.com/RafaelCRdeLima/SUCURI.