Tetrad and Newman–Penrose Computation

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

A tetrad calculation replaces coordinate components by projections onto an orthonormal or null frame. This can expose algebraic structure that a coordinate array hides, but it introduces frame gauge freedom: two correct computations may report different spin coefficients or Weyl scalars because their null legs differ by boosts, spins, exchanges, or orientations. Reliable symbolic work therefore declares the frame and tests its normalization before computing curvature.

Required background. Connections and Covariant Derivatives supplies frame derivatives and connection coefficients, while Curvature supplies the Weyl and Ricci projections being computed.

Helpful background. Formulations of General Relativity provides the continuum field equations behind the projections; Ordinary Differential Equations supports later transport and perturbation systems; and Spacetime Symmetries and Killing Fields helps choose symmetry-adapted frames without confusing adaptation with uniqueness.

Specifying the null frame

Use the ordered contravariant tetrad $(\ell^\mu,n^\mu,m^\mu,\bar m^\mu)$ with

\ell\!\cdot n=-1, \qquad m\!\cdot\bar m=1, \qquad \ell^2=n^2=m^2=\bar m^2=0,

and all other inner products zero. Use the site’s Weyl-scalar definition

\Psi_2=-C_{\alpha\beta\gamma\delta} \ell^\alpha m^\beta\bar m^\gamma n^\delta.

For comparison with a package defining $\Psi_2^{\rm alt}=+C(\ell,m,\bar m,n)$ , the translation is

\Psi_2^{\rm alt}=-\Psi_2.

Keep this map with the imported output; the site convention card defines all five Weyl scalars used here. These choices form a package: changing metric signature, Riemann sign, tetrad order, or the sign in the definition of $\Psi_2$ changes downstream formulas. The original null-tetrad construction and its signature-dependent normalization appear in Newman and Penrose 1962, § II, Eq. 2, pp. 566–567. The algorithm should first construct the $4\times4$ Gram matrix of tetrad inner products and compare it with the declared null metric. Only a zero matrix residual permits curvature projection. Frame-based computer algebra and its translation between tensor and spinor objects are set out in García-Parrado and Martín-García 2012, § 2 and Eqs. 1–11.

Schwarzschild in a Kinnersley-type frame

Let $M$ be the physical mass, so $GM$ is its geometrized length in the inherited $c=1$ units. For $f=1-2GM/r$ in the exterior Schwarzschild chart, choose

\ell^\mu=\left(f^{-1},1,0,0\right), \qquad n^\mu=\frac12\left(1,-f,0,0\right),

m^\mu=\frac{1}{\sqrt2\,r} \left(0,0,1,i\csc\theta\right).

For example, $\ell\cdot n=-f(f^{-1})(1/2)+f^{-1}(1)(-f/2)=-1$ , and the two angular terms each contribute $1/2$ to $m\cdot\bar m$ . Direct contraction gives the full required normalization on $r\gt2GM$ and $0\lt\theta\lt\pi$ . The poles in $m^\mu$ are frame-coordinate singularities of the spherical angular basis; they are not curvature singularities. With the site convention above, the vacuum Ricci scalars vanish and

\Psi_2=+\frac{GM}{r^3}, \qquad \Psi_0=\Psi_1=\Psi_3=\Psi_4=0.

This pattern encodes Petrov type D in a principal frame. It is stronger than checking only one scalar because an incorrect nonprincipal tetrad can make other $\Psi_i$ nonzero while leaving invariant combinations unchanged. The implementation should compute the projections directly from the Weyl tensor, not insert the expected answer as a replacement rule.

Gauge transformations and weighted quantities

Under a real boost and a spin rotation,

\ell\rightarrow A\ell, \quad n\rightarrow A^{-1}n, \quad m\rightarrow e^{i\chi}m, \quad \bar m\rightarrow e^{-i\chi}\bar m,

the normalization is unchanged. Individual Newman–Penrose quantities transform with boost and spin weights; $\Psi_2$ is invariant under this restricted transformation, whereas $\Psi_0$ and $\Psi_4$ are not. A symbolic environment should record whether $A$ and $\chi$ are constant or coordinate-dependent, because derivatives of a local transformation enter spin coefficients. The implementation architecture and validation examples in García-Parrado and Martín-García 2012, §§ 3–4 illustrate why frame definitions and conversion rules must remain explicit.

