The modular Hamiltonian of a region—the logarithm of the reduced density matrix of the vacuum—is, for most regions and most theories, a non-local operator with no closed form. The known exceptions are precious. For the Rindler wedge of any Wightman theory, the Bisognano–Wichmann theorem identifies the modular Hamiltonian with \(2\pi\) times the boost generator . For a causal diamond, locality of the modular Hamiltonian is known only for conformally invariant fields: Hislop and Longo established the geometric modular action for the massless scalar , and Casini, Huerta and Myers exhibited the generator explicitly , \[\begin{equation} K\;=\; 2\pi \int_{r<R} d^3x\; f(r)\, T_{00}(x), \qquad f(r)= \frac{R^2 - r^2}{2R}, \label{eq:CHM} \end{equation}\] via the conformal map that takes the diamond to \(\mathbb{R}\times\mathbb{H}^3\) and its modular flow to a genuine Killing time-translation.
The geometric character of \(\eqref{eq:CHM}\) has acquired a structural role in quantum gravity. The crossed-product constructions of Type II von Neumann algebras—initiated for the de Sitter static patch by Chandrasekaran, Longo, Penington and Witten , following , and extended to general subregions by Jensen, Sorce and Speranza —identify the modular flow of a causal region with the proper time of an observer who carries a clock along it. The identification is physically clean precisely when the flow is geometric: a flow that acts as a time-translation is one an observer can follow. For a finite diamond and conformal matter, the Casini–Huerta–Myers (CHM) map supplies this geometric structure exactly.
The graviton breaks the pattern. The Fierz–Pauli action is not conformally invariant—among free massless gauge fields in four dimensions, the Maxwell field and the conformally coupled scalar are conformal field theories, the graviton is not—so the CHM map does not act on it, and the route that establishes \(\eqref{eq:CHM}\) is unavailable. Since the gravitational field is the one field a theory of quantum gravity cannot treat as optional, this looks like a structural obstruction: diamond modular theory, and the Type II constructions built on it, would apply to every field except the one that matters most. The prevailing expectation, shaped by the failure of the conformal route at its very first step, has accordingly been that the method does not extend to the physical spin-2 field.
This paper shows that the expectation conflates a sufficient condition with a necessary one. Conformal invariance of the action is one way to establish geometric modular flow; it is not the only way. The modular flow of a region is intrinsic to the pair (algebra, state). The right question is therefore not “does the CHM map act on the graviton action?” but
Does the gauge-invariant graviton algebra of a ball, in the Minkowski vacuum, have geometric modular flow?
We show that the answer is yes, and that every load-bearing ingredient is already in the published literature; what has been missing is their assembly.
Three results combine.
Dimensional reduction of the graviton ball algebra. Benedetti and Casini decomposed the linearized graviton in tensor spherical harmonics on the \(t=0\) ball and fixed the gauge so that the two dynamical modes for each angular momentum \((\ell, m)\) decouple, each with the dynamics of a massless-scalar spherical mode; the \(\ell = 0, 1\) towers are absent (they carry no radiative degrees of freedom), and the gauge-fixed degrees of freedom inside the sphere represent gauge-invariant operators localized in the same region. The graviton ball algebra is two scalar towers with \(\ell \geq 2\).
Block-diagonality. The Minkowski vacuum is Gaussian and rotation invariant, so its one-particle data restricted to the ball—the symplectic form together with the correlators—are block-diagonal over \((\ell, m, \text{parity})\), and by Araki’s modular theory of quasi-free states so is the modular operator. The modular flow does not mix the towers (Section 3.3).
Per-mode geometric flow without conformal invariance. The dimensionally reduced angular modes are half-line theories with a centrifugal potential; they are not conformally invariant. Nevertheless, Huerta and van der Velde showed that for an interval \((0, R)\) attached to the origin, the modular Hamiltonian of each reduced mode is local in the energy density with the parent weight \(f(r)\) of \(\eqref{eq:CHM}\), and that this local operator generates the modular flow—the surviving portion of the parent conformal symmetry suffices. Geometric modular flow, one non-conformal tower at a time.
Chaining the three: the vacuum modular flow of the gauge-invariant graviton ball algebra is the CHM conformal-Killing flow restricted to the \(\ell \geq 2\) double tower, with generator \[\begin{equation} K_{\mathrm{grav}} \;=\; 2\pi \sum_{P = \pm}\,\sum_{\ell \geq 2}\sum_{|m| \le \ell}\, \int_0^R dr\; f(r)\; T_{00}^{(\ell, m, P)}(r), \label{eq:main-result} \end{equation}\] where \(T_{00}^{(\ell,m,P)}\) is the energy density of the reduced one-dimensional mode of parity \(P\), Eq. \(\eqref{eq:Hell}\). The graviton’s diamond modular flow is geometric—for a field that is not a conformal field theory.
The mode equivalence in step 1 has an independent, elementary anchor: it is the flat-space degeneration of black-hole perturbation theory. The odd-parity (Regge–Wheeler ) and even-parity (Zerilli ) master potentials on a Schwarzschild background of mass \(M\) are distinct from each other and from the massless scalar’s potential at \(M \neq 0\); at \(M = 0\) all three collapse to the same centrifugal potential \(\ell(\ell+1)/r^2\), and the common master solutions are the Riccati–Bessel functions \(r\,j_\ell(\omega r)\)—the scalar spherical modes. We verify this degeneration symbolically for general \(\ell\) and at machine precision on a radial lattice (Section 6), where the per-mode equivalence is confirmed at the level of symplectic spectra, trace invariants, and entanglement entropies.
The assembly is validated on a case where the answer is independently known: the Maxwell field, whose gauge-invariant sphere content is two scalar towers with \(\ell \geq 1\) , is a CFT, so its geometric diamond flow follows from both the conformal route and the reduction route. The two agree (Proposition [prop:maxwell]), and the graviton then differs from Maxwell only in where its towers start. The resulting helicity ladder—two scalars \(\supset\) Maxwell (\(\ell \geq 1\)) \(\supset\) graviton (\(\ell \geq 2\))—is verified mode by mode on the lattice.
The geometric flow carries three consequences (Section 5). The local (Tolman) temperature of the diamond, \(T_{\mathrm{loc}}(r) = 1/(2\pi f(r))\), applies to the graviton towers as to conformal matter: the diamond’s thermality is helicity-blind. The crossed product of the graviton diamond algebra by its modular flow is a Type II\(_\infty\) von Neumann algebra with a canonical trace , and the geometric character of the flow places the physical observer-clock identification of on the same footing for the graviton as for conformal matter. And any centrally localized excitation of energy \(\Delta E\) accumulates modular phase at the proper-time rate \(\Delta E/\hbar\), independently of the arbitrary diamond radius \(R\)—a statement that now holds uniformly across helicity, and that bears on proposed gravitational-decoherence mechanisms in which a modular clock sets the rate (Section 5.3).
Equally important is what this paper does not claim (Section 7). The result concerns the modular flow; the entropy bookkeeping of the graviton diamond is subtler. The graviton is not globally equivalent to two conformal scalars—the universal logarithmic coefficient of its sphere entropy is \(-61/45\) , far from twice the scalar’s \(-1/90\) —because the missing \(\ell = 0, 1\) towers and the algebra’s edge structure contribute; edge modes and the center of the algebra affect entropy constants and are the subject of an active literature , but not the bulk flow. The observer clock of the crossed product remains transcribed from the de Sitter construction rather than built intrinsically for the diamond—a residue our result renders no longer graviton-specific at the flow level, though graviton-specific subtleties of the transcription itself (constraint structure, edge sector) remain possible. The trace-anomaly normalization of the crossed-product trace is not computed here. And \(\eqref{eq:main-result}\) is a mode-sum statement on the gauge-invariant algebra: we do not claim a covariant local gauge-invariant energy density for the graviton, which does not exist in the standard sense. All results are at linearized order in the gravitational field.
Every ingredient we use is published: for the reduction, for the per-mode modular Hamiltonian, for the parent geometric flow, for the crossed product, for the Maxwell reduction. The contribution of this paper is the assembly and its consequence: to our knowledge, the statement that the physical graviton’s diamond modular flow is geometric—and that the Type II construction therefore extends to it—has not been drawn before, and the prevailing expectation, shaped by the failure of the conformal route, has been the opposite. Corroborating evidence that the hyperbolic-cylinder frame captures graviton bulk physics despite non-conformality appears in the work of David and Mukherjee , who found that the Benedetti–Casini sphere entropy is reproduced by a Kaluza–Klein tower on \(S^1 \times AdS_3\), with edge modes accounting for the difference in a companion analysis . The flat-space degeneration of the Regge–Wheeler and Zerilli equations is classical ; its role as a first-principles anchor for diamond modular theory appears to be new.
Section 2 fixes conventions and states the problem. Section 3 assembles the main result and the Maxwell validation. Section 4 constructs the crossed product. Section 5 develops the consequences, Section 6 the numerical and symbolic verification, and Section 7 the honest residues and the general spin-\(s\) outlook.
Let \(D\) be the causal diamond over the \(t=0\) ball \(B_R = \{r < R\}\) in Minkowski space, \(\mathcal{A}(D)\) the von Neumann algebra of observables localized in \(D\), and \(\Omega\) the Minkowski vacuum. Tomita–Takesaki theory assigns to the pair \((\mathcal{A}(D), \Omega)\) a modular operator \(\Delta\) and the one-parameter modular automorphism group \(\sigma_s = \mathrm{Ad}\,\Delta^{is}\); the modular Hamiltonian is \(K= -\ln \Delta\). For a general region and theory, \(\sigma_s\) acts non-locally. The flow is called geometric when it acts by point transformations—when there is a vector field \(\xi\) on \(D\) such that \(\sigma_s\) moves operators along the flow of \(\xi\).
For conformally invariant fields the diamond flow is geometric : \(\xi\) is the conformal Killing vector preserving \(D\), \[\begin{equation} \xi \;=\; \frac{\pi}{R}\Big[\,\tfrac{1}{2}\big(R^2 - t^2 - r^2\big)\,\partial_t \;-\; t\,r\,\partial_r\Big], \label{eq:CKV} \end{equation}\] normalized so that \(\xi^t|_{t=0} = \pi f(r)\) with \(f(r)= (R^2 - r^2)/(2R)\); the generator on the \(t = 0\) slice is the CHM expression \(\eqref{eq:CHM}\) , whose flow has lapse \(2\pi f(r)\), so \(\sigma_s\) is the flow of \(\eqref{eq:CKV}\) at parameter \(2s\). Two features of this structure matter below; we work in units \(\hbar = c = k_B = 1\) except where stated.
First, the weight \(f(r)\) vanishes at the edge \(r = R\) and peaks at the center, \(f(0) = R/2\); the associated local temperature \[\begin{equation} T_{\mathrm{loc}}(r) \;=\; \frac{1}{2\pi f(r)} \label{eq:Tloc} \end{equation}\] is the temperature registered by a static detector that follows the modular flow at radius \(r\).
Second, the route to \(\eqref{eq:CHM}\) runs through conformal symmetry, and its geometry is explicit : the Weyl rescaling \(\hat g = \Omega^{-2}\eta\) with \[\begin{equation} \Omega^2 \;=\; \frac{[(R+t)^2 - r^2]\,[(R-t)^2 - r^2]}{R^4} \label{eq:weyl} \end{equation}\] maps \(D\) to the static cylinder \(\mathbb{R}\times\mathbb{H}^3\); the \(t = 0\) ball becomes the hyperbolic space \(\mathbb{H}^3\) of curvature radius \(R/2\), with the diamond boundary \(r = R\) sent to its infinity; and \(\xi\) becomes an exact Killing vector of constant norm, \(|\xi|^2_{\hat g} = -\pi^2 R^2/4\)—a globally static time-translation with no interior horizon. The Tolman relation \(T_{\mathrm{loc}}(r)\,\Omega(r) = \mathrm{const}\) exhibits the diamond temperature profile \(\eqref{eq:Tloc}\) as the blueshifted image of the uniform cylinder temperature. For a conformal field the algebra and vacuum transport along this map, and the diamond inherits the cylinder’s thermal structure.
The transport in the last step requires conformal covariance of the field. The linearized graviton does not have it: the Fierz–Pauli action is not Weyl invariant, the graviton is not a CFT, and the CHM map simply does not act on the theory. Every step of the chain “conformal invariance \(\Rightarrow\) cylinder frame \(\Rightarrow\) Killing flow \(\Rightarrow\) geometric modular Hamiltonian” fails at the first link.
If geometric modular flow required conformal invariance, that would be the end of the story, and the modular theory of the diamond would not apply to the gravitational field. The Type II crossed-product constructions , which need the flow to be a time-translation an observer can follow, would inherit the same restriction. It is this expectation—natural, given the failure of the conformal route at its first step—that the present paper revisits.
The key observation is that the implication only runs one way. The modular flow of \((\mathcal{A}, \Omega)\) is determined by the algebra and the state; conformal invariance is a tool for computing it, not a precondition for its geometric character. There is a second route to the diamond modular Hamiltonian that never invokes conformal invariance of the full theory: dimensional reduction.
A massless scalar in four dimensions, decomposed in spherical harmonics, is an infinite tower of half-line theories: for each \((\ell, m)\) the rescaled radial field \(\psi_{\ell m}(r) = r\,\phi_{\ell m}(r)\), regular at the origin (\(\psi \sim r^{\ell+1}\)), has Hamiltonian \(H_\ell = \int_0^\infty dr\; T_{00}^{(\ell)}(r)\) with reduced energy density \[\begin{equation} T_{00}^{(\ell)}(r) \;=\; \frac{1}{2}\Big[\pi_{\ell m}^2 + (\partial_r \psi_{\ell m})^2 + \frac{\ell(\ell+1)}{r^2}\,\psi_{\ell m}^2\Big]. \label{eq:Hell} \end{equation}\] This explicit density is the object whose locality the whole construction turns on. The centrifugal potential breaks one-dimensional conformal invariance for every \(\ell \geq 1\): the reduced theories are not CFTs. Nevertheless, the ball algebra of the parent scalar is the direct sum of the interval algebras of the modes, the vacuum reduces mode by mode, and the parent modular Hamiltonian \(\eqref{eq:CHM}\) therefore decomposes into mode contributions. Huerta and van der Velde analyzed the reduced theories directly and showed that for the interval \((0, R)\) attached to the origin each mode’s modular Hamiltonian is local, \[\begin{equation} K_\ell \;=\; 2\pi \int_0^R dr\; f(r)\; T_{00}^{(\ell)}(r), \label{eq:Kell} \end{equation}\] with the same weight \(f(r)\) as the parent, and—the nontrivial part—that \(\eqref{eq:Kell}\) generates the modular flow of the reduced mode. The symmetries that survive the reduction (the subgroup of the parent conformal group commuting with the angular decomposition) are enough; full conformal invariance of the reduced theory is not needed.
This is the mechanism we exploit for the graviton: if the gauge-invariant graviton ball algebra decomposes into towers whose reduced data—radial Hamiltonian, origin boundary condition, and vacuum—coincide with those of \(\eqref{eq:Hell}\), then \(\eqref{eq:Kell}\) applies tower by tower, and the graviton’s diamond modular flow is geometric with no appeal to conformal invariance of the graviton action. The next section establishes exactly this decomposition.
Benedetti and Casini analyzed the linearized graviton \(h_{\mu\nu}\) on the \(t = 0\) ball: decompose \(h_{\mu\nu}\) in tensor spherical harmonics (Regge–Wheeler conventions ) and fix the gauge so that, for each angular momentum \((\ell, m)\), the two dynamical (radiative) modes decouple, one per parity sector. Their results, in the form we need:
For each \((\ell, m)\) with \(\ell \geq 2\) there are exactly two dynamical modes (one even-parity, one odd-parity), and each has the dynamics of a massless-scalar spherical mode: the reduced radial Hamiltonians are those of the dimensionally reduced scalar, Eq. \(\eqref{eq:Hell}\), at the same \(\ell\), with the same regularity condition \(\psi \sim r^{\ell+1}\) at the origin.
The \(\ell = 0\) and \(\ell = 1\) sectors contain no radiative modes. This is the angular-momentum form of the Arnowitt–Deser–Misner (ADM) count : of the ten components of \(h_{\mu\nu}\), four are gauge and four are fixed by constraints, leaving the two helicities, which appear first at \(\ell = 2\). The low-\(\ell\) sectors are constraint/gauge data (mass, momentum, center of mass, angular momentum), not local degrees of freedom.
The gauge-fixed field degrees of freedom inside the sphere represent gauge-invariant operators of the theory localized in the same region: the reduction is a statement about the local algebra, not merely about mode functions.
Consequently the sphere entanglement entropy of the graviton equals that of a pair of massless scalars with the \(\ell = 0, 1\) contributions removed, with universal logarithmic coefficient \(-61/45\).
Of these, (R3) is the load-bearing statement for modular theory: it identifies the gauge-invariant graviton ball algebra as \[\begin{equation} \mathcal{A}_{\mathrm{grav}}(B_R)\;\cong\; \bigotimes_{P = \pm}\;\bigotimes_{\ell \geq 2,\, |m| \leq \ell} \mathcal{A}_\ell^{(m, P)}\big((0,R)\big), \label{eq:algebra-decomp} \end{equation}\] where each factor is the interval algebra of a reduced scalar mode \(\eqref{eq:Hell}\), and the Minkowski vacuum restricts to the corresponding mode vacuum factor by factor. (The algebra tensor-factorizes over the blocks; correspondingly, the one-particle structures introduced in Section 3.3 decompose as a direct sum.)
Statement (R1) can look like an accident of gauge fixing. It is not; it is the flat-space limit of a classical fact of black-hole perturbation theory, and exhibiting it this way makes the per-mode equivalence checkable by elementary means.
On a Schwarzschild background of mass \(M\), the two parities of the graviton obey master equations \(-\psi'' + V(r)\psi = \omega^2 \psi\) (primes denoting tortoise-coordinate derivatives) with the Regge–Wheeler and Zerilli potentials \[\begin{align} V_{\mathrm{RW}}&= \Big(1 - \frac{2M}{r}\Big)\left[\frac{\ell(\ell+1)}{r^2} - \frac{6M}{r^3}\right], \label{eq:VRW}\\[2pt] V_{\mathrm{Z}}&= \Big(1 - \frac{2M}{r}\Big)\, \frac{2\lambda^2(\lambda+1)r^3 + 6\lambda^2 M r^2 + 18\lambda M^2 r + 18 M^3}{r^3\,(\lambda r + 3M)^2}, \qquad \lambda = \tfrac{(\ell-1)(\ell+2)}{2}, \label{eq:VZ} \end{align}\] while a massless scalar obeys the same equation with \(V_{\mathrm{s}}= (1 - 2M/r)\,[\ell(\ell+1)/r^2 + 2M/r^3]\). At \(M \neq 0\) the three potentials are pairwise distinct (the relation between \(\eqref{eq:VRW}\) and \(\eqref{eq:VZ}\) is isospectrality via a Darboux transformation , not equality). At \(M = 0\) they degenerate: \[\begin{equation} V_{\mathrm{RW}}\big|_{M=0} \;=\; V_{\mathrm{Z}}\big|_{M=0} \;=\; V_{\mathrm{s}}\big|_{M=0} \;=\; \frac{\ell(\ell+1)}{r^2}, \label{eq:degeneration} \end{equation}\] where the Zerilli case uses the identity \((\ell-1)(\ell+2) + 2 = \ell(\ell+1)\); the tortoise coordinate degenerates to \(r\); and the common master solutions are the Riccati–Bessel functions \(\psi = r\,j_\ell(\omega r)\), regular at the origin (\(\psi \sim r^{\ell+1}\))—precisely the scalar spherical modes of \(\eqref{eq:Hell}\), with the same origin boundary condition. Both graviton parities, mode by mode in \((\ell, m)\), are scalar modes in flat space. We verify \(\eqref{eq:degeneration}\) symbolically for general \(\ell\), and the resulting lattice couplings agree to the last floating-point bit while being pairwise distinct at \(M \neq 0\) (Section 6)—the degeneration is exact and the comparison is not vacuous.
Let \(\omega\) be a quasi-free (Gaussian) state of a free bosonic field, invariant under spatial rotations and parity, and let the region algebra decompose as in \(\eqref{eq:algebra-decomp}\) into \((\ell, m, P)\) blocks generated by variables localized in the region. Then the modular operator of \((\mathcal{A}(B_R), \omega)\) is the second quantization of a direct sum over the blocks: the modular flow acts within each tower.
Argument. For a free Bose field, the local von Neumann algebra is generated by Weyl operators over the real symplectic space \((S_A, \sigma_A)\) of initial data supported in \(B_R\), and the quasi-free state is specified by a covariance on \(S_A\)—concretely, by the pair of restricted correlators \((X_A, P_A)\) of Appendix 10. By Araki’s theory of quasi-free states and the Tomita–Takesaki structure of free-field algebras , the modular operator of the pair (Weyl algebra of \((S_A, \sigma_A)\), quasi-free state) is the second quantization of a one-particle modular operator built from \((\sigma_A, X_A, P_A)\) alone. Now use the symmetry twice: the ball is rotation and parity invariant, so the restriction to \(B_R\) commutes with the \(SO(3) \times\) parity action; and the symplectic form, the vacuum correlators, and hence the full one-particle data are invariant under that action. Both \(\sigma_A\) and \((X_A, P_A)\) are therefore simultaneously block-diagonal over the \((\ell, m, P)\) decomposition, the one-particle modular operator decomposes as a direct sum of the blocks’ modular operators, and its second quantization factorizes over the towers. The one hypothesis doing real work is that the block variables are themselves localized in the ball—which is exactly statement (R3) of the Benedetti–Casini reduction. ◻
The vacuum modular flow of the gauge-invariant linearized-graviton algebra \(\mathcal{A}_{\mathrm{grav}}(B_R)\) of a ball \(B_R\) in Minkowski space is geometric: it acts within each \((\ell, m, P)\) block as the modular flow of the corresponding reduced scalar mode, which is generated by the local operator \(\eqref{eq:Kell}\); equivalently, the full generator is the mode sum \[\begin{equation} K_{\mathrm{grav}} \;=\; 2\pi \sum_{P = \pm}\sum_{\ell \geq 2}\sum_{|m| \leq \ell} \int_0^R dr\; f(r)\; T_{00}^{(\ell, m, P)}(r), \qquad f(r)= \frac{R^2 - r^2}{2R}, \label{eq:thm-main} \end{equation}\] with \(T_{00}^{(\ell,m,P)}\) the reduced energy density \(\eqref{eq:Hell}\): the Hislop–Longo/CHM conformal-Killing flow restricted to the \(\ell \geq 2\) double tower.
Assembly. By the Benedetti–Casini reduction (R1)–(R3), \(\mathcal{A}_{\mathrm{grav}}(B_R)\) decomposes as \(\eqref{eq:algebra-decomp}\) with the vacuum reducing mode by mode to the scalar mode vacua; the flat master-equation degeneration \(\eqref{eq:degeneration}\) supplies the per-mode identity independently. By Lemma [lem:block] the modular operator is block-diagonal over the towers. On each block, the (algebra, state) pair coincides with that of a reduced scalar mode on \((0, R)\)—the same radial Hamiltonian \(\eqref{eq:Hell}\), the same regularity condition \(\psi \sim r^{\ell+1}\) at the origin, and the same vacuum—and this is precisely the setting of the Huerta–van der Velde theorem , which holds for every angular mode \(\ell\) of the reduced scalar and hence in particular for \(\ell \geq 2\): the modular Hamiltonian is the local expression \(\eqref{eq:Kell}\) and generates the flow. Since the modular operator is intrinsic to the (algebra, state) pair, each graviton block inherits it. Summing blocks gives \(\eqref{eq:thm-main}\). ◻
Three remarks delimit the claim. First, the theorem is a statement at linearized order; interactions are outside its scope, as everywhere in free-field modular theory. Second, Eq. \(\eqref{eq:thm-main}\) is a mode sum over the gauge-invariant towers. We do not claim it can be rewritten as \(2\pi\int d^3x\, f(r)\, T^{\mathrm{grav}}_{00}(x)\) for a covariant, gauge-invariant local energy density: no such density exists for the graviton in the standard sense, and the mode-sum form is what the algebra supports—and all that the applications below require. Third, the statement concerns the bulk tower algebra \(\eqref{eq:algebra-decomp}\); the boundary (edge-mode/center) sector, which contributes to entropy constants, is discussed in Section 7.
The assembly can be tested against a case where the answer is independently known. The Maxwell field is a conformal field theory in four dimensions, so its diamond modular flow is geometric by the classic conformal route . But Maxwell also admits a Benedetti–Casini-type reduction: Casini and Huerta showed that the gauge-invariant Maxwell content of a sphere is two towers of massless-scalar spherical modes with the \(\ell = 0\) mode removed (towers \(\ell \geq 1\); one per helicity), with sphere entropy equal to two scalars minus their \(\ell = 0\) contributions and universal logarithmic coefficient \(-16/45\).
Applied to the Maxwell field, the assembly of Theorem [thm:main]—reduction , Lemma [lem:block], and the per-mode theorem —yields geometric modular flow with generator \(2\pi \sum \int f\, T_{00}^{(\ell, m, P)}\) over the \(\ell \geq 1\) double tower. This agrees with the result of the conformal route, which is available for Maxwell because it is a CFT.
The two routes are logically independent: the conformal route transports the wedge result through the CHM map and never decomposes in modes; the reduction route never invokes conformal invariance of the action. Their agreement on Maxwell is therefore a genuine validation of the reduction route on the one gauge field where both apply—and the graviton then differs from Maxwell only in where its towers start (\(\ell \geq 2\) instead of \(\ell \geq 1\), two absent low-\(\ell\) sectors instead of one). The resulting helicity ladder \[\begin{equation} \text{two scalars } (\ell \geq 0) \;\supset\; \text{Maxwell } (\ell \geq 1) \;\supset\; \text{graviton } (\ell \geq 2) \label{eq:ladder} \end{equation}\] is verified on the lattice at machine precision in Section 6: the per-mode data are identical along the ladder, and the successive entropy deficits are exactly the dropped towers (\(2 S_0\) from scalars to Maxwell, \(6 S_1\) from Maxwell to graviton, per mode multiplicities).
The local algebra of a causal diamond in a QFT vacuum is a Type III\(_1\) factor : it admits no trace, no density matrices, and no canonical normalization of entropy. The crossed-product construction repairs this. For a Type III\(_1\) factor, the crossed product by its own modular automorphism group, \[\begin{equation} \widehat{\mathcal{A}} \;=\; \mathcal{A}_{\mathrm{grav}}(B_R)\rtimes_{\sigma} \mathbb{R}, \label{eq:crossed} \end{equation}\] is a Type II\(_\infty\) von Neumann algebra, with a faithful semifinite trace constructed from the dual weight and the added generator; this is Takesaki duality in the form now standard in the quantum-gravity literature . The step is geometry-blind: it requires only the modular flow \(\sigma\), which every (algebra, state) pair possesses, and in particular it never needed the CHM map. The existence of a Type II algebra with a trace for the graviton diamond was therefore never in doubt.
What the conformal route supplied—and what was thought to fail for the graviton—is the physical content layered on top of \(\eqref{eq:crossed}\) by Chandrasekaran, Longo, Penington and Witten and by Jensen, Sorce and Speranza : the crossed-product generator is identified with the energy of a gravitating observer, the modular flow with the observer’s proper time, and the renormalized trace with a finite von Neumann (generalized) entropy. This identification is clean precisely when the modular flow is a geometric time-translation that an observer can carry a clock along. For the de Sitter static patch the flow is the boost Killing field and the identification is immediate; for the flat diamond and conformal matter, the CHM map turns the conformal-Killing flow into the constant-norm Killing time of \(\mathbb{R}\times\mathbb{H}^3\), and the construction transcribes.
Theorem [thm:main] supplies exactly the missing input for the graviton: its diamond modular flow is geometric—the same conformal-Killing flow \(\eqref{eq:CKV}\), acting tower by tower. The geometry of the flow (the Weyl rescaling \(\eqref{eq:weyl}\) to \(\mathbb{R}\times\mathbb{H}^3\), the constant-norm Killing character of \(\xi\) there, the Tolman relation between the diamond and cylinder temperatures) is a statement about the flow, not about the matter content, and applies unchanged. The observer-clock identification of therefore carries over to the graviton diamond with the same standing it has for conformal matter.
Two clarifications keep the bookkeeping honest. First, “the same standing” includes the same caveat: for finite diamonds the CLPW observer is transcribed from the de Sitter construction—its bounded-below Hamiltonian and clock are carried over by structural analogy with the Killing case—rather than constructed intrinsically. Our result removes the flow-level obstruction that was specific to the graviton; whether the transcription itself harbors further graviton-specific subtleties (the constraint structure of the low-\(\ell\) sector, the edge modes) is left open. Second, the trace normalization (equivalently the additive constant of the generalized entropy) involves the trace anomaly of the field content; for the graviton this constant is tied to its anomaly coefficients and is not computed here.
That the cylinder frame captures graviton bulk physics despite the absence of conformal invariance is corroborated independently. David and Mukherjee computed the entanglement entropy of \(D=4\) gravitons across a sphere from partition functions on the hyperbolic cylinder—a Kaluza–Klein tower of transverse-traceless massive spin-2 fields on \(S^1 \times AdS_3\)—and found agreement with the direct Benedetti–Casini mode computation for the bulk contribution; the difference is carried by edge modes, analyzed in (see also for the Rindler edge sector). In our language: the graviton’s tower algebra transports to the cylinder frame mode by mode, exactly as Theorem [thm:main] asserts, while the edge sector—which does not affect the flow on the tower algebra—accounts for the entropy bookkeeping that the naive cylinder computation misses.
Because each graviton tower carries the same reduced (algebra, state) pair as the corresponding scalar mode, the local temperature seen by the modular flow, \[\begin{equation} T_{\mathrm{loc}}(r) \;=\; \frac{1}{2\pi f(r)} \;=\; \frac{R}{\pi (R^2 - r^2)}, \qquad T_{\mathrm{loc}}(0) = \frac{1}{\pi R}, \label{eq:Tloc-grav} \end{equation}\] applies verbatim to graviton excitations. A finite causal diamond is a thermal object in a precise modular sense, and that thermality does not distinguish helicity 0 from helicity 1 or 2: a graviton wavepacket near the center of a diamond of radius \(R\) experiences the same local temperature \(1/(\pi R)\) as a scalar or photon excitation. The Tolman product \(T_{\mathrm{loc}}(r)\,\Omega(r)\) is constant across the diamond, with \(\Omega\) the Weyl factor \(\eqref{eq:weyl}\) of the cylinder map—the diamond temperature profile is the blueshifted image of the uniform cylinder temperature, now for the physical spin-2 field as well.
Consider an excitation of energy \(\Delta E\) localized near the center of the diamond, \(r \ll R\), in either a matter tower or a graviton tower (units \(\hbar = c = 1\) in the intermediate steps). Restricted to centrally localized energy, the generator \(\eqref{eq:thm-main}\) reads \(K\approx 2\pi f(0)\, H_{\mathrm{loc}} = \pi R\, H_{\mathrm{loc}}\). Two consequences follow at once: the excitation shifts the (dimensionless) modular Hamiltonian by \(\Delta K= \pi R\, \Delta E\); and conjugation by \(e^{-i Ks}\) acts on centrally localized operators as a time translation with lapse \(2\pi f(0)\), so proper time advances at \(d\tau/ds = 2\pi f(0) = \pi R\) per unit modular parameter \(s\). The physical phase rate is their ratio: \[\begin{equation} \text{(modular phase rate in proper time)} \;=\; \frac{\Delta K}{d\tau/ds} \;=\; \Delta E \;\;\xrightarrow{\;\text{restoring units}\;}\;\; \frac{\Delta E}{\hbar}, \label{eq:phase-rate} \end{equation}\] independent of the arbitrary diamond radius \(R\). Both the modular gap and the modular-to-proper-time conversion scale as \(\pi R\); the observable rate does not. The diamond one draws around a physical system leaves no imprint on the local clock—an infrared-safety statement—and, by Theorem [thm:main], the clock that ticks in \(\eqref{eq:phase-rate}\) is the same clock for matter, photons, and gravitons. The modular time of a causal diamond is universal across helicity.
A long-standing hypothesis, going back to Diósi and Penrose , holds that a spatial superposition of a mass \(M\) with branch separation \(d\) loses coherence at a rate set by the gravitational self-energy of the difference between its branches, \(\Gamma \sim \EG/\hbar\) with \(\EG = GM^2/d\) for point-like branches; at the benchmark \(M = 1~\mu\mathrm{g}\), \(d = 1~\mathrm{mm}\) this gives \(\EG/\hbar \approx 6.3 \times 10^{8}\,\mathrm{s}^{-1}\), a decoherence time of \(1.6\) ns. One program pursuing a microphysical derivation of this rate rests on the identification of diamond modular flow with physical time: the relative modular phase between branches, evaluated through Eq. \(\eqref{eq:phase-rate}\) with \(\Delta E = \EG\), reproduces the Diósi–Penrose rate up to an \(O(1)\) coefficient. We do not re-derive or evaluate that program here; we record what the present result contributes to it.
The contribution concerns the propagating graviton sector. In such a derivation the interaction energy \(\EG\) resides in the constrained (non-radiative, \(\ell = 0, 1\)) sector, fixed by the matter configuration, and for static sources the transverse-traceless towers remain in their vacuum. But the total diamond algebra contains the TT towers as well, and if their modular flow had been non-geometric, the identification “modular flow \(=\) physical time” for the total algebra—matter, constraint sector, and radiation—would have carried an unquantified defect. Theorem [thm:main] closes this specific failure mode: the TT towers flow geometrically with the same weight \(f(r)\) as everything else, so the phase rate \(\eqref{eq:phase-rate}\) holds uniformly across the diamond’s field content; and since for static sources both branches carry the same TT vacuum, we expect no additional modular gap from the radiative sector at leading order in the source velocities. (The latter statement is the static limit, not a theorem; time-dependent sources radiate and would contribute at higher order.)
The epistemic status deserves emphasis: nothing here establishes the modular-time identification itself, which remains a hypothesis at the full operator level within that program’s own accounting ; the present result removes one specific, previously unbounded failure mode (a non-geometric radiative sector) and leaves the remaining open items—the intrinsic observer construction and the trace normalization—common to all field content.
It is worth displaying explicitly what the per-mode equivalence does not imply, because the distinction is exactly where edge physics lives. At any fixed lattice regularization (Section 6), the graviton partial entropy over the towers satisfies the finite-sum identity \[\begin{equation} S_{\mathrm{grav}}(n) \;=\; 2 \sum_{\ell \geq 2} (2\ell + 1)\, S_\ell^{\mathrm{scalar}}(n), \label{eq:tower-sum} \end{equation}\] which is not twice the full scalar entropy: the deficit is the \(\ell = 0, 1\) contribution, \(2\,[\,S_0 + 3 S_1\,] > 0\). The helicity ladder \(\eqref{eq:ladder}\) organizes this bookkeeping: Maxwell drops only the \(\ell = 0\) pair (deficit \(2 S_0\)), the graviton additionally drops the \(\ell = 1\) towers (further deficit \(6 S_1\)). In the continuum these deficits, together with the boundary/edge sector, are why the universal logarithmic coefficients differ from the naive two-scalar value \(2 \times (-1/90) = -1/45\) : \(-16/45\) for Maxwell and \(-61/45\) for the graviton . The flow statement of Theorem [thm:main] lives on the towers and is exact; the entropy bookkeeping of the full diamond involves, in addition, the absent low-\(\ell\) (constraint) sectors and the edge modes, and is the subject of Section 7.
All quantitative claims in this paper are machine-verified: the master-potential degeneration by exact computer algebra (ten symbolic identities, residuals identically zero), and the per-mode equivalence and its consequences by a lattice computation comprising 42 numerical checks. The code is available from the author and will be provided as ancillary files.
With \(V_{\mathrm{RW}}\), \(V_{\mathrm{Z}}\), \(V_{\mathrm{s}}\) as in Section 3.2, computer algebra verifies, for symbolic \(\ell\): (i) \(V_{\mathrm{RW}}|_{M=0} = V_{\mathrm{Z}}|_{M=0} = V_{\mathrm{s}}|_{M=0} = \ell(\ell+1)/r^2\) exactly, the Zerilli case reducing via \((\ell-1)(\ell+2)+2 = \ell(\ell+1)\); (ii) the three potentials are pairwise distinct at a generic point with \(M \neq 0\); (iii) the tortoise measure degenerates, \(dr_*/dr \to 1\); and (iv) \(\psi = r\,j_\ell(\omega r)\) solves the common flat master equation for \(\ell = 2, 3, 4\). All ten identities hold as rational expressions, not merely to numerical tolerance.
The radial lattice method of Srednicki (reviewed in ) discretizes each reduced mode on \(N\) sites and computes the reduced state on the inner \(n\) sites from ground-state correlators; the entanglement data are the symplectic eigenvalues \(\nu_k\) of the restricted covariance (Appendix 10). We run two discretizations: Srednicki’s, for contact with the literature, and a direct half-line discretization of the master equation, in which the scalar and the two graviton parities are built from their own potential formulas \(V_{\mathrm{s}}\), \(V_{\mathrm{RW}}|_{M\to 0}\), \(V_{\mathrm{Z}}|_{M\to 0}\)—so that the per-mode comparison tests the degeneration \(\eqref{eq:degeneration}\), not a shared code path.
The verified statements:
Baseline. The scalar tower reproduces Srednicki’s area law: \(S(n)/(n + \tfrac{1}{2})^2 = 0.2948\) and \(0.2952\) at \(n = 10, 20\) (\(N = 60\)), consistent with the classic coefficient \(\approx 0.295\) ; all symplectic eigenvalues respect the Heisenberg bound \(\nu_k \geq \tfrac12\); per-mode entropies decrease monotonically with \(\ell\).
Per-mode identity (the crux). For \(\ell = 2, \dots, 6\) and both parities, the lattice couplings built from the graviton master potentials at \(M \to 0\) agree with the scalar’s to the last floating-point bit (\(\max|\Delta K| = 0\)); the trace invariants \(\mathrm{Tr}\,[(X_A P_A)^k]\), \(k = 1,\dots,4\), agree to better than \(10^{-15}\) relative; sorted symplectic spectra and per-mode entanglement entropies agree to better than \(10^{-12}\).
Non-vacuousness. At \(M \neq 0\) (lattice units \(M = 0.1\)) the Regge–Wheeler, Zerilli, and scalar couplings are pairwise distinct at the \(10^{-2}\) level: the \(M \to 0\) degeneration is a real limit of genuinely different theories, not a definitional identity.
Mode bookkeeping. Two helicities per \((\ell, m)\) for \(\ell \geq 2\) and none for \(\ell = 0, 1\): the tower count obeys \(\sum_{\ell=2}^{L} 2(2\ell+1) = 2[(L+1)^2 - 4]\) (e.g. \(234\) at \(L = 10\)), the ADM count \(10 - 4 - 4 = 2\), and the two-full-scalars-minus-eight-modes identity.
Tower sums and the deficit. The finite-sum identity \(\eqref{eq:tower-sum}\) holds to \(10^{-12}\) relative with each side computed through its own potential formulas; the deficit from two full scalar towers equals the \(\ell = 0, 1\) contribution and is positive; the per-mode identity is stable across lattice sizes \(N = 40, 60, 100\).
The helicity ladder (Maxwell benchmark). The mode counts obey \(\sum_{\ell=1}^{L} 2(2\ell+1) = 2[(L+1)^2 - 1]\) for Maxwell, with ladder deficits of exactly \(2\) modes (the \(\ell = 0\) pair, scalars \(\to\) Maxwell) and \(6\) modes (the \(\ell = 1\) towers, Maxwell \(\to\) graviton) per \(\ell_{\max}\); the entropy ladder \(S_{2\,\mathrm{scalars}} > S_{\mathrm{Maxwell}} > S_{\mathrm{grav}}\) holds with the successive differences equal to \(2 S_0\) and \(6 S_1\) as finite-sum identities at \(10^{-12}\) relative—the lattice form of the nesting \(\eqref{eq:ladder}\).
The lattice locks the algebraic content of Theorem [thm:main] and Proposition [prop:maxwell]: the per-mode (algebra, state) equivalence, in the strong form of bitwise-equal couplings and machine-precision spectral data, plus the mode bookkeeping and tower sums along the helicity ladder. The geometricity of each tower’s modular flow—the statement that \(\eqref{eq:Kell}\) generates it—is the Huerta–van der Velde theorem , an analytic result we apply, not re-prove; direct lattice reconstruction of the modular kernel is known to converge slowly and is not used as evidence here. The division of labor is deliberate: machine precision where the new content lies (the equivalence), published proof where the literature already carries the load (the flow).
At linearized order, on the gauge-invariant tower algebra of a ball: the graviton’s vacuum modular flow is geometric (Theorem [thm:main]), and the same assembly reproduces the independently known Maxwell case (Proposition [prop:maxwell]); the Type II\(_\infty\) crossed product with canonical trace exists and its physical observer-clock reading has the same standing as for conformal matter (Section 4); the diamond’s local temperature and IR-safe central phase rate \(\Delta E/\hbar\) are helicity-blind (Section 5). The inputs are published theorems ; the assembly, the flat-space master-equation anchor, and the consequences are the contribution here. The moral we would draw for the broader program of subregion algebras in gravity is simple: conformal invariance of the action was never the right hypothesis for geometric diamond modular flow—mode-wise reducibility to theories with known flow is enough, and the physical graviton has it.
The assembly uses three inputs, of which only one is field-specific: a decomposition of the gauge-invariant ball content into scalar towers, localized in the ball. This suggests a general statement.
Let a free massless field on Minkowski space admit a Benedetti–Casini-type reduction: a decomposition of its gauge-invariant ball content into towers of massless-scalar spherical modes (any starting \(\ell\), any multiplicities), localized in the ball, with the vacuum reducing mode by mode. Then its vacuum diamond modular flow is geometric—the CHM conformal-Killing flow restricted to its towers—by Lemma [lem:block] and the per-mode theorem of .
The hypothesis is established for \(s = 0\) (trivially), \(s = 1\) , and \(s = 2\) , with towers starting at \(\ell = s\)—the angular-momentum threshold for radiative multipoles of a helicity-\(s\) field. For higher spins and \(p\)-forms the corresponding tower structure underlies the hyperbolic-cylinder entropy computations of David and Mukherjee , and we expect the reduction to hold with towers at \(\ell \geq s\), but a flat-ball derivation at the algebra level (the analogue of statement (R3)) has to our knowledge not been written down for \(s \geq 3\); the proposition is stated conditionally for this reason. If the hypothesis holds generally, geometric modular flow in a diamond is the rule for free massless field content, not a peculiarity of CFTs—the role of conformal invariance was only ever to make the computation easy.
Edge modes and the center. The decomposition \(\eqref{eq:algebra-decomp}\) is the bulk tower algebra. Gauge theories and gravity assign a sphere additional boundary data—edge modes, or a center of the algebra—which contribute to entropy constants and to the universal logarithmic coefficient, and for the graviton are actively studied . Our flow statement is for the algebra modulo its center; the edge sector does not carry a flow obstruction, but its contribution to the crossed-product trace normalization is open.
The observer is transcribed, not constructed. For finite diamonds, the bounded-below observer Hamiltonian and clock of are carried over from the de Sitter Killing case by structural analogy. This was the state of the art for conformal matter before our result; our result removes the flow-level graviton-specific obstruction, while graviton-specific subtleties of the transcription itself (constraints, edge sector) remain possible and open.
Trace-anomaly normalization. The additive normalization of the crossed-product trace (equivalently of the generalized entropy) involves the anomaly coefficients of the field content and is not computed here. For the graviton this interacts with the edge sector; the two open items are likely linked.
No covariant local density. Equation \(\eqref{eq:thm-main}\) is a mode sum. A covariant, gauge-invariant local graviton energy density does not exist in the standard sense, and we make no claim to one; applications that would require it (rather than the mode sum) are outside the result.
Linearized order. Everything here is free-field modular theory. Graviton self-interaction enters at the next order in the gravitational coupling and is not addressed.
For decoherence programs. The result removes the radiative-sector failure mode of modular-clock decoherence mechanisms (Section 5.3) but does not by itself establish the modular-time identification on which such mechanisms rest; that identification remains a hypothesis at the full operator level .
Three directions seem ripe. First, the intrinsic observer: with the flow now known to be geometric for all the field content of a tabletop diamond—scalar, Maxwell, graviton—the remaining work in promoting the transcribed CLPW clock to a construction is the same problem for matter and gravity, and can be attacked once rather than per species. Second, the edge sector: the split between a helicity-blind bulk flow and a field-specific edge/anomaly bookkeeping suggests that the crossed-product trace normalization—and through it the additive entropy constant—may be computable by matching the cylinder partition functions of against the canonical Takesaki trace, a calculation that would also bear on the \(O(1)\) coefficient in modular-clock decoherence rates. Third, and most simply, the hypothesis of Proposition [prop:spin-s] deserves a direct flat-ball proof for \(s \geq 3\): if it holds, the modular geometry of finite regions is universal across the massless spectrum.
The modular flow of a causal diamond was known to be geometric for conformally invariant fields and expected to be inaccessible for the one field that gravitational physics cannot do without. We have shown that the expectation rested on mistaking a computational route for a physical requirement. The gauge-invariant content of the linearized graviton in a ball is two towers of massless-scalar spherical modes with \(\ell \geq 2\)—a published reduction with an elementary anchor in the flat-space degeneration of the Regge–Wheeler and Zerilli equations, verified here symbolically and at machine precision on the lattice. Rotational invariance keeps the vacuum modular operator block-diagonal over the towers, and each tower, though manifestly non-conformal, inherits a local modular Hamiltonian with the conformal-Killing weight \(f(r)\) that generates its flow . The graviton’s diamond modular flow is therefore geometric, its crossed product is Type II\(_\infty\) with a canonical trace, and the modular clock of a finite diamond—including its IR-safe central phase rate \(\Delta E/\hbar\)—is blind to helicity. The same assembly, run on the Maxwell field where the conformal route is also available, reproduces the known answer: the validation and the generalization come from the same mechanism, and the helicity ladder two scalars \(\supset\) Maxwell \(\supset\) graviton is exact mode by mode.
What remains open is stated plainly: the edge-mode/center sector and the trace-anomaly normalization of the crossed-product trace; the intrinsic (rather than transcribed) construction of the diamond observer; the flat-ball reduction hypothesis for spins \(s \geq 3\); and everything beyond linearized order. None of these is a graviton-specific obstruction at the level of the flow any longer. That, we think, is the structural point: the modular theory of finite regions, and the Type II algebras built from it, extend to the physical gravitational field, and the problems that remain are the problems of the subject as a whole.
For completeness we record the master equations and the degeneration used in Section 3.2. Perturbations of a Schwarzschild background of mass \(M\) (units \(G = c = 1\)) decompose into odd- and even-parity sectors per \((\ell, m)\); each sector reduces to a single master field \(\psi(t, r)\) obeying \[\begin{equation} \frac{\partial^2 \psi}{\partial t^2} - \frac{\partial^2 \psi}{\partial r_*^2} + V(r)\,\psi = 0, \qquad \frac{dr_*}{dr} = \Big(1 - \frac{2M}{r}\Big)^{-1}, \end{equation}\] with \(V = V_{\mathrm{RW}}\) (odd parity, Regge–Wheeler ), \(V = V_{\mathrm{Z}}\) (even parity, Zerilli ), as given in Eqs. \(\eqref{eq:VRW}\)–\(\eqref{eq:VZ}\), while a minimally coupled massless scalar has \(V = V_{\mathrm{s}}= (1 - 2M/r)[\ell(\ell+1)/r^2 + 2M/r^3]\).
Setting \(M = 0\): the tortoise coordinate becomes \(r_* = r\); \(V_{\mathrm{RW}}\to \ell(\ell+1)/r^2\) and \(V_{\mathrm{s}}\to \ell(\ell+1)/r^2\) by inspection; and for Zerilli, \[\begin{equation} V_{\mathrm{Z}}\big|_{M = 0} = \frac{2\lambda^2(\lambda+1) r^3}{r^3\,\lambda^2 r^2} = \frac{2(\lambda+1)}{r^2} = \frac{(\ell-1)(\ell+2) + 2}{r^2} = \frac{\ell(\ell+1)}{r^2}, \end{equation}\] using \(\lambda = (\ell-1)(\ell+2)/2\). The common flat master equation \(-\psi'' + [\ell(\ell+1)/r^2]\psi = \omega^2\psi\) is the Riccati–Bessel equation, solved by \(\psi = r\,j_\ell(\omega r)\), regular at the origin (\(\psi \sim r^{\ell+1}\))—precisely the radial mode functions of the dimensionally reduced massless scalar \(\eqref{eq:Hell}\).
At \(M \neq 0\) the three potentials are pairwise distinct; \(V_{\mathrm{RW}}\) and \(V_{\mathrm{Z}}\) are isospectral via a Darboux transformation but not equal, and neither equals \(V_{\mathrm{s}}\). The flat-space identity of the graviton parities with the scalar is thus a genuine degeneration at \(M = 0\), not an artifact of a shared definition—which is what makes its lattice verification (bitwise equality of couplings built from the three different formulas at \(M = 0\), against \(O(10^{-2})\) differences at \(M = 0.1\)) a meaningful check.
The master fields above arise from the Regge–Wheeler gauge-fixing of \(h_{\mu\nu}\) on Schwarzschild (tensor-harmonic conventions as in ); the variables of arise from a gauge fixing adapted to the \(t = 0\) ball in flat space. At \(M = 0\) both procedures land on the same reduced dynamics—two decoupled scalar towers with \(\ell \geq 2\)—as they must, since the gauge-invariant content per \((\ell, m)\) is two radiative modes and the reduced dynamics is gauge-independent. We use the master-equation route as the anchor because its \(M \to 0\) limit makes the equivalence elementary and machine-checkable.
Each reduced mode is a half-line theory \(H = \tfrac12\int dr\,[\pi^2 + (\partial_r\psi)^2 + V(r)\psi^2]\), discretized on \(N\) sites at unit spacing as \(H = \tfrac12 \pi^{\mathsf T}\pi + \tfrac12 \psi^{\mathsf T} K \psi\). We use two schemes. (i) Srednicki’s : \(K_{11} = \tfrac94 + \ell(\ell+1)\), \(K_{jj} = [(j-\tfrac12)^2 + (j+\tfrac12)^2 + \ell(\ell+1)]/j^2\), \(K_{j,j+1} = -(j+\tfrac12)^2/[j(j+1)]\), the standard benchmark for the area law. (ii) A direct V-scheme for the master fields: \(K_{jj} = 2 + V(j)\), \(K_{j,j+1} = -1\), Dirichlet at both ends, with \(V\) taken from the mode’s own potential formula—\(V_{\mathrm{s}}\), \(V_{\mathrm{RW}}|_{M\to 0}\), or \(V_{\mathrm{Z}}|_{M\to 0}\)—so that scalar and graviton couplings are computed through different expressions and their equality is informative.
Ground-state correlators are \(X = \tfrac12 K^{-1/2}\), \(P = \tfrac12 K^{1/2}\). Restricting to the inner \(n < N\) sites (\(X_A\), \(P_A\)), the symplectic eigenvalues \(\nu_k \geq \tfrac12\) are the positive square roots of the eigenvalues of \(X_A^{1/2} P_A X_A^{1/2}\), and the mode entropy is \(S = \sum_k [(\nu_k + \tfrac12)\ln(\nu_k + \tfrac12) - (\nu_k - \tfrac12)\ln(\nu_k - \tfrac12)]\) . The modular operator of the reduced quasi-free state is determined by the pair \((X_A, P_A)\) together with the symplectic form , which is what makes the per-mode equality of \((X_A, P_A)\) an equality of modular data.
The numerical suite comprises 42 checks in five blocks: (A) scalar baseline—positivity of \(K_\ell\), the Heisenberg bound \(\nu_k \geq \tfrac12\), the area-law coefficient \(S/(n+\tfrac12)^2 \in [0.290, 0.300]\) at \(n = 10, 20\) (\(N = 60\), summed to angular-momentum convergence), monotonicity of \(S_\ell\) in \(\ell\); (B) per-mode identity—bitwise equality of \(K\) from the three potential formulas at \(M = 0\) for \(\ell = 2,\dots,6\) and both parities, trace invariants \(\mathrm{Tr}[(X_A P_A)^k]\), \(k \leq 4\), at \(10^{-15}\), symplectic spectra and entropies at \(10^{-12}\), and pairwise distinctness at \(M = 0.1\); (C) mode bookkeeping—the tower count \(2[(L+1)^2 - 4]\), the ADM count, the missing-modes count; (D) tower sums—the finite-sum identity \(\eqref{eq:tower-sum}\) at \(10^{-12}\) with independently computed sides, positivity of the \(\ell = 0, 1\) deficit, and stability across \(N = 40, 60, 100\); (E) the helicity ladder—Maxwell mode counts \(2[(L+1)^2 - 1]\), the ladder deficits (\(2\) modes scalars \(\to\) Maxwell, \(6\) modes Maxwell \(\to\) graviton), and the entropy-ladder identities \(S_{2\,\mathrm{scalars}} - S_{\mathrm{Maxwell}} = 2 S_0\), \(S_{\mathrm{Maxwell}} - S_{\mathrm{grav}} = 6 S_1\) at \(10^{-12}\), with strict ordering. The symbolic companion verifies the degeneration \(\eqref{eq:degeneration}\) at general symbolic \(\ell\) with identically vanishing residuals (ten identities).
99
J. J. Bisognano and E. H. Wichmann, On the duality condition for quantum fields, J. Math. Phys. 17, 303 (1976). doi:10.1063/1.522898
P. D. Hislop and R. Longo, Modular structure of the local algebras associated with the free massless scalar field theory, Commun. Math. Phys. 84, 71 (1982). doi:10.1007/BF01208372
H. Casini, M. Huerta, and R. C. Myers, Towards a derivation of holographic entanglement entropy, J. High Energy Phys. 05 (2011) 036. doi:10.1007/JHEP05(2011)036
V. Benedetti and H. Casini, Entanglement entropy of linearized gravitons in a sphere, Phys. Rev. D 101, 045004 (2020), arXiv:1908.01800. doi:10.1103/PhysRevD.101.045004
M. Huerta and G. van der Velde, Modular Hamiltonian of the scalar in the semi infinite line: dimensional reduction for spherically symmetric regions, J. High Energy Phys. 06 (2023) 097, arXiv:2301.00294. doi:10.1007/JHEP06(2023)097
H. Casini and M. Huerta, Entanglement entropy of a Maxwell field on the sphere, Phys. Rev. D 93, 105031 (2016), arXiv:1512.06182. doi:10.1103/PhysRevD.93.105031
J. R. David and J. Mukherjee, Hyperbolic cylinders and entanglement entropy: gravitons, higher spins, \(p\)-forms, J. High Energy Phys. 01 (2021) 202, arXiv:2005.08402. doi:10.1007/JHEP01(2021)202
J. R. David and J. Mukherjee, Entanglement entropy of gravitational edge modes, J. High Energy Phys. 08 (2022) 065, arXiv:2201.06043. doi:10.1007/JHEP08(2022)065
A. Blommaert and S. Colin-Ellerin, Gravitons on the edge, J. High Energy Phys. 03 (2025) 116, arXiv:2405.12276. doi:10.1007/JHEP03(2025)116
T. Regge and J. A. Wheeler, Stability of a Schwarzschild singularity, Phys. Rev. 108, 1063 (1957). doi:10.1103/PhysRev.108.1063
F. J. Zerilli, Effective potential for even-parity Regge–Wheeler gravitational perturbation equations, Phys. Rev. Lett. 24, 737 (1970). doi:10.1103/PhysRevLett.24.737
K. Martel and E. Poisson, Gravitational perturbations of the Schwarzschild spacetime: A practical covariant and gauge-invariant formalism, Phys. Rev. D 71, 104003 (2005). doi:10.1103/PhysRevD.71.104003
S. Chandrasekhar and S. Detweiler, The quasi-normal modes of the Schwarzschild black hole, Proc. R. Soc. Lond. A 344, 441 (1975). doi:10.1098/rspa.1975.0112
S. Chandrasekhar, The Mathematical Theory of Black Holes, Oxford University Press, Oxford (1983).
M. Fierz and W. Pauli, On relativistic wave equations for particles of arbitrary spin in an electromagnetic field, Proc. R. Soc. Lond. A 173, 211 (1939). doi:10.1098/rspa.1939.0140
R. Arnowitt, S. Deser, and C. W. Misner, The dynamics of general relativity, in Gravitation: An Introduction to Current Research, ed. L. Witten, Wiley, New York (1962); reprinted in Gen. Relativ. Gravit. 40, 1997 (2008). doi:10.1007/s10714-008-0661-1
H. Araki, A lattice of von Neumann algebras associated with the quantum theory of a free Bose field, J. Math. Phys. 4, 1343 (1963). doi:10.1063/1.1703912
J.-P. Eckmann and K. Osterwalder, An application of Tomita’s theory of modular Hilbert algebras: Duality for free Bose fields, J. Funct. Anal. 13, 1 (1973). doi:10.1016/0022-1236(73)90062-1
K. Fredenhagen, On the modular structure of local algebras of observables, Commun. Math. Phys. 97, 79 (1985). doi:10.1007/BF01206179
R. Haag, Local Quantum Physics: Fields, Particles, Algebras, 2nd ed., Springer, Berlin (1996).
M. Takesaki, Duality for crossed products and the structure of von Neumann algebras of type III, Acta Math. 131, 249 (1973). doi:10.1007/BF02392041
E. Witten, Gravity and the crossed product, J. High Energy Phys. 10 (2022) 008, arXiv:2112.12828. doi:10.1007/JHEP10(2022)008
V. Chandrasekaran, R. Longo, G. Penington, and E. Witten, An algebra of observables for de Sitter space, J. High Energy Phys. 02 (2023) 082, arXiv:2206.10780. doi:10.1007/JHEP02(2023)082
K. Jensen, J. Sorce, and A. J. Speranza, Generalized entropy for general subregions in quantum gravity, J. High Energy Phys. 12 (2023) 020, arXiv:2306.01837. doi:10.1007/JHEP12(2023)020
M. Srednicki, Entropy and area, Phys. Rev. Lett. 71, 666 (1993). doi:10.1103/PhysRevLett.71.666
H. Casini and M. Huerta, Entanglement entropy in free quantum field theory, J. Phys. A 42, 504007 (2009). doi:10.1088/1751-8113/42/50/504007
H. Casini and M. Huerta, Entanglement entropy for the \(n\)-sphere, Phys. Lett. B 694, 167 (2010), arXiv:1007.1813. doi:10.1016/j.physletb.2010.09.054
L. Diósi, A universal master equation for the gravitational violation of quantum mechanics, Phys. Lett. A 120, 377 (1987). doi:10.1016/0375-9601(87)90681-5
R. Penrose, On gravity’s role in quantum state reduction, Gen. Relat. Gravit. 28, 581 (1996). doi:10.1007/BF02105068
M. Sperzel, The Quantum-Geometric Correspondence: Axioms, Gravitational Decoherence, and the Wheeler–DeWitt Rate, companion paper, in preparation (2026).