Weyl Scalars: Key Aspects in Gravitational Physics

• Weyl scalars are five complex quantities that decompose spacetime curvature into radiative, Coulombic, and longitudinal components, providing key insights into gravitational fields.
• Their radial falloff hierarchy, established by the peeling theorem, supports precise gravitational wave extraction in numerical relativity, with deviations around 4%.
• Canonical methods leverage Weyl scalars as intrinsic observables to label spacetime points, influencing analyses in binary black holes and compact object studies.

The Weyl scalars are five complex scalar quantities, denoted $\Psi_0$ through $\Psi_4$, derived by projecting the Weyl tensor onto a null tetrad in four-dimensional spacetime. They constitute a natural decomposition of the algebraically independent, conformally invariant part of the curvature, isolating radiative, Coulombic, and longitudinal features of the gravitational field. The framework generalizes across exact solutions, perturbation theory, and numerical relativity, serving as the backbone of modern gravitational radiation theory, Petrov algebraic classification, and canonical approaches to general covariance.

1. Formal Definition and Algebraic Structure

Given a vacuum spacetime with Weyl tensor $C_{abcd}$ and a Newman–Penrose (NP) null tetrad $\{\ell^a, n^a, m^a, \bar m^a\}$, satisfying

$$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$

the Weyl scalars are defined by the explicit contractions: \begin{align*} \Psi_0 &= C_{abcd}\,\ell^a\,m^b\,\ell^c\,m^d, \ \Psi_1 &= C_{abcd}\,\ell^a\,n^b\,\ell^c\,m^d, \ \Psi_2 &= C_{abcd}\,\ell^a\,m^b\,\bar m^c\,n^d, \ \Psi_3 &= C_{abcd}\,\ell^a\,n^b\,\bar m^c\,n^d, \ \Psi_4 &= C_{abcd}\,n^a\,\bar m^b\,n^c\,\bar m^d. \end{align*} These quantities are the spin-frame components of the totally symmetric Weyl spinor $\Psi_{ABCD}$, foundational in the NP and two-spinor formalism (Hinder et al., 2011, Han, 2022).

In static, spherically symmetric or axisymmetric contexts, the full contraction $\mathcal W = C^{abcd}C_{abcd}$—sometimes referred to as the Weyl scalar in the scalar-invariant sense—also plays a central role, particularly in compact object phenomenology (Danarianto et al., 5 Oct 2025).

2. The Peeling Theorem and Radial Falloff

In asymptotically flat spacetimes, the peeling theorem establishes a strict radial decay hierarchy for the Weyl scalars along outgoing null directions. Assuming the existence of a conformal compactification $\tilde g_{ab} = \Omega^2 g_{ab}$ with $\Omega = 0$ at future null infinity $\mathscr I^+$ and suitable regularity of the rescaled curvature, the theorem states: $\Psi_4$0 along outgoing radial null geodesics parametrized by $\Psi_4$1. Thus,

$\Psi_4$2

This hierarchy underpins both exact and perturbative descriptions of gravitational radiation and sources the asymptotic structure of the field (Hinder et al., 2011). In numerical relativity (NR), direct simulation of binary black hole (BBH) spacetimes has confirmed the predicted falloffs to within $\Psi_4$3, validating the theoretical picture even with approximate null rays and coordinate-based tetrads.

3. Physical Interpretation in the Petrov and NP Formalisms

The $\Psi_4$4 each encode distinct aspects of the spacetime curvature:

• $\Psi_4$5: Represents the Coulombic (Newtonian) part—mass monopole and spin dipole. Survives in stationary, axisymmetric fields and is nonzero for Petrov type D solutions.
• $\Psi_4$6 and $\Psi_4$7: Encode transverse gravitational radiation. $\Psi_4$8 represents incoming radiation at $\Psi_4$9; $C_{abcd}$0 is outgoing. At large $C_{abcd}$1, energy flux and Bondi news at null infinity are constructed from $C_{abcd}$2.
• $C_{abcd}$3 and $C_{abcd}$4: Carry longitudinal field components—less dominant at large distances. These vanish in algebraically special spacetimes except for specific degenerate types.

In the spinor formalism, the multiplicity and structure of the principal spinors in the decomposition $C_{abcd}$5 directly reflect the Petrov type: type N, D, or III (Han, 2022).

4. Canonical and Covariant Constructions

The canonical (Hamiltonian) formulation of general relativity supports the explicit construction of (dynamically gauge-invariant) Weyl scalars as functions of phase-space variables. In the 3+1 Arnowitt–Deser–Misner (ADM) or York canonical basis, scalar invariants

$C_{abcd}$6

are functions only of the spatial metric $C_{abcd}$7 and their conjugate momenta $C_{abcd}$8. They serve as intrinsic, diffeomorphism-invariant coordinate fields $C_{abcd}$9, generalizing the proposal of Bergmann and Komar for labeling spacetime points by curvature (Watson et al., 2024). This enables a fully covariant gauge-fixing and is instrumental in canonical quantization programs.