The curvature weights can be checked by counting legs in the defining contractions:

\Psi_j\longrightarrow A^{2-j}e^{i(2-j)\chi}\Psi_j, \qquad j=0,1,2,3,4.

There are no derivatives in a curvature projection, so this rule holds pointwise for a local spin or boost as well. It does not apply unchanged to connection coefficients, whose definitions differentiate the tetrad. Because all weighted scalars except $\Psi_2$ vanish in the principal Schwarzschild frame, their unchanged zeros alone cannot test the powers of $A$ or the phases; check the symbolic contraction rule or use a deliberately nonprincipal frame with a nonzero weighted scalar.

Discrete exchanges also matter. Interchanging $\ell$ and $n$ swaps the roles of incoming and outgoing quantities; exchanging $m$ and $\bar m$ complex-conjugates spin-weighted projections. A comparison between packages must translate these operations before interpreting a sign disagreement as an error.

Spin coefficients provide another diagnostic because they depend on both the connection and derivatives of the frame. Compute them from their defining contractions, then test the Newman–Penrose identities appropriate to vacuum rather than trusting a package-specific name. Their values need not remain fixed under a local boost or spin; the predicted inhomogeneous derivative terms must appear. This test catches a common implementation error in which a tetrad transformation is applied algebraically to curvature scalars but not differentiated when the connection coefficients are rebuilt.

Curvature reconstruction as a benchmark

The Schwarzschild benchmark should reconstruct a coordinate-independent quantity from the frame output. With only real $\Psi_2$ nonzero,

R_{\alpha\beta\gamma\delta}R^{\alpha\beta\gamma\delta} =48\Psi_2^2 =\frac{48(GM)^2}{r^6}.

Compare this with a direct coordinate contraction made before the tetrad was introduced. That is an independent representation check: it exercises the frame projection and the coordinate curvature pipeline differently. Also apply several symbolic boosts $A(r)\gt 0$ and spins $\chi(r,\theta)$ ; normalization, $\Psi_2$ , and the reconstructed invariant must remain unchanged, while quantities with nonzero weights must transform as predicted.

The environment record pins the computer-algebra and frame packages, complex-conjugation assumptions, chart domain, tetrad order, Weyl-scalar signs, and exact input metric. A reproducible sequence constructs the metric and tetrad, prints the Gram residual, computes all Ricci and Weyl scalars, applies a declared boost-spin test, reconstructs the Kretschmann scalar, and emits exact residual expressions plus versions. Precision refinement at generic points provides a secondary convergence check for an independent numerical contraction; discrepancies should fall below a stated tolerance when precision rises.

Principal failure modes are accepting a tetrad before normalization, treating $m$ and $\bar m$ as independent real vectors, silently conjugating parameters, and comparing weighted scalars without translating frame gauge. A result confined to the equatorial plane also misses polar frame singularities and angular sign errors. No single package output is sufficient evidence unless the independent coordinate invariant agrees.

Check your understanding

A boost changes a computed $\Psi_4$ by $A^{-2}$ but leaves $\Psi_2$ fixed. Does the changed $\Psi_4$ show that one of the two calculations is wrong?

Solution

Not by itself. $\Psi_4$ has nonzero boost weight, so its numerical value depends on the normalization of the null legs. Translate both tetrads to the same $A$ , verify their Gram matrices, and compare a boost-invariant curvature combination. A disagreement after that translation is meaningful; a raw disagreement before it is gauge.

References

García-Parrado Gómez-Lobo, A., and Martín-García, J. M. (2012), “Spinors: A Mathematica Package for Doing Spinor Calculus in General Relativity,” Computer Physics Communications 183, 2214–2225. DOI.
Newman, E. T., and Penrose, R. (1962), “An Approach to Gravitational Radiation by a Method of Spin Coefficients,” Journal of Mathematical Physics 3, 566–578. DOI.