---
title: Fully Gravitational Wave Function
url: https://www.emergentmind.com/topics/fully-gravitational-wave-function
type: topic
---

# Fully Gravitational Wave Function

Fully gravitational wave function denotes a family of constructions in which gravitational radiation or gravitational states are formulated using fully dynamical geometric data rather than a fixed background, a special gauge choice, or externally frozen boundary data. The literature uses the expression, and closely related ones, in several technically distinct settings: curvature-based detector response in classical general relativity, second-order nonlinear wave equations built from the full metric, Hartle–Hawking states with dynamical boundaries in AdS, representation-theoretic bulk states in JT supergravity, and operational waveform models for fully-precessing eccentric binaries [1310.2871], [2010.08837], [2605.13970], [2405.09289], [2402.06804]. This suggests that the unifying idea is not a single canonical object but a programmatic requirement: gravity is to be treated as fully dynamical at the level where the wave description is defined.

## 1. Terminological scope and common theme

In the detector-response literature, the fully gravitational viewpoint replaces gauge-dependent metric perturbations by the Riemann curvature tensor and measured worldline data, so that the response is “separately gauge invariant and has a clear physical interpretation” [1310.2871]. In nonlinear perturbation theory, the same phraseology points to a description in which \(h_{\mu\nu}\) is not treated as a field on a passive Minkowski background, but as part of the geometry itself, with indices raised and lowered by the full metric and the harmonic gauge imposed consistently to second order [2010.08837].

In quantum gravity, the phrase becomes more literal. For the Hartle–Hawking wave function in AdS, a fully gravitational construction is one in which the extra spacetime boundary induced by an open spatial slice is itself dynamical and integrated over; the contrasting partially frozen construction fixes that boundary as in AdS/CFT [2605.13970]. In JT supergravity, “gravitational wavefunctions” are specific mixed parabolic matrix elements of a positive semigroup version of the relevant supergroup, chosen so that asymptotic AdS boundary conditions are enforced directly at the quantum level [2405.09289].

A further, operationally distinct usage appears in waveform modeling. The EFPE and EFPE\(_{\rm ME}\) families provide efficient frequency-domain inspiral waveforms for compact binaries that are both spin-precessing and eccentric, and the synthesis of that work describes the result as a “fully gravitational-wave function” in the sense of a waveform model that retains both effects simultaneously [2402.06804]. The coexistence of these usages indicates that the term has no single universal definition across gravitational physics.

## 2. Curvature-based classical description

A central classical formulation holds that gravitational-wave detectors do not fundamentally measure metric perturbations \(h_{ab}\); they measure the effect of spacetime curvature on light propagation [1310.2871]. The detector class considered is broad: any detector that measures the phase of a remote clock by means of an electromagnetic signal propagating along null geodesics. In that framework the observable is written in terms of intrinsic clock evolution, the accelerations of clock and receiver, a kinematic light-time term, and a curvature term containing the null-path integral
\[
\int_0^{\lambda_R} R(\mathbf k,\mathbf U,\mathbf k,\mathbf U)\, d\lambda .
\]

In the short-wavelength or high-frequency limit, the Riemann tensor separates into a background curvature contribution and a wave curvature contribution. The perturbation \(R^{(1)}_{abcd}\) is then gauge invariant up to terms suppressed by \(\delta^2\), with \(\delta=(\omega L)^{-1}\ll 1\), and the gravitational-wave response reduces to the integral of the wave curvature along the unperturbed null geodesic:
\[
\left(\frac{d^2\phi}{d\tau^2}\right)^{(1)}
=
\frac{d\phi/dt}{\mathbf U\cdot\left(\frac{d\phi}{d\tau}\right)^2}
\int_0^{\lambda_R} R^{(1)}(\mathbf k,\mathbf U,\mathbf k,\mathbf U)\, d\lambda
+\mathcal O(\delta).
\]
This applies, according to the paper, to ground-based interferometers, LISA-like detectors, spacecraft Doppler tracking, and pulsar timing arrays [1310.2871].

A related fully gravitational characterization of radiation in full nonlinear GR is formulated at null infinity rather than in the bulk. In asymptotically Minkowskian spacetimes, the radiative degrees of freedom are not encoded in the universal structure \((q_{ab},n^a)\) of \(\scri\), but in the derivative operator \(\mathcal D\) induced on \(\scri\). The gauge-invariant radiative data are carried by the trace-free shear \(\sigma_{ab}\), and the Bondi news satisfies
\[
N_{ab}=2\,\mathcal L_n \sigma^\circ_{ab}=\dot{\sigma}^\circ_{ab},
\]
while the radiative Weyl scalar obeys
\[
\Psi^\circ_4=-\ddot{\bar{\sigma}}^\circ .
\]
In the linearized limit this becomes
\[
\Psi^\circ_4=-\frac12\left(\ddot h^\circ_+ - i\,\ddot h^\circ_\times\right),
\]
so the familiar strain polarizations arise as the asymptotic shadow of the nonlinear shear [2201.11634].