In the York framework, null tetrads built from the canonical variables lead to the Hamiltonian analogues of the NP Weyl scalars, and the radar-tensor formalism supplies manifestly four-scalar components under general coordinate transformations (Lusanna et al., 2014).

5. Computation and Practical Extraction in Numerical Relativity

Standard NR codes operationalize extraction of the Weyl scalars using coordinate-based null tetrads constructed via Gram–Schmidt orthonormalization from the unit normal $\{\ell^a, n^a, m^a, \bar m^a\}$0, radial vector $\{\ell^a, n^a, m^a, \bar m^a\}$1, and two orthonormal angular vectors. The corresponding NP null tetrad is

$\{\ell^a, n^a, m^a, \bar m^a\}$2

This construction is robust under most numerical schemes, though special care is required near axis singularities and in the presence of numerical noise, particularly for the rapidly decaying $\{\ell^a, n^a, m^a, \bar m^a\}$3 and $\{\ell^a, n^a, m^a, \bar m^a\}$4 (Hinder et al., 2011, Iozzo et al., 2020).

Current best practices involve electric-magnetic (E/B) decompositions of $\{\ell^a, n^a, m^a, \bar m^a\}$5 on extraction spheres, projection onto analytic or numerically precomputed $\{\ell^a, n^a, m^a, \bar m^a\}$6 fields, mode decomposition in spin-weighted harmonics, and multi-radius extrapolation to asymptotic null infinity. This pipeline underlies wave extraction in collaborations such as SXS and in codes like SpEC (Iozzo et al., 2020). Convergence, numerical stability, and validation against analytic solutions (e.g., Kerr) are now routine.

6. Weyl Scalars as Observables and Intrinsic Coordinates

Weyl scalars, as pointwise curvature invariants, are simultaneously local probes of the gravitational field and candidates for intrinsic spacetime labeling variables. In canonical language, any complete set of four functionally independent scalars built from the Weyl tensor (and thus from $\{\ell^a, n^a, m^a, \bar m^a\}$7) can serve as coordinate functions $\{\ell^a, n^a, m^a, \bar m^a\}$8—intrinsic coordinates—in covariance-preserving gauge fixings (Watson et al., 2024, Lusanna et al., 2014). However, a subtlety arises: the algebraic (Bergmann–Komar) invariants constructed from the Weyl tensor are not always Dirac observables, as their explicit form typically depends on gauge variables (lapse, shift, conformal gauge), precluding direct identification as predictive observables except in constraint-satisfying (on-shell) settings (Lusanna et al., 2014).

7. Applications in Compact Objects and Binary Mergers

The scalar invariant $\{\ell^a, n^a, m^a, \bar m^a\}$9—the quadratic contraction of the Weyl tensor—has emerged as a robust, nearly equation-of-state–independent diagnostic of interior and surface tidal distortion in neutron stars and quark stars: $$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$0 with vanishing central value, surface limit $$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$1, and tight log–log fits to normalized moment of inertia and tidal deformability across microphysical models (Danarianto et al., 5 Oct 2025). These relations enable direct mapping of gravitational-wave constraints into spacetime-curvature bounds inside neutron stars, demonstrating the astrophysical impact of Weyl scalar diagnostics in strong-field gravity.

In binary black hole spacetimes, numerical confirmation of the radial falloff exponents of all five Weyl scalars along near–null directions provides a definitive test of the peeling theorem and supports the standard interpretation of $$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$2 as the gravitational wave observable in the wave zone (Hinder et al., 2011). Zone analysis further quantifies where higher-$$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$3 Weyl scalars contribute non-negligibly to curvature, refining the division into near, transition, and radiation zones crucial for waveform modeling.

---

References:

• (Hinder et al., 2011): Falloff of the Weyl scalars in binary black hole spacetimes
• (Han, 2022): Weyl double copy and massless free fields in curved spacetimes
• (Iozzo et al., 2020): Extending Gravitational Wave Extraction Using Weyl Characteristic Fields
• (Watson et al., 2024): An Intrinsic Coordinate Reference Frame Procedure I: Tensorial Canonical Weyl Scalars
• (Danarianto et al., 5 Oct 2025): $$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$4Love$$\ell_a \ell^a = n_a n^a = m_a m^a = 0, \qquad \ell_a n^a = -1, \quad m_a \bar m^a = +1,$$5Curvature: Exploring compact stars' quasi-universal relation with curvature scalars
• (Lusanna et al., 2014): Hamiltonian Expression of Curvature Tensors in the York Canonical Basis: II) The Weyl Tensor, Weyl Scalars, the Weyl Eigenvalues and the Problem of the Observables of the Gravitational Field