---
title: Parity-Even Cubic Weyl Operator
url: https://www.emergentmind.com/topics/parity-even-cubic-weyl-operator
type: topic
---

# Parity-Even Cubic Weyl Operator

Searching arXiv for Cano's Teukolsky-based Love number analysis and related context.
arxiv_search.query({"search_query":"au:Cano black holes Love numbers Teukolsky cubic Weyl", "start": 0, "max_results": 10})
arxiv_search.query({"search_query":"all:Teukolsky Love numbers cubic Weyl Cano", "start": 0, "max_results": 10})
arxiv_search.query({"search_query":"all:\"Weyl\" AND all:\"Love number\" AND au:Cano", "start": 0, "max_results": 10})
The parity-even cubic Weyl operator is the cubic contraction of the Weyl tensor,
\[
C^3 \equiv C_{\mu\nu}{}^{\rho\sigma}C_{\rho\sigma}{}^{\alpha\beta}C_{\alpha\beta}{}^{\mu\nu},
\]
also written as \(W^3\) in cosmological applications. In the four-dimensional gravitational effective field theory considered for Schwarzschild black holes, it is the unique parity-even cubic Weyl invariant, while in inflationary perturbation theory it appears as a parity-conserving higher-derivative graviton interaction with distinctive bispectral selection rules [2606.16909, 1108.0175]. Across these settings, the operator is introduced as a cutoff-suppressed correction to Einstein gravity and serves as a probe of finite-size effects, higher-curvature response, and tensor non-Gaussianity.

## 1. Algebraic definition and symmetry properties

In the black-hole EFT formulation, the bulk action is
\[
S_{\rm bulk}=\frac{1}{16\pi G}\int d^4x\,\sqrt{-g}\,\Bigl[R+\lambda_e\Lambda^{-4}\,
C_{\mu\nu}{}^{\rho\sigma}C_{\rho\sigma}{}^{\alpha\beta}C_{\alpha\beta}{}^{\mu\nu}\Bigr],
\]
with cutoff \(\Lambda\) and Wilson coefficient \(\lambda_e\) [2606.16909]. The same cubic structure appears in the cosmological literature as
\[
L_{W^3}(\tau)=\frac{f(\tau)}{\Lambda^2}\,
W^{\mu\nu}{}_{\rho\sigma}W^{\rho\sigma}{}_{\alpha\beta}W^{\alpha\beta}{}_{\mu\nu},
\qquad
f(\tau)=\Bigl(\frac{\tau}{\tau_*}\Bigr)^A,
\]
where the Weyl tensor is
\[
W_{\mu\nu\rho\sigma}\equiv R_{\mu\nu\rho\sigma}
-\frac{1}{2}\bigl(g_{\mu[\rho}R_{\sigma]\nu}-g_{\nu[\rho}R_{\sigma]\mu}\bigr)
+\frac{R}{6}g_{\mu[\rho}g_{\sigma]\nu}
\]
[1108.0175].

By construction, \(W^3\) is parity-even. In the cosmological treatment it is explicitly described as parity-conserving and is contrasted with the parity-odd operator \(\tilde W W^2\), whose observational imprint occupies complementary CMB multipole sectors [1108.0175]. This separation is structural rather than conventional: the parity assignment follows from the tensor contraction itself.

## 2. Effective-field-theory role and perturbative control

For Schwarzschild black holes, the relevant small parameter is
\[
\epsilon_e\equiv \lambda_e(\Lambda r_s)^{-4}\ll 1,
\qquad r_s=2GM,
\]
so the cubic Weyl correction is treated perturbatively around the general-relativistic background [2606.16909]. The background metric ansatz to \(O(\epsilon_e)\) is
\[
ds^2=-A(r)dt^2+\frac{dr^2}{B(r)}+r^2d\Omega^2,
\]
with
\[
A(r)=f(r)\,[1+\epsilon_e a(r)],\qquad
B(r)=f(r)\,[1+\epsilon_e b(r)],\qquad
f(r)=1-\frac{r_s}{r}.
\]

In the inflationary setting, the same operator is introduced on an exact de Sitter background with conformal time \(\tau\), and its strength is suppressed by \(\Lambda^{-2}\) together with a time-dependent coupling \(f(\tau)=(\tau/\tau_*)^A\) [1108.0175]. The normalization differs from the black-hole EFT normalization. This suggests that the operator’s physical interpretation is stable across applications, while the precise power of the cutoff and the coupling convention are context-dependent.

