Gravitational Trilinear Kernel in Tau-Function Gravity

• Gravitational trilinear kernel is a unique multilinear operator that encodes the second-derivative nonlinearities in stationary axisymmetric Einstein equations.
• It employs GL(2)-covariant, Hirota-type trilinear operators to unify multidimensional soliton hierarchies with gravitational geometry.
• Its formulation in the tau-function framework clarifies the algebraic structure of Tomimatsu–Sato spacetimes and related integrable gravitational metrics.

A gravitational trilinear kernel refers to a specific multilinear operator structure encoding the core nonlinearity of stationary axisymmetric gravitational field equations—specifically, the Ernst equation—in the tau-function formalism. The kernel provides a universal, GL(2)-covariant trilinear encapsulation of the second-derivative sector and underpins both the integrability properties of certain metrics (notably the Tomimatsu–Sato spacetimes) and the algebraic unification of multidimensional soliton hierarchies with gravitational geometry. Its foundational role, uniqueness, and explicit construction have been established in recent analyses of integrable structures in stationary axisymmetric general relativity and in the multilinear realization of the Yu–Toda–Fukuyama hierarchy (Fukuyama, 2 Dec 2025, Fukuyama, 16 Jan 2026).

1. Tau-Function Formalism and the Ernst Equation

The tau-function formalism introduces a projective parametrization of the complex Ernst potential $\mathcal{E}$, which encodes stationary axisymmetric vacuum gravitational metrics: $$\mathcal{E} = \frac{\tau_1}{\tau_0}, \qquad \tau_0 \neq 0$$ The Ernst equation (with flat two-dimensional base $(x, y)$) is expressed in reduced real form as

$$\mathcal{E} \, \Delta \mathcal{E} - |\nabla \mathcal{E}|^2 = 0$$

Substituting the $\tau$-ratio ansatz and clearing denominators yields a fourth-degree polynomial numerator $N(\tau_0, \tau_1)$, which decomposes canonically into cubic and quartic (gradient envelope) parts: $$N(\tau_0,\tau_1) = N_{\rm cubic}(\tau_0,\tau_1) + N_{\rm quartic}(\tau_0,\tau_1)$$ where

$$N_{\rm cubic} = \tau_1 \tau_0^2 \Delta \tau_1 - \tau_1^2 \tau_0 \Delta \tau_0$$

$$N_{\rm quartic} = -\tau_0^2 |\nabla \tau_1|^2 + \tau_1^2 |\nabla \tau_0|^2$$

Here, all second-derivative (nonlinear) terms reside exclusively in $N_{\rm cubic}$, isolating the sector of interest for integrability analyses (Fukuyama, 2 Dec 2025).

2. Definition and Algebraic Structure of the Trilinear Kernel

The trilinear kernel leverages $Z_3$-symmetric Hirota-type operators acting on function triples. Introducing the primitive third root of unity $\omega = e^{2\pi i/3}$, the trilinear derivative (in direction $x$) is: $$T_x[a, b, c] = (\partial_{x_1} + \omega \partial_{x_2} + \omega^2 \partial_{x_3})\, a(x_1) b(x_2) c(x_3) \big|_{x_1=x_2=x_3=x}$$ The gravitational trilinear kernel for the projective pair $(\tau_0,\tau_1)$ is the unique minimal GL(2)-covariant, homogeneous, and antisymmetric object constructed from these Hirota operators: $$\mathcal{Y}(\tau_0, \tau_1) = \partial_x[\tau_1 T_x(\tau_0,\tau_0,\tau_1) - \tau_0 T_x(\tau_1,\tau_1,\tau_0)] + \partial_y[\tau_1 T_y(\tau_0,\tau_0,\tau_1) - \tau_0 T_y(\tau_1,\tau_1,\tau_0)]$$ This structure ensures genuine three-way interference and projects out bilinear (solitonic) terms, emphasizing its multilinear character (Fukuyama, 2 Dec 2025, Fukuyama, 16 Jan 2026).

3. Connection to Integrable Hierarchies and the YTF Kernel