A recurrent misconception addressed by these formulations is that a detector’s operation is fundamentally explained by a gauge-fixed metric perturbation, for example in TT gauge. The curvature-based account instead treats TT-gauge formulas as special cases and relocates the invariant content to curvature transport and null-infinity radiative data [1310.2871], [2201.11634].

## 3. Nonlinear wave descriptions based on the full metric

A second use of the fully gravitational idea appears in second-order perturbation theory. The key change is to write
\[
g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu},
\]
but to raise and lower indices with the full inverse metric expanded consistently to second order,
\[
g^{\mu\nu}=\eta^{\mu\nu}-h^{\mu\nu}+h^\mu{}_\rho h^{\rho\nu},
\qquad O(h^2),
\]
and only then discard terms of order \(h^3\) and above [2010.08837].

The motivation is that once nonlinear curvature terms are retained, the perturbation contributes to the geometry, so traces, contractions, and gauge restrictions must be defined at the same perturbative order. The Ricci tensor is expanded to \(O(h^2)\) using Christoffel symbols built from the full metric, vacuum is imposed through \(R_{\mu\nu}=0\), and the full harmonic gauge condition
\[
\partial_\mu\!\left(\sqrt{-g}\,g^{\mu\nu}\right)=0
\]
is also expanded consistently to second order. The resulting equation is a nonlinear vacuum wave equation of schematic form
\[
\Box h_{\mu\nu}=\text{quadratic terms in }h\text{ and its derivatives},
\]
with self-couplings and derivative couplings that are absent in the usual linear plane-wave treatment [2010.08837].

Within that framework, the perturbation ansatz
\[
h_{\mu\nu}(x)=a_{\mu\nu}(x)e^{ikx}+a_{\mu\nu}^*(x)e^{-ikx}
\]
leads to the claim that non-trivial solutions are necessarily non-plane-wave modes with non-zero trace when the amplitude varies. The paper further argues that allowed nonlinear solutions generically contain both longitudinal and transverse components, and that the second-order harmonic gauge does not preserve the standard linear TT reduction. These claims are specific to that second-order full-metric construction and depart from the familiar transverse-traceless plane-wave picture [2010.08837].

This creates an important distinction within the literature. Null-infinity treatments identify radiative content through shear, news, and peeling in full GR [2201.11634], whereas the second-order full-metric perturbative program emphasizes nonlinear trace structure and longitudinal-transverse mixing [2010.08837]. The two are not formulated at the same level of description, and conflating them obscures what each is designed to capture.

## 4. Single-equation formulations on curved backgrounds

A different route toward a fully gravitational description seeks a single master equation for perturbations on curved backgrounds. In type D spacetimes, one starts from the Newman–Penrose formalism, where the radiative gravitational variables are the extreme Weyl scalars \(\Psi_0\) and \(\Psi_4\), and Teukolsky’s equations provide two decoupled equations for them [2012.03203].

The proposed unification introduces a transformation function \(E_{p-s}\) and a perturbation variable \(\psi_p=E_{p-s}\Phi_p\), allowing the spin-\(\pm s\) equations to be rewritten in a common form. The geometrical content is then packaged into the “spin-coefficient connection”
\[
\Gamma_p = 2p\,d^\mu S_\mu,
\qquad
S_\mu = (-\gamma,\,-\epsilon-\rho,\,\alpha,\,\beta+\tau).
\]
Using NP identities, the unified perturbation equation becomes
\[
\Big[(\nabla^\mu+p\,\Gamma^\mu)(\nabla_\mu+p\,\Gamma_\mu)-4p^2\Psi_2+\mathcal E R\Big]\Phi_p
=
2\kappa^2 T_p,
\]
with source-free form
\[
\Big[(\nabla^\mu+p\,\Gamma^\mu)(\nabla_\mu+p\,\Gamma_\mu)-4p^2\Psi_2+2R\Big]\Phi_p=0.
\]

In this formulation, the kinetic operator describes propagation with a built-in spin connection, the term \(-4p^2\Psi_2\) encodes the background tidal field in type D spacetimes, and the Ricci scalar term captures scalar-curvature effects. The same structural equation is claimed to govern massless fields of nonzero spin \(s<2\), including Weyl, Maxwell, and Rarita–Schwinger fields, which the paper interprets as analogue models of gravitational waves [2012.03203].