The EFT logic is similar in both cases: \(W^3\) encodes higher-derivative corrections beyond the Einstein-Hilbert term, and the perturbative expansion isolates its leading observable consequences. In the black-hole problem those consequences appear in static tidal response; in the cosmological problem they appear in a graviton three-point function.

## 3. Static even-parity quadrupole sector on Schwarzschild

The black-hole analysis focuses on static even-parity \(\ell=2\) perturbations in Regge-Wheeler gauge,
\[
\delta g_{tt}=-A(r)\eta H_0(r)P_2(\mu),\qquad
\delta g_{rr}=B(r)^{-1}\eta H_2(r)P_2(\mu),
\]
\[
\delta g_{\theta\theta}=r^2\eta K(r)P_2(\mu),\qquad
\delta g_{\phi\phi}=r^2\sin^2\theta\,\eta K(r)P_2(\mu),
\]
with first-order \(\epsilon_e\) expansions
\[
H_0=-H(r)+X_0(r)+O(\epsilon_e^2),\qquad
H_2=H(r)+X_2(r)+O(\epsilon_e^2),
\]
\[
K=K_{\rm GR}[H](r)+X_K(r)+O(\epsilon_e^2).
\]
Here \(H(r)\) and \(K_{\rm GR}[H]\) satisfy the \(\ell=2\) Zerilli/Regge-Wheeler static equations in pure GR [2606.16909].

Substituting these fields into the bulk action and expanding to quadratic order in the perturbation amplitude \(\eta\) and to first order in \(\epsilon_e\) yields a one-dimensional reduced radial Lagrangian,
\[
L_2^{\rm radial}
=
L_{2,EH}^{(0)}[H_0,H_2,K]
+\epsilon_e\,L_{2,EH}^{(1)}[H_0,H_2,K;a,b]
+r_s^4\,L_{2,C^3}^{(0)}[H_0,H_2,K].
\]
The three pieces have distinct origins: the GR quadratic action on Schwarzschild, the correction induced by the \(O(\epsilon_e)\) change in the background, and the Weyl-cubic term evaluated on the unperturbed background but to quadratic order in the GR tidal perturbation [2606.16909].

A key technical feature is that higher-order radial derivatives, up to fourth order, appear only at \(O(\epsilon_e)\). They are then eliminated by order reduction using the zeroth-order GR tidal equations. With the generalized Euler-Lagrange operator
\[
\frac{\delta L}{\delta q}
=
\sum_{k=0}^{N}(-1)^k\frac{d^k}{dr^k}
\left(\frac{\partial L}{\partial q^{(k)}(r)}\right),
\qquad q\in\{H_0,H_2,K\},\quad N\le 4,
\]
the \(O(\epsilon_e)\) field equations become a linear inhomogeneous system for \(X_0,X_2,X_K\). One equation is purely algebraic in \(X_2\), fixing
\[
X_2(r)=\mathcal{C}_2[X_0,X_K;r].
\]
After substitution, the remaining system closes on the two-vector
\[
X(r)\equiv [X_0(r),X_K(r)]^T,
\]
which satisfies
\[
\frac{1}{r}\frac{dX}{dr}=A_2(r)\,X+b_2(r),
\]
with \(A_2(r)\) and \(b_2(r)\) rational functions of \(r/r_s\) [2606.16909]. The explicit radial integrand is described as bulky and was generated with Mathematica; the full expressions for \(A_2(r)\), \(b_2(r)\), and the cubic invariant on Schwarzschild are provided in the supplemental Mathematica notebook.

## 4. Boundary conditions and the fixed-quadrupole response benchmark

Near the horizon, with \(r=r_s+\rho\), the residue data of the homogeneous system are
\[
A_{\rm hor}=\operatorname{Res}_{r=r_s}A_2=
\begin{bmatrix}
-1 & 0\\
0 & 0
\end{bmatrix},
\qquad
b_{\rm hor}=(-144,0)^T,
\]
so the two local exponents are \(\{-1,0\}\) [2606.16909]. Regularity on the future horizon removes the \((-1)\) mode and leaves one free homogeneous datum,
\[
q_0\equiv X_K(r_s).
\]