Yu–Toda–Fukuyama (YTF) established that the universal trilinear kernel in multidimensional soliton hierarchies generates all commuting flows: $$\bigl(36\,T_x^2 T_t + T_x^4 T_z^* + 8\,T_x^3 T_x^* T_z + 9\,T_z^3\bigr)\,\tau\cdot\tau\cdot\tau = 0$$ When generalized to the projective $\tau$-function formulation required by gravity, the gravitational trilinear kernel $\mathcal{Y}(\tau_0,\tau_1)$ is the unique GL(2)-covariant projection of the YTF kernel, with vanishing precisely equivalent to the Ernst equation (Fukuyama, 16 Jan 2026). Particular reductions correspond to specific integrable flows, while the full kernel encodes multidimensional soliton dynamics.

4. Tomimatsu–Sato Geometries and Degeneracy

For the $\delta=2$ Tomimatsu–Sato (TS) spacetime family, the explicit polynomials $A(\xi,\eta)$, $B(\xi,\eta)$ define

$\tau_0 = A + B, \qquad \tau_1 = A - B$

with

$\mathcal{E} = \tau_1 / \tau_0$

The cubic sector matches the trilinear YTSF kernel: $$N_{\rm cubic}(\tau_0,\tau_1) = \kappa \mathcal{Y}(\tau_0,\tau_1) + N_{\leq 1}(\tau_0,\tau_1),\qquad \kappa = -4 p^2 q^2$$ where $N_{\leq 1}$ contains only first-derivative terms. Symbolic calculation verifies that all second-derivative content is captured by $\mathcal{Y}$, demonstrating a gravitational realization of the YTSF integrable structure. In generic cases, the trilinear kernel is nonzero; TS metrics are degenerate submanifolds where three-way interference collapses, and the kernel reduces to Hirota bilinear form (Fukuyama, 2 Dec 2025, Fukuyama, 16 Jan 2026).

5. Physical Interpretation and Integrability

The quartic gradient envelope $N_{\rm quartic}$ serves as a universal first-derivative renormalization of the conformal factor, leaving integrability unaffected. The core nonlinearities of stationary axisymmetric vacuum gravity are encoded entirely in the trilinear kernel. Integrability is characterized by a dual-layer hierarchy:

• Discrete index: Bilinear Hirota/Toda-molecule hierarchy in soliton number ($\delta$)
• Continuous base: Trilinear kernel in Weyl-base coordinates ($\xi,\eta$)

The gravitational trilinear kernel embodies the minimal closed algebra governing multidimensional, non-factorizable soliton dynamics. Degenerate submanifolds correspond to solution sectors (e.g., TS, Kerr-NUT, Plebański–Demiański geometries) with enhanced algebraic constraints eliminating genuine trilinear interaction.

6. Generalizations and Theoretical Implications

For multipole index $\delta>2$, the tau-functions are constructed via the semi-infinite Toda-molecule equation. The gravitational trilinear kernel $\mathcal{Y}(\tau_0,\tau_1)$ remains universal across the TS hierarchy. Empirical strategies include constructing discrete $\tau$ towers, splitting the Ernst numerator, and verifying annihilation of second-derivative monomials by the trilinear kernel for $\delta\geq 3$. The uniqueness and necessity of the trilinear kernel in stationary gravity and multidimensional soliton hierarchies establish a unified framework linking projective geometry, multidimensional integrability, and solution classification (Fukuyama, 2 Dec 2025, Fukuyama, 16 Jan 2026).

7. Relation to Lagrangian and Field-Theoretic Trilinear Structures

Light-cone gauge formulations of the Einstein–Hilbert Lagrangian reveal a trilinear graviton interaction vertex that is totally symmetric under permutation of legs and can be derived directly from the double-copy of Yang–Mills theories. The cubic interaction kernel in field-theoretic language provides the fundamental building block for tree-level graviton amplitudes and is structurally analogous to the projective trilinear kernel seen in stationary gravity when mapped to tau-function language, although the explicit representations may differ (Beneke et al., 2021).

---

In summary, the gravitational trilinear kernel provides the minimal, unique, GL(2)-covariant multilinear operator capturing all nonlinear second-derivative structure of the stationary axisymmetric vacuum Einstein equations in tau-function coordinates. It unifies soliton integrable hierarchies in multiple dimensions with gravitational geometry, subsumes standard bilinear hierarchies as degenerate cases, and manifests explicitly in the algebraic structure of Tomimatsu–Sato and related spacetimes.