This program is not a quantum-mechanical wave function in the Hartle–Hawking or JT sense. Rather, it is a unified classical field equation whose “fully gravitational” character lies in making the background geometry explicit through \(\Psi_2\), \(R\), and \(\Gamma_\mu\), and in avoiding dependence on a specific coordinate system or a specific metric form within the type D class [2012.03203].

## 5. Fully gravitational Hartle–Hawking constructions in AdS

For the Hartle–Hawking wave function in AdS, open spatial slices require an additional spacetime boundary \(\mathcal B_-\) meeting the slice \(\Sigma\) at a corner \(C=\Sigma\cap\mathcal B_-\). This produces two inequivalent ensembles. The fully gravitational wave function integrates over the geometry of \(\mathcal B_-\),
\[
\Psi_{\rm HH}[\gamma_{ab}]
=
\int_{\mathcal B_-}\! h_{ij}
\int_{\partial_-=\Sigma\cup\mathcal B_-}\mathcal D g_{\mu\nu}\,
e^{-I_E[g_{\mu\nu},h_{ij},\gamma_{ab}]},
\]
or equivalently fixes the metric on \(\Sigma\) while summing over bulk geometries with \(\partial_-\supset\Sigma\). The partially frozen version instead holds \(\mathcal B_-\) fixed, as in an AdS/CFT Dirichlet ensemble [2605.13970].

In AdS\(_3\) Einstein gravity, the fully gravitational action includes the bulk Einstein–Hilbert term, Gibbons–Hawking–York terms on \(\Sigma\) and \(\mathcal B_-\), a tension term on \(\mathcal B_-\), and a Hayward corner term at \(C\). Varying the metric on \(\mathcal B_-\) yields a Neumann-type condition that reduces in \(D=3\) to
\[
K_{ij}\big|_{\mathcal B_-}=T\,h_{ij},
\]
so the boundary geometry is not fixed externally but determined by the saddle-point equations. Explicit minisuperspace analyses are carried out in AdS\(_3\) Einstein gravity and AdS\(_2\) Jackiw–Teitelboim gravity, and in both cases real saddles exist only below a critical boundary size, beyond which the relevant saddle becomes complex [2605.13970].

The one-loop distinction is especially sharp. For the hyperbolic-ball partition function with a dynamical boundary, the trace operator on the boundary sphere has \(D+1\) negative modes, giving
\[
Z_{\rm hyperbolic\,ball}^{\text{1-loop}}
=
(\mp i)^{D+1}\times(\text{positive number}),
\]
whereas with Dirichlet boundary conditions the one-loop answer remains real and positive:
\[
Z_{\rm hyperbolic\,ball,\,DBC}^{\text{1-loop}}
=
(\text{positive number}).
\]
A converse de Sitter calculation with a fixed equator also yields a real positive answer after nontrivial phase cancellation. The paper’s conclusion is that the phase problem is controlled by whether the path integral is fully dynamical or partially frozen, not by AdS versus dS alone [2605.13970].

## 6. Gravitational wavefunctions in JT supergravity

In supersymmetric JT gravity, gravitational wavefunctions are constructed representation-theoretically rather than by solving a minisuperspace differential equation in isolation. For \(\mathcal N=2\) JT supergravity, the relevant structure is the positive semigroup
\[
\mathrm{OSp}^+(2|2,\mathbb R)\subset \mathrm{OSp}(2|2,\mathbb R),
\]
retaining only the hyperbolic conjugacy class so as to match smooth asymptotically AdS Euclidean bulk geometries [2405.09289].

The basic objects are mixed parabolic matrix elements
\[
\langle \nu,\psi_1,\psi_2 \mid g \mid \lambda,\psi_3,\psi_4\rangle,
\]
where the left and right boundary states are Whittaker-type eigenvectors of the parabolic generators. Because the fermionic generators anticommute, the construction uses Majorana/Clifford variables \(\psi_i\) satisfying
\[
\{\psi_i,\psi_j\}=2\delta_{ij},
\]
so that the fermionic eigenvalue problem is encoded in a Clifford algebra. After fixing the parabolic data, the bulk dependence reduces to the Cartan variables, and the nontrivial wavefunction is the Whittaker function
\[
\Psi_{j,q}(\phi,\theta)
=
\langle \nu,\psi_1,\psi_2 \mid e^{2\phi H}e^{2i\theta Z}\mid \lambda,\psi_3,\psi_4\rangle,
\]
which the paper computes explicitly in terms of modified Bessel functions \(K_{2j-1}\), \(K_{2j}\), and \(K_{2j+1}\) [2405.09289].