At asymptotic infinity, a Laurent expansion,
\[
X_0=\sum_n a_n r^n,\qquad X_K=\sum_n k_n r^n,
\]
reveals a free growing-branch constant proportional to \(\epsilon_e\), denoted \(\alpha\), with \(\alpha=r_s a_1\). The large-\(r\) behavior is
\[
X_0\sim -\frac{24+\alpha}{r_s^2}r^2+\cdots+192\,r_s^3 r^{-3}+\cdots,
\]
\[
X_K\sim \frac{24+\alpha}{r_s^2}r^2+\cdots-480\,r_s^3 r^{-3}+\cdots,
\]
\[
X_2\sim \frac{24+\alpha}{r_s^2}r^2+\cdots-480\,r_s^3 r^{-3}+\cdots.
\]
The \(r^2\) terms renormalize the applied tidal field, while the \(r^{-3}\) terms define the induced response. Pure GR has \(\alpha=0\) but no decaying branch [2606.16909].

Matching the horizon and asymptotic expansions gives the degeneracy relation
\[
\alpha=144+2q_0.
\]
The remaining horizon datum therefore only shifts the growing solution. Imposing the no-tidal-renormalization convention,
\[
\alpha+24=0\quad\Rightarrow\quad \alpha=-24,\qquad q_0=-84,
\]
removes this ambiguity and fixes the decaying branch:
\[
X_0^{(-3)}=192\,r_s^3,\qquad
X_K^{(-3)}=X_2^{(-3)}=-480\,r_s^3.
\]
Restoring the metric-level definitions gives
\[
H_2^{(1)}[-3]=K^{(1)}[-3]=-480\,r_s^3.
\]

Calibrating the spatial sector at fixed \(\ell=2\) against the associated-Legendre branches \(P_2^2\) and \(Q_2^2\), and using
\[
Q_2^2\sim \frac{r_s^3}{5}r^{-3},
\]
the fixed-quadrupole response amplitude becomes
\[
\Delta(B/A)_{\ell=2}^{\rm fix}=-2400\,\epsilon_e.
\]
Equivalently, the scalar fixed-\(\ell\) quotient is
\[
\Delta k_{2,\rm sc}^{\rm fix}=-20\,\epsilon_e
\]
[2606.16909]. In the terminology of that work, these are the primary metric-sector results at \(O(\epsilon_e)\).

## 5. Relation to Love numbers and interpretive caveats

A central interpretive issue is that the scalar quantity \(\Delta k_{2,\rm sc}^{\rm fix}\) is not, by itself, the analytically continued, gauge-invariant electric Love number [2606.16909]. At fixed integer \(\ell\), the separation between physical response, tidal-field redefinition, and gauge artifacts is not unique. In the metric derivation, this ambiguity is encoded in the freedom \(q_0\) or, equivalently, \(\alpha\), which reintroduces a pure growing solution unless a convention is imposed.

The canonical extraction of a gauge-invariant Love number \(k_\ell^+\) requires three further steps: working in a gauge-invariant master-variable framework, such as Teukolsky; performing an analytic continuation in \(\ell\), \(\ell\to 2+\delta\), and isolating the finite, non-running part as \(\delta\to 0\); and matching normalizations carefully [2606.16909]. A Teukolsky-based analysis cited there finds, for the same parity-even \(C^3\) operator,
\[
k_2^+=28\,\lambda_{ev}/M^4,\qquad
k_2^-=-20\,\lambda_{ev}/M^4.
\]
The same source notes that a rough identification \(\lambda_{ev}M^4\simeq 16\epsilon_e\), obtained by setting \(r_s=2M\), would suggest
\[
\Delta k_{2,\rm sc}^{\rm fix}\simeq -1.25\,\lambda_{ev}/M^4,
\]
which differs from the canonical \(k_2^+\). The mismatch is expected because the fixed-\(\ell\) metric quotient and the analytically continued gauge-invariant Love number are different observables.

Within this framework, the parity-even cubic Weyl operator breaks the exact GR cancellation that makes the static \(\ell=2\) Love numbers of a Schwarzschild black hole vanish, and it induces a genuine decaying response at \(O(\epsilon_e)\) [2606.16909]. The correction enters at relative order \((r_s/\Lambda)^4\) and remains small so long as \(r_s\ll \Lambda^{-1}\).

## 6. Primordial graviton bispectra and CMB signatures