After imposing the Brown–Henneaux choice \(\lambda=\nu=1/\epsilon\) with \(\epsilon\to0\), the resulting wavefunction matches the \(\mathcal N=2\) Liouville minisuperspace Hamiltonian derived from the quadratic Casimir,
\[
\Big(-\partial_\ell^2-\frac14\partial_\theta^2+e^{-\ell}
+ie^{-\ell/2-i\theta}\,\overline{\psi}_-\psi_+
+ie^{-\ell/2+i\theta}\,\psi_-\overline{\psi}_+\Big)\Psi_{k,q}(\ell,\theta)
=
E(k,q)\Psi_{k,q}(\ell,\theta),
\]
with \(E(k,q)=k^2+q^2\) and \(j=ik\). The representation-theoretic method also fixes the normalization of the fermionic multiplet and yields the density of states
\[
\rho(k,q)=\frac{k\sinh(2\pi k)}{\pi E(k,q)}.
\]

The thermofield-double or Hartle–Hawking state exhibits the strongest “fully gravitational” feature in this setting: the no-boundary condition does not permit an arbitrary fermionic multiplet component, but selects a unique linear combination,
\[
\ket{\mathrm{TFD}}
=
\sum_{q\in\mathbb Z/2}\int_0^\infty dk\;
\rho(k,q)e^{-\beta E(k,q)/2}\,
\big(i\ket{H}+\ket{L}\big),
\]
with a \(U(1)\) holonomy fixing the relative phase when present. Here the geometry, boundary condition, and fermionic structure are tied together at the level of group representation theory [2405.09289].

## 7. Operational waveform constructions for fully-precessing eccentric binaries

A more operational usage appears in gravitational-wave modeling for compact binary inspirals. The Efficient Fully-Precessing Eccentric family, EFPE, is a fast frequency-domain model for spin-precessing, eccentric compact-binary inspirals; EFPE\(_{\rm ME}\) extends it to moderate eccentricities by retaining the same precession dynamics and Fourier-transform framework but replacing the small-\(e\) amplitude expansion with a full harmonic-series representation [2402.06804].

The model includes inspiral dynamics only, with no merger or ringdown, spin-induced orbital precession, orbital eccentricity described by the time eccentricity \(e\), and orbit-averaged PN radiation reaction. The source-frame geometry is organized using Euler angles \(\theta_L\), \(\phi_z\), and \(\zeta\), and the waveform is built by twisting orbital-frame modes into an inertial frame:
\[
h^{lm}=\sum_{m'=-l}^{l} D^l_{m'm}(\phi_z,\theta_L,\zeta)\,H^{lm'}.
\]
The orbital-frame modes are decomposed into harmonics of the mean anomaly,
\[
\hat H^{lm}(t)=\sum_{j=-\infty}^{\infty}N_j^{lm}e^{-i(j\ell+m\lambda)},
\]
and the Fourier transform is computed with shifted uniform asymptotics rather than plain stationary phase, because precession introduces stationary-phase pathologies [2402.06804].

The distinction between EFPE and EFPE\(_{\rm ME}\) lies in the amplitude sector. Earlier EFPE amplitudes used a post-circular expansion truncated at about \(O(e^6)\), giving a regime of validity around \(e<0.3\) at the reference time four years before merger. EFPE\(_{\rm ME}\) instead constructs the amplitudes as full Fourier series in the mean anomaly and evaluates enough harmonics to meet a target error tolerance, remaining accurate up to \(e\le 0.8\). In the LISA band, the comparison shows no significant difference between EFPE and EFPE\(_{\rm ME}\) for \(e_0\le0.5\), but significant deviations for \(e_0>0.5\), especially for total masses below about \(10^5\,M_\odot\). The extension is, however, computationally expensive: at \(e_0=0.8\), EFPE\(_{\rm ME}\) can be about 1000 times slower than EFPE [2402.06804].

This usage is conceptually different from curvature-based detector response, null-infinity radiative data, or Hartle–Hawking/JT states. It nevertheless shares the same methodological tendency: the waveform is constructed so that the relevant gravitational structure—here simultaneous precession and eccentricity—is retained rather than frozen or linearized away. Across the literature, that is the strongest common thread linking the otherwise disparate meanings of “fully gravitational wave function.”

Source: https://www.emergentmind.com/topics/fully-gravitational-wave-function