In the cosmological application, the parity-even cubic Weyl operator generates a graviton interaction Hamiltonian obtained from
\[
H_{\rm int}(\tau)=-\int d^3x\,L_{W^3}(\tau),
\]
after inserting the all-orders expansion of the metric and retaining cubic terms in the transverse-traceless graviton field \(\gamma_{ij}\) [1108.0175]. The graviton bispectrum is then computed at tree level with the in-in formalism,
\[
\langle \gamma_{\lambda_1}(k_1)\gamma_{\lambda_2}(k_2)\gamma_{\lambda_3}(k_3)\rangle
=
i\int_{-\infty}^{0}d\tau\,
\langle 0|[H_{W^3}(\tau),\gamma_{\lambda_1}(k_1)\gamma_{\lambda_2}(k_2)\gamma_{\lambda_3}(k_3)]|0\rangle.
\]

The resulting three-point function is parity-even and, in the notation of the inflationary analysis, scales as
\[
B_{\lambda_1\lambda_2\lambda_3}(k_i)
\propto
(H/M_{\rm pl})^6(H/\Lambda)^2
\cos\!\Bigl(\frac{\pi A}{2}\Bigr)\Gamma(6+A)\,
S_{W^3}^{(A)}(k_i),
\]
where \(S_{W^3}^{(A)}\) is the reduced shape function [1108.0175]. In exact de Sitter, no slow-roll suppression appears. The operator therefore produces a purely parity-even primordial bispectrum even in that limit.

Projected into the CMB, the reduced bispectrum obeys a parity-even selection rule:
\[
\ell_1+\ell_2+\ell_3=\text{even}.
\]
In particular, the temperature bispectrum \(b^{III}_{\ell\ell\ell}\) is nonzero only when \(3\ell\) is even, i.e. \(\ell\) even [1108.0175]. By contrast, the parity-odd operator \(\tilde W W^2\) contributes only when \(\ell_1+\ell_2+\ell_3\) is odd. This cleanly separates parity-even and parity-odd cubic Weyl interactions at the level of CMB multipole configurations.

At large scales, the temperature bispectrum from \(W^3\) is estimated numerically as
\[
|b^{III}_{\ell\ell\ell}|
\sim
\ell^{-4}\times 3.2\times 10^{-2}
\left(\frac{\rm GeV}{\Lambda}\right)^2
\left(\frac{r}{0.1}\right)^4,
\]
and comparison with an equilateral-type nonlinearity bound \(f_{\rm NL}^{\rm eq}<300\) yields
\[
\Lambda\gtrsim 3\times 10^6\ {\rm GeV},
\]
assuming \(r=0.1\) [1108.0175]. In this setting, the parity-even cubic Weyl operator functions as a higher-derivative source of tensor non-Gaussianity with a specific angular-parity signature.

## 7. Scope, significance, and common points of confusion

The operator’s significance depends on which observable is being extracted. In the black-hole problem, the central result is a fixed-quadrupole metric-sector benchmark,
\[
\Delta(B/A)_{\ell=2}^{\rm fix}=-2400\,\epsilon_e,
\]
obtained after removing tidal renormalization by convention [2606.16909]. In the cosmological problem, the same parity-even cubic Weyl structure is important because it produces a parity-even primordial graviton bispectrum and contributes only to even \(\ell_1+\ell_2+\ell_3\) CMB bispectra [1108.0175].

A common misunderstanding is to treat all “Love-number-like” quotients extracted at fixed \(\ell\) as gauge-invariant tidal Love numbers. The black-hole analysis explicitly rejects that identification: \(\Delta k_{2,\rm sc}^{\rm fix}\) is a scalar fixed-\(\ell\) conversion of a metric branch ratio and not, by itself, the analytically continued gauge-invariant electric Love number [2606.16909]. Another possible confusion is to equate “parity-even” with the absence of observationally distinctive parity structure. The CMB analysis shows the opposite: parity conservation forces the signal into even total multipole configurations, while the parity-odd cubic Weyl interaction occupies the odd sector [1108.0175].

Taken together, these results establish the parity-even cubic Weyl operator as a concrete higher-curvature deformation with two sharply defined uses. In black-hole EFT it provides a reproducible metric-space benchmark for static tidal response. In primordial cosmology it provides a parity-conserving source of graviton non-Gaussianity with specific CMB selection rules.

Source: https://www.emergentmind.com/topics/parity-even-cubic-weyl-operator