Home / Quantum Gravity / Locality, lattices and gravity

Exact universal graviton coupling and periodic dispersion

Status
Theorem 2 (no exact cone) is proved for finite-dimensional, real-analytic, periodic Hermitian or unitary symbols. Theorem 1 (locking) is proved under the assumed hypotheses (B) and (S1)–(S4). Proposition 4 concerns an algebraic necessary condition and Proposition 5 a spin-s model. Two steps are sketches: the complexification in Lemma A and the Weinberg–Witten argument of section 6.
Excludes
Lattice Hamiltonian and Floquet gravitons with a finite-dimensional, real-analytic, periodic, Hermitian or unitary Bloch symbol whose band and amplitudes satisfy (B) and (S1)–(S4): a relativistic low-energy band, simple soft poles, exact soft gauge invariance, non-zero energy coupling, and hard graviton scattering with distinct velocities. Not excluded: soft factors universal to finite order only, and models violating any hypothesis.
Revised

On a lattice the graviton dispersion is periodic in the wave vector, so it cannot be Ω(k)=∣k∣\Omega(k)=\lvert k\rvert everywhere. The question is whether the graviton can still couple to energy and momentum as Weinberg’s soft theorem demands, with one coupling for every species and an exactly gauge-invariant soft factor. Under hypotheses on the soft expansion of the amplitude of a hard graviton-graviton scattering process (section 2), which are assumed and not derived from a lattice action, the answer is no. Exact universality forces the graviton band to be an exact cone, Ω(k)=Q(k)\Omega(k)=\sqrt{Q(k)} with QQ a positive quadratic form, on a neighbourhood of k=0k=0 (Theorem 1). No band of a lattice Hamiltonian or Floquet symbol that is a finite Hermitian or unitary matrix, periodic and real-analytic in the wave vector, is an exact cone on any open set (Theorem 2).

The necessary condition used in Theorem 1 can be met only to finite order in the lattice spacing: if Ω2−∣k∣2\Omega^2-\lvert k\rvert^2 begins at degree d≥3d\ge3, the condition fails at relative order (Ea)d−2(Ea)^{d-2} in the energy EE and the lattice spacing aa, and no admissible coupling does better (Proposition 4). It is open whether a lattice amplitude exists whose soft factor is gauge invariant to that order. In a model where the one-particle space is a spin-ss representation, a local analytic energy density has off-forward matrix elements of order one between neighbouring momenta, so the covariance hypothesis of the Weinberg–Witten theorem is not met (Proposition 5).

Weinberg derived the universality of the soft graviton coupling, and with it the equality of gravitational and inertial mass, from the pole structure of amplitudes with a soft graviton in a Lorentz-covariant S-matrix [Weinberg 1964] [Weinberg 1965]. Weinberg and Witten showed that a massless particle of helicity above one has no Lorentz-covariant conserved energy-momentum tensor [Weinberg–Witten]. Porrati extended that argument to theories without a gauge-invariant stress tensor and concluded that massless particles interacting with gravity have spin at most two [Porrati], and Benincasa and Cachazo excluded spin above two from consistency conditions on four-particle amplitudes [Benincasa–Cachazo]. All of these assume exact Poincaré covariance. For supersymmetry, no lattice field theory in infinite volume with finitely many flavours has translation invariance, locality and the Leibniz rule together [Kato–Sakamoto–So]; the obstruction below has a similar shape, with lattice translations and locality on one side and an exact soft Ward identity on the other. The page A conserved lattice momentum and graviton scattering shows that, in classical Hamiltonian lattice theories under the hypotheses listed there, a conserved lattice momentum with an analytic symbol is incompatible with low-energy graviton scattering. That is the charge side of the obstruction; here the dispersion is locked through the soft theorem, and two lemmas of that page are restated in section 3. The lattice spacing is aa throughout, set to 1 where a statement concerns a symbol periodic in 2π2\pi, and the low-energy speed of light is 1.

1. Weinberg’s argument on the light cone

Consider a process with hard particles i=1,…,ni=1,\dots,n of momenta kik_i and energies Ωi\Omega_i, incoming (σi=+1\sigma_i=+1) or outgoing (σi=−1\sigma_i=-1), and add one soft graviton of momentum K=(K0,K)K=(K_0,\mathbf K) and polarisation εμν\varepsilon_{\mu\nu}. The amplitude has poles in KK from the propagators of the hard legs next to the emission. In the continuum the leading term is

An+1≃κ∑iσipiμpiνεμνpi⋅K An,pi=(Ωi,ki).(1)\mathcal A_{n+1}\simeq\kappa\sum_i\sigma_i\frac{p_i^\mu p_i^\nu\varepsilon_{\mu\nu}}{p_i\cdot K}\,\mathcal A_n ,\qquad p_i=(\Omega_i,k_i). \tag{1}

Under a gauge transformation εμν→εμν+Kμξν+ξμKν\varepsilon_{\mu\nu}\to\varepsilon_{\mu\nu}+K_\mu\xi_\nu+\xi_\mu K_\nu with K⋅ξ=0K\cdot\xi=0, the expression (1) changes by 2κ ξ⋅∑iσipi An2\kappa\,\xi\cdot\sum_i\sigma_ip_i\,\mathcal A_n, and this vanishes because energy-momentum is conserved. The coupling κ\kappa is the same for every species because it multiplies the conserved four-vector ∑iσipi\sum_i\sigma_ip_i and nothing else. If two species had couplings κA≠κB\kappa_A\ne\kappa_B, the variation would be proportional to κA∑Api+κB∑Bpi\kappa_A\sum_Ap_i+\kappa_B\sum_Bp_i, which does not vanish in general. This is the universality of the gravitational coupling.

For a photon the same reasoning, with εμ→εμ+Kμ\varepsilon_\mu\to\varepsilon_\mu+K_\mu, gives ∑iσiei=0\sum_i\sigma_ie_i=0 for the charges eie_i. The soft charge is then an internal scalar, conserved on any lattice with an exact U(1)U(1) whatever the dispersion. For the graviton the conserved soft charge is the four-vector ∑iσipi\sum_i\sigma_ip_i, which is tied to the kinematics of the hard legs. On a lattice the question is therefore whether the dispersion allows the charge that gauge invariance requires to be the lattice energy-momentum. Under the hypotheses of section 2, it does so only for an exact cone (Theorem 1).

2. The soft expansion off the cone

Let the graviton band be a function Ω(k)\Omega(k) on a punctured neighbourhood of the origin. In a Hamiltonian lattice theory Ω\Omega is an eigenvalue of the Bloch Hamiltonian H(k)H(k). In a Floquet theory with time step 1, e−iΩ(k)e^{-i\Omega(k)} is an eigenvalue of the unitary Bloch matrix W(k)W(k), and near k=0k=0 the branch with Ω→0\Omega\to0 is used, so that energies are conserved as real numbers with no wrap-around. A hard leg with momentum kik_i that emits a soft graviton has the denominator Ω(ki+K)−Ω(ki)−K0≃−si⋅K\Omega(k_i+\mathbf K)-\Omega(k_i)-K_0\simeq-\mathsf s_i\cdot K, where

si⋅K=K0−vi⋅K,vi=∇Ω(ki).(2)\mathsf s_i\cdot K=K_0-v_i\cdot\mathbf K,\qquad v_i=\nabla\Omega(k_i). \tag{2}

The covector si=(1,−vi)\mathsf s_i=(1,-v_i) is the normal to the on-shell variety {K0=Ω(K)}\{K_0=\Omega(\mathbf K)\} at kik_i. On the cone it is proportional to pip_i, and (2) reduces to −pi⋅K/Ωi-p_i\cdot K/\Omega_i.

The lattice theories considered here have no covariant S-matrix from which the soft expansion could be derived, so we take its form as hypothesis. The first hypothesis concerns the band.

(B) After a linear change of the spatial axes, Ω(k)=∣k∣ ω~(∣k∣,k^)\Omega(k)=\lvert k\rvert\,\tilde\omega(\lvert k\rvert,\hat k) for 0<∣k∣<ρ00<\lvert k\rvert<\rho_0, where ω~\tilde\omega is real-analytic in a neighbourhood of [0,ρ0)×S2[0,\rho_0)\times S^2 and ω~(0,k^)=1\tilde\omega(0,\hat k)=1. Then Ω=∑m≥1Ωm\Omega=\sum_{m\ge1}\Omega_m near k=0k=0, with Ωm\Omega_m real-analytic on R3∖{0}\mathbb R^3\setminus\{0\} and positively homogeneous of degree mm, and Ω1=∣k∣\Omega_1=\lvert k\rvert. This says that the dispersion is relativistic at low energy and has a Taylor-type expansion in the scale. It holds, for example, when Ω2\Omega^2 is real-analytic at 00 with Ω2=∣k∣2+O(k3)\Omega^2=\lvert k\rvert^2+O(k^3).

Write η=diag(−1,1,1,1)\eta=\mathrm{diag}(-1,1,1,1) for the metric whose null cone K0=∣K∣K_0=\lvert\mathbf K\rvert is the leading cone of the band, and let si♯=−(1,vi)\mathsf s_i^\sharp=-(1,v_i) be the vector with η(si♯,x)=si⋅x\eta(\mathsf s_i^\sharp,x)=\mathsf s_i\cdot x for all xx; then η(si♯,si♯)=∣vi∣2−1\eta(\mathsf s_i^\sharp,\mathsf s_i^\sharp)=\lvert v_i\rvert^2-1. The remaining hypotheses concern the soft factor.

(S1) For a hard 2→22\to2 graviton process at small momenta with A4≠0\mathcal A_4\ne0, the amplitude with an extra soft graviton of polarisation εμν\varepsilon_{\mu\nu} and momentum KK on the leading cone is, as K→0K\to0,

A5=[∑iσiεμνJiμνsi⋅K+O(K0)]A4,\mathcal A_5=\Big[\sum_i\sigma_i\frac{\varepsilon_{\mu\nu}J_i^{\mu\nu}}{\mathsf s_i\cdot K}+O(K^0)\Big]\mathcal A_4 ,

where Jiμν=Jμν(ki)J_i^{\mu\nu}=J^{\mu\nu}(k_i) is a symmetric tensor that depends on the leg ii only, acts as a multiple of the identity on the helicity index, and has the form J^(∣k∣,k^)\hat J(\lvert k\rvert,\hat k) with J^\hat J real-analytic near [0,ρ0)×S2[0,\rho_0)\times S^2. In general relativity Jiμν=−κΩi si♯μsi♯νJ_i^{\mu\nu}=-\kappa\Omega_i\,\mathsf s_i^{\sharp\mu}\mathsf s_i^{\sharp\nu}, which is (1) rewritten with pi⋅K=−Ωi si⋅Kp_i\cdot K=-\Omega_i\,\mathsf s_i\cdot K.

(S2) The pole term of (S1) is invariant under εμν→εμν+Kμξν+ξμKν\varepsilon_{\mu\nu}\to\varepsilon_{\mu\nu}+K_\mu\xi_\nu+\xi_\mu K_\nu for every KK on the leading cone and every ξ\xi with η(K,ξ)=0\eta(K,\xi)=0 (indices lowered with η\eta). The variation of the pole term is of order K0K^0 and that of the next term is of order KK, so the invariance of the amplitude implies that of the pole term separately.

The next hypothesis needs the set of low-energy collisions. Let Λ={q=(q1,…,q4)∈(R3)4: q1+q2=q3+q4}\Lambda=\{q=(q_1,\dots,q_4)\in(\mathbb R^3)^4:\ q_1+q_2=q_3+q_4\} be the space of momentum-conserving quartets, with σ=(+,+,−,−)\sigma=(+,+,-,-), and for ρ>0\rho>0 let

Mρ={q∈Λ: ∑iσiΩ(ρqi)=0}\mathcal M_\rho=\Big\{q\in\Lambda:\ \sum_i\sigma_i\Omega(\rho q_i)=0\Big\}

be the rescaled energy shell. Its limit as ρ→0\rho\to0 is the massless resonant set M0={q∈Λ: ∣q1∣+∣q2∣=∣q3∣+∣q4∣}\mathcal M_0=\{q\in\Lambda:\ \lvert q_1\rvert+\lvert q_2\rvert=\lvert q_3\rvert+\lvert q_4\rvert\}, and M0∗\mathcal M_0^\ast is its part on which ∣q1∣+∣q2∣>∣q1+q2∣\lvert q_1\rvert+\lvert q_2\rvert>\lvert q_1+q_2\rvert.

(S3) There are an open set U⊂ΛU\subset\Lambda meeting M0∗\mathcal M_0^\ast and a number ρ1>0\rho_1>0 such that for 0<ρ<ρ10<\rho<\rho_1 the quartets ρq\rho q with q∈Mρ∩Uq\in\mathcal M_\rho\cap U have A4(ρq)≠0\mathcal A_4(\rho q)\ne0 and pairwise distinct velocities viv_i.

(S4) The tensor JJ is not identically a multiple of ημν\eta^{\mu\nu}. (By Lemma A this means that the coefficient β\beta of s♯s♯\mathsf s^\sharp\mathsf s^\sharp in JJ is not identically zero: the graviton couples to energy and momentum.)

In general relativity Ω=∣k∣\Omega=\lvert k\rvert, s♯=−p/Ω\mathsf s^\sharp=-p/\Omega is null, and β=−κΩ\beta=-\kappa\Omega; the hypotheses hold in that case. We do not assert that they hold for a lattice graviton; Theorem 1 states what they imply.

3. Exact universality locks the dispersion

Let J^i\hat J_i be the operator (J^iξ)μ=Jiμνηνρξρ(\hat J_i\xi)^\mu=J_i^{\mu\nu}\eta_{\nu\rho}\xi^\rho, which is self-adjoint for η\eta because JiJ_i is symmetric. Lemma A says what (S1) and (S2) leave of JiJ_i, and what condition on the hard momenta results.

Lemma A (the soft vertex). Assume (S1) and (S2), and that for the process considered the velocities viv_i are pairwise distinct with ∣vi∣≠1\lvert v_i\rvert\ne1. Then Jiμν=αiημν+βi si♯μsi♯νJ_i^{\mu\nu}=\alpha_i\eta^{\mu\nu}+\beta_i\,\mathsf s_i^{\sharp\mu}\mathsf s_i^{\sharp\nu} for scalars αi,βi\alpha_i,\beta_i that depend on kik_i only, and

∑iσiβi si♯=0,that is∑iσi β(ki) (1,∇Ω(ki))=0.(3)\sum_i\sigma_i\beta_i\,\mathsf s_i^\sharp=0,\qquad\text{that is}\qquad\sum_i\sigma_i\,\beta(k_i)\,\big(1,\nabla\Omega(k_i)\big)=0. \tag{3}

Moreover, assume (B) and suppose that the decomposition holds, with scalars α(k),β(k)\alpha(k),\beta(k), for all kk in a non-empty open subset of a small punctured ball 0<∣k∣<ρ0′0<\lvert k\rvert<\rho_0'. Then β\beta and α\alpha extend to functions of (∣k∣,k^)(\lvert k\rvert,\hat k) that are real-analytic near ∣k∣=0\lvert k\rvert=0, and the decomposition of JJ holds at every kk of the punctured ball, including where ∣v∣=1\lvert v\rvert=1.

Proof. The leg ii is the only term of the gauge variation with a pole at si⋅K=0\mathsf s_i\cdot K=0, so its numerator must vanish there. This fixes JiJ_i up to two scalars, and what remains of the variation is (3).

With the operator J^i\hat J_i the variation in (S2) is

∑iσiη(K,J^iξ)si⋅K=0for real null K, η(K,ξ)=0.(4)\sum_i\sigma_i\frac{\eta(K,\hat J_i\xi)}{\mathsf s_i\cdot K}=0\qquad\text{for real null }K,\ \eta(K,\xi)=0. \tag{4}

Multiplying by ∏j(sj⋅K)\prod_j(\mathsf s_j\cdot K) turns (4) into a polynomial identity in (K,ξ)(K,\xi) that holds on the real points of the complex variety V={(K,ξ)∈C4×C4: η(K,K)=0, η(K,ξ)=0}\mathcal V=\{(K,\xi)\in\mathbb C^4\times\mathbb C^4:\ \eta(K,K)=0,\ \eta(K,\xi)=0\}. We need it on the complex points, because the real null vectors in the hyperplane si⋅K=0\mathsf s_i\cdot K=0 may be absent when ∣vi∣<1\lvert v_i\rvert<1. The passage from real to complex points is the standard fact that a polynomial vanishing on a real-analytic submanifold of an irreducible complex variety, of real dimension equal to the complex dimension of the variety, vanishes on the whole variety. Here V\mathcal V is irreducible (a rank-3 bundle over the irreducible complex null cone minus the origin; the part with K=0K=0 has dimension 4, less than 6, and is not a component), and its real points K≠0K\ne0 null, ξ\xi real with η(K,ξ)=0\eta(K,\xi)=0 form a manifold of real dimension 6, the complex dimension of V\mathcal V. This step is a sketch of the standard argument, without proof.

Fix ii. Because ∣vi∣≠1\lvert v_i\rvert\ne1, si♯\mathsf s_i^\sharp is not null, and the hyperplane si⋅K=0\mathsf s_i\cdot K=0 is a three-dimensional space on which η\eta is non-degenerate. Its complex null vectors form a quadric cone CiC_i that spans it, and Ci∖{0}C_i\setminus\{0\} is connected. Because the velocities are distinct, the hyperplanes sj⋅K=0\mathsf s_j\cdot K=0 for j≠ij\ne i do not contain CiC_i, so at a generic point of CiC_i only the term ii of the multiplied identity survives. It gives η(K,J^iξ)=0\eta(K,\hat J_i\xi)=0 for K∈CiK\in C_i and all ξ\xi orthogonal to KK. By self-adjointness η(J^iK,ξ)=0\eta(\hat J_iK,\xi)=0 for all ξ⊥K\xi\perp K, so J^iK\hat J_iK is a multiple of KK: J^iK=λ(K)K\hat J_iK=\lambda(K)K on CiC_i. The eigenvalue λ(K)\lambda(K) is continuous on the connected set Ci∖{0}C_i\setminus\{0\} and takes values in the finite spectrum of J^i\hat J_i, so it is a constant λi\lambda_i, and J^i=λi\hat J_i=\lambda_i on the span of CiC_i, the hyperplane (si♯)⊥(\mathsf s_i^\sharp)^\perp. A self-adjoint operator that is λi\lambda_i on (si♯)⊥(\mathsf s_i^\sharp)^\perp maps the η\eta-orthogonal complement of that hyperplane, the line through the non-null vector si♯\mathsf s_i^\sharp, into itself: J^isi♯=μisi♯\hat J_i\mathsf s_i^\sharp=\mu_i\mathsf s_i^\sharp. Hence J^i=λi+(μi−λi) si♯ η(si♯,⋅)/η(si♯,si♯)\hat J_i=\lambda_i+(\mu_i-\lambda_i)\,\mathsf s_i^\sharp\,\eta(\mathsf s_i^\sharp,\cdot)/\eta(\mathsf s_i^\sharp,\mathsf s_i^\sharp), which is the stated form with αi=λi\alpha_i=\lambda_i and βi=(μi−λi)/η(si♯,si♯)\beta_i=(\mu_i-\lambda_i)/\eta(\mathsf s_i^\sharp,\mathsf s_i^\sharp).

Insert this form into (4). The αi\alpha_i term gives αiη(K,ξ)=0\alpha_i\eta(K,\xi)=0. The βi\beta_i term gives βi η(K,si♯) η(si♯,ξ)/(si⋅K)=βi η(si♯,ξ)\beta_i\,\eta(K,\mathsf s_i^\sharp)\,\eta(\mathsf s_i^\sharp,\xi)/(\mathsf s_i\cdot K)=\beta_i\,\eta(\mathsf s_i^\sharp,\xi), because η(K,si♯)=si⋅K\eta(K,\mathsf s_i^\sharp)=\mathsf s_i\cdot K. So η(y,ξ)=0\eta(y,\xi)=0 for y=∑iσiβisi♯y=\sum_i\sigma_i\beta_i\mathsf s_i^\sharp, every real null KK and every ξ⊥K\xi\perp K. Then yy is parallel to every real null vector, hence y=0y=0. Since si♯=−(1,vi)\mathsf s_i^\sharp=-(1,v_i) this is (3).

For the last statement, take a covector n=(1,k^)\mathsf n=(1,\hat k), which is null for the inverse metric, and contract the decomposition with nμnν\mathsf n_\mu\mathsf n_\nu. Where the decomposition holds, β=Jμνnμnν/(s♯ ⁣⋅n)2\beta=J^{\mu\nu}\mathsf n_\mu\mathsf n_\nu/(\mathsf s^\sharp\!\cdot\mathsf n)^2 with s♯ ⁣⋅n=−(1+v⋅k^)\mathsf s^\sharp\!\cdot\mathsf n=-(1+v\cdot\hat k). For small ∣k∣\lvert k\rvert, v→k^v\to\hat k and the denominator tends to −2-2, so the right-hand side is real-analytic in (∣k∣,k^)(\lvert k\rvert,\hat k) near ∣k∣=0\lvert k\rvert=0, and we take it as the definition of β\beta on the whole punctured ball. Likewise α=(Jμνημν−β η(s♯,s♯))/4\alpha=(J^{\mu\nu}\eta_{\mu\nu}-\beta\,\eta(\mathsf s^\sharp,\mathsf s^\sharp))/4 is analytic. The difference J−αη−βs♯s♯J-\alpha\eta-\beta\mathsf s^\sharp\mathsf s^\sharp is then real-analytic on the punctured ball, which is connected, and vanishes on the open set where the decomposition holds, so it vanishes identically. This is the decomposition at every kk, including where ∣v∣=1\lvert v\rvert=1. In particular β=∑m≥0βm\beta=\sum_{m\ge0}\beta_m with βm\beta_m real-analytic on R3∖{0}\mathbb R^3\setminus\{0\} and homogeneous of degree mm.

For the exact cone Ω=∣k∣\Omega=\lvert k\rvert, (3) reads ∑iσiβ(ki)(1,k^i)=0\sum_i\sigma_i\beta(k_i)(1,\hat k_i)=0, which is energy-momentum conservation when β∝∣k∣\beta\propto\lvert k\rvert; this is Weinberg’s case. On a general band, (3) says that the vector β (1,∇Ω)\beta\,(1,\nabla\Omega) is conserved in every hard process, a condition on β\beta and Ω\Omega together.

The proof of Theorem 1 expands (3) in the scale of the hard momenta, and at each order it needs to know which functions are conserved in every massless collision. The answer is in two lemmas of A conserved lattice momentum and graviton scattering, restated here with Δg(q)=g(q1)+g(q2)−g(q3)−g(q4)\Delta g(q)=g(q_1)+g(q_2)-g(q_3)-g(q_4).

Lemma B (massless collisions). Parts (i) and (ii) are Lemmas 3 and 4 of the page cited above, where the sets Λ\Lambda, Mρ\mathcal M_\rho, M0\mathcal M_0, M0∗\mathcal M_0^\ast appear with the same meaning (the scale is called λ\lambda there); part (iii) follows from (i) and (B) and is proved below.

(i) M0∗\mathcal M_0^\ast is a connected real-analytic submanifold of Λ\Lambda of dimension 8. On it all qi≠0q_i\ne0, and the gradient along Λ\Lambda of F0(q)=∣q1∣+∣q2∣−∣q3∣−∣q4∣F_0(q)=\lvert q_1\rvert+\lvert q_2\rvert-\lvert q_3\rvert-\lvert q_4\rvert does not vanish.

(ii) Let g:R3∖{0}→Rℓg:\mathbb R^3\setminus\{0\}\to\mathbb R^{\ell} be real-analytic and positively homogeneous of a real degree mm. If Δg\Delta g vanishes on a non-empty relatively open subset of M0∗\mathcal M_0^\ast, then g≡0g\equiv0 unless m∈{0,1}m\in\{0,1\}. In degree 1, the linear functions and ∣q∣\lvert q\rvert satisfy Δg=0\Delta g=0 on all of M0\mathcal M_0; in degree 0, so do the constants.

(iii) Assume (B). Every q∗∈M0∗q^\ast\in\mathcal M_0^\ast is a limit of points qρ∈Mρq_\rho\in\mathcal M_\rho with qρ−q∗=O(ρ)q_\rho-q^\ast=O(\rho).

Part (ii) rests on the quartet (3z,z;2z,2z)(3z,z;2z,2z), which lies in the closure of M0∗\mathcal M_0^\ast and gives, for homogeneous gg of degree mm, Δg=(3m+1−2m+1) g(z)\Delta g=(3^m+1-2^{m+1})\,g(z). After division by 2m2^m the factor is (3/2)m+2−m−2(3/2)^m+2^{-m}-2, which is strictly convex in mm and vanishes only at m=0,1m=0,1. For part (iii), set Fρ(q)=ρ−1∑iσiΩ(ρqi)F_\rho(q)=\rho^{-1}\sum_i\sigma_i\Omega(\rho q_i), so that Mρ={Fρ=0}\mathcal M_\rho=\{F_\rho=0\}. By (B), Fρ→F0F_\rho\to F_0 in C1C^1 with an error O(ρ)O(\rho) near points where all qi≠0q_i\ne0, and by (i) the gradient of F0F_0 does not vanish there, so the implicit function theorem gives (iii).

Theorem 1 (locking of the dispersion). Assume (B) and (S1)–(S4). Then Ω(k)=∣k∣\Omega(k)=\lvert k\rvert for 0<∣k∣<ρ00<\lvert k\rvert<\rho_0, in the coordinates of (B): the band is an exact cone.

Proof. We expand (3) in the scale ρ\rho of the hard momenta. At each order the conserved vector is a collision invariant of massless kinematics, and by Lemma B(ii) the only homogeneous invariants are the energy-momentum itself and constants. This forces β=γΩ\beta=\gamma\Omega for a constant γ\gamma, and then removes every correction to the cone.

Step 0. The function ∣∇Ω∣2−1\lvert\nabla\Omega\rvert^2-1 is real-analytic on the punctured ball, and we first suppose that it vanishes identically. Its part of degree m−1m-1, m≥2m\ge2, is 2k^⋅∇Ωm+∑∇Ωm1⋅∇Ωm22\hat k\cdot\nabla\Omega_m+\sum\nabla\Omega_{m_1}\cdot\nabla\Omega_{m_2} summed over m1+m2=m+1m_1+m_2=m+1 with m1,m2≥2m_1,m_2\ge2. If Ω2=⋯=Ωm−1=0\Omega_2=\dots=\Omega_{m-1}=0 the sum vanishes, and Euler’s relation k⋅∇Ωm=mΩmk\cdot\nabla\Omega_m=m\Omega_m gives 2m Ωm/∣k∣=02m\,\Omega_m/\lvert k\rvert=0. By induction Ωm=0\Omega_m=0 for all m≥2m\ge2, so Ω=∣k∣\Omega=\lvert k\rvert near 00 and, by analytic continuation, on 0<∣k∣<ρ00<\lvert k\rvert<\rho_0.

Otherwise the function is not identically zero on the connected punctured ball, and its zero set ZZ is nowhere dense. Shrink UU to a neighbourhood of a point of M0∗∩U\mathcal M_0^\ast\cap U. For small ρ\rho, Mρ\mathcal M_\rho is then a smooth hypersurface in UU, and the projection to the momentum qiq_i of any one leg is a submersion. For leg 1, take q1q_1 fixed and q4=q1+q2−q3q_4=q_1+q_2-q_3: then ∂Fρ/∂q3→−q^3+q^4\partial F_\rho/\partial q_3\to-\hat q_3+\hat q_4, which is non-zero because ∣q3∣+∣q4∣=∣q1∣+∣q2∣>∣q1+q2∣=∣q3+q4∣\lvert q_3\rvert+\lvert q_4\rvert=\lvert q_1\rvert+\lvert q_2\rvert>\lvert q_1+q_2\rvert=\lvert q_3+q_4\rvert on M0∗\mathcal M_0^\ast, so q3q_3 and q4q_4 are not positively parallel. The other legs are treated in the same way. The preimage of the nowhere dense set ZZ under the open map q↦ρqiq\mapsto\rho q_i is nowhere dense, so the quartets with ∣vi∣≠1\lvert v_i\rvert\ne1 for all four legs are open and dense in Mρ∩U\mathcal M_\rho\cap U. By (S3) their velocities are pairwise distinct. Lemma A applies there, and (3) holds on a dense subset of Mρ∩U\mathcal M_\rho\cap U, hence by continuity on all of it, for each small ρ\rho. The leg momenta ρqi\rho q_i of the quartets with ∣vi∣≠1\lvert v_i\rvert\ne1 fill a non-empty open set of the punctured ball, so by the last statement of Lemma A the decomposition of JJ holds at every kk with 0<∣k∣<ρ0′0<\lvert k\rvert<\rho_0'.

Step 1. Put Ψ(k)=β(k) (1,∇Ω(k))=∑mΨm(k)\Psi(k)=\beta(k)\,(1,\nabla\Omega(k))=\sum_m\Psi_m(k) with Ψm\Psi_m homogeneous of degree mm and real-analytic off the origin. For q∈Mρ∩Uq\in\mathcal M_\rho\cap U, (3) reads ∑iσiΨ(ρqi)=0\sum_i\sigma_i\Psi(\rho q_i)=0. Momentum is conserved exactly, ∑iσiqi=0\sum_i\sigma_iq_i=0 on Λ\Lambda, and by Lemma B(iii) every q∗∈M0∗∩Uq^\ast\in\mathcal M_0^\ast\cap U is a limit of points qρ∈Mρ∩Uq_\rho\in\mathcal M_\rho\cap U with qρ−q∗=O(ρ)q_\rho-q^\ast=O(\rho).

Step 2. The lowest order fixes the coupling. Let m∗m_\ast be the lowest degree of β\beta, so that βm∗≢0\beta_{m_\ast}\not\equiv0 by (S4) and Lemma A. Since ∇Ω=q^+O(ρ)\nabla\Omega=\hat q+O(\rho), we have Ψ(ρq)=ρm∗βm∗(q)(1,q^)+O(ρm∗+1)\Psi(\rho q)=\rho^{m_\ast}\beta_{m_\ast}(q)(1,\hat q)+O(\rho^{m_\ast+1}), and the limit ρ→0\rho\to0 gives ∑iσiβm∗(qi)(1,q^i)=0\sum_i\sigma_i\beta_{m_\ast}(q_i)(1,\hat q_i)=0 on M0∗∩U\mathcal M_0^\ast\cap U. This extends to all of M0∗\mathcal M_0^\ast by analyticity and Lemma B(i). The energy component is a homogeneous collision invariant of degree m∗m_\ast, so m∗∈{0,1}m_\ast\in\{0,1\} by Lemma B(ii).

Write βm∗(q)=b(q)∣q∣\beta_{m_\ast}(q)=b(q)\lvert q\rvert and pi=(∣qi∣,qi)p_i=(\lvert q_i\rvert,q_i) for the null four-momenta. Then ∑iσib(qi) pi=0\sum_i\sigma_ib(q_i)\,p_i=0, and also ∑iσipi=0\sum_i\sigma_ip_i=0. For a quartet in M0∗\mathcal M_0^\ast with q1,q2,q3q_1,q_2,q_3 linearly independent, p1,p2,p3p_1,p_2,p_3 are independent and p4=p1+p2−p3p_4=p_1+p_2-p_3, so the four vectors span a three-dimensional space and satisfy exactly one linear relation up to scale. Comparing the two relations gives b(q1)=b(q2)=b(q3)=b(q4)b(q_1)=b(q_2)=b(q_3)=b(q_4).

Given two linearly independent non-zero vectors t1,t2t_1,t_2 there is such a quartet with q1=t1q_1=t_1 and q3=t2q_3=t_2. Take q4=q2+t1−t2q_4=q_2+t_1-t_2, so that momentum is conserved, and require q2=xq_2=x to lie on the surface ∣x∣−∣x+t1−t2∣=∣t2∣−∣t1∣\lvert x\rvert-\lvert x+t_1-t_2\rvert=\lvert t_2\rvert-\lvert t_1\rvert, which expresses energy conservation. The right-hand side lies strictly between ±∣t1−t2∣\pm\lvert t_1-t_2\rvert because t1,t2t_1,t_2 are independent, so the surface is a sheet of a hyperboloid with foci 00 and t2−t1t_2-t_1 (a plane if ∣t1∣=∣t2∣\lvert t_1\rvert=\lvert t_2\rvert). It is non-empty, because on the focal line the left-hand side runs continuously from −∣t1−t2∣-\lvert t_1-t_2\rvert to ∣t1−t2∣\lvert t_1-t_2\rvert; it is smooth because the gradient of ∣x∣−∣x+t1−t2∣\lvert x\rvert-\lvert x+t_1-t_2\rvert vanishes only on the focal line, where the left-hand side equals ±∣t1−t2∣\pm\lvert t_1-t_2\rvert, and it is invariant under rotations about the focal line, which lies in the plane spanned by t1,t2t_1,t_2. A two-dimensional surface is not contained in a line, so it has a point off the focal line, and the circle of that point under the rotations meets the plane in at most two points. Hence there is a point xx off the plane spanned by t1,t2t_1,t_2, and it completes the quartet. So b(t1)=b(t2)b(t_1)=b(t_2) for independent t1,t2t_1,t_2. For any two vectors there is a third one independent of both, so bb is constant, b≡γb\equiv\gamma, and γ≠0\gamma\ne0 because βm∗≢0\beta_{m_\ast}\not\equiv0. Then βm∗=γ∣q∣\beta_{m_\ast}=\gamma\lvert q\rvert has degree one, so m∗=1m_\ast=1 and β1=γΩ1\beta_1=\gamma\Omega_1.

Step 3. The energy component gives β=γΩ\beta=\gamma\Omega. On Mρ\mathcal M_\rho the energy is conserved exactly, ∑iσiΩ(ρqi)=0\sum_i\sigma_i\Omega(\rho q_i)=0, so the energy component of Step 1 gives ∑iσi(β−γΩ)(ρqi)=0\sum_i\sigma_i(\beta-\gamma\Omega)(\rho q_i)=0. The function g=β−γΩ=∑m≥2gmg=\beta-\gamma\Omega=\sum_{m\ge2}g_m has no term of degree one. If g≢0g\not\equiv0, let gmg_m be its lowest non-zero term. Then ρm∑iσigm(qρ,i)+O(ρm+1)=0\rho^m\sum_i\sigma_ig_m(q_{\rho,i})+O(\rho^{m+1})=0, so ∑iσigm(qi∗)=0\sum_i\sigma_ig_m(q^\ast_i)=0 on M0∗∩U\mathcal M_0^\ast\cap U, and Lemma B(ii) gives gm≡0g_m\equiv0, a contradiction. Hence β=γΩ\beta=\gamma\Omega as series.

Step 4. The momentum components give the cone. Now Ψj=γΩ ∂jΩ=γ2∂j(Ω2)\Psi^j=\gamma\Omega\,\partial_j\Omega=\tfrac\gamma2\partial_j(\Omega^2). Write Ω2=∣k∣2+R\Omega^2=\lvert k\rvert^2+R with R=∑d≥3RdR=\sum_{d\ge3}R_d, RdR_d homogeneous of degree dd and real-analytic off the origin. Exact momentum conservation gives ∑iσi∂j∣ki∣2=2∑iσiki,j=0\sum_i\sigma_i\partial_j\lvert k_i\rvert^2=2\sum_i\sigma_ik_{i,j}=0, so (3) becomes ∑iσi∂jR(ρqi)=0\sum_i\sigma_i\partial_jR(\rho q_i)=0 on Mρ∩U\mathcal M_\rho\cap U. Suppose R≢0R\not\equiv0 and let RdR_d be its lowest non-zero term. Then ρd−1∑iσi∂jRd(qρ,i)+O(ρd)=0\rho^{d-1}\sum_i\sigma_i\partial_jR_d(q_{\rho,i})+O(\rho^d)=0, so ∂jRd\partial_jR_d is a homogeneous collision invariant of degree d−1≥2d-1\ge2 on M0∗∩U\mathcal M_0^\ast\cap U and vanishes by Lemma B(ii). A homogeneous function of positive degree with zero gradient vanishes, so Rd≡0R_d\equiv0, a contradiction. So R≡0R\equiv0 and Ω2=∣k∣2\Omega^2=\lvert k\rvert^2 near 00, and by analytic continuation on the punctured ball.

The energy row fixes the coupling to be proportional to the energy. The momentum rows say that Ω∇Ω=12∇Ω2\Omega\nabla\Omega=\tfrac12\nabla\Omega^2 is conserved in every hard process, which for a band that is not the cone cannot hold on all the collisions that the lattice allows at low energy.

4. A periodic symbol has no exact cone

Theorem 2 (no exact cone). Let A(k)A(k) be an N×NN\times N matrix function on R3\mathbb R^3, real-analytic and 2π2\pi-periodic in each coordinate. Let QQ be a positive-definite quadratic form and M≥0M\ge0.

(a) If AA is Hermitian (a lattice Hamiltonian), no continuous eigenvalue branch Ω\Omega of A(k)A(k) satisfies Ω(k)2=Q(k)+M2\Omega(k)^2=Q(k)+M^2 on a non-empty open set.

(b) If AA is unitary (a Floquet or cellular-automaton symbol), no branch with e−iΩ(k)∈spec A(k)e^{-i\Omega(k)}\in\mathrm{spec}\,A(k) satisfies Ω(k)2=Q(k)+M2\Omega(k)^2=Q(k)+M^2 on a non-empty open set.

In both cases the open set is taken not to contain k=0k=0 when M=0M=0.

Proof. A determinant identity on a ball spreads by analyticity to all of space, and periodicity then places an unbounded sequence in a finite spectrum. On a small ball in the open set, Ω≠0\Omega\ne0 and is continuous, so Ω=χQ+M2\Omega=\chi\sqrt{Q+M^2} with a fixed sign χ=±1\chi=\pm1. Put ϕ0(k)=Q(k)+M2\phi_0(k)=\sqrt{Q(k)+M^2}, real-analytic on R3\mathbb R^3 (on R3∖{0}\mathbb R^3\setminus\{0\} if M=0M=0).

(a) The function Ξ(k)=det⁡(χϕ0(k)−A(k))\Xi(k)=\det(\chi\phi_0(k)-A(k)) is real-analytic on the connected set R3∖{0}\mathbb R^3\setminus\{0\} and vanishes on the ball, so Ξ≡0\Xi\equiv0: χϕ0(k)∈spec A(k)\chi\phi_0(k)\in\mathrm{spec}\,A(k) for all k≠0k\ne0. By periodicity χϕ0(k+2πne1)∈spec A(k)\chi\phi_0(k+2\pi ne_1)\in\mathrm{spec}\,A(k) for every integer nn. But ϕ0(k+2πne1)→∞\phi_0(k+2\pi ne_1)\to\infty and the spectrum of a fixed N×NN\times N matrix is bounded. Contradiction.

(b) The same argument with Ξ(k)=det⁡(e−iχϕ0(k)−A(k))\Xi(k)=\det(e^{-i\chi\phi_0(k)}-A(k)) gives e−iχϕn∈spec A(k)e^{-i\chi\phi_n}\in\mathrm{spec}\,A(k) for all nn, where ϕn=ϕ0(k+2πne1)\phi_n=\phi_0(k+2\pi ne_1). The sign χ\chi is the same for all nn because it was fixed in Ξ\Xi. The spectrum has at most NN points on the unit circle, so χϕn\chi\phi_n lies, modulo 2π2\pi, in a finite set Θ\Theta of phases, and the increments δn=ϕn+1−ϕn\delta_n=\phi_{n+1}-\phi_n satisfy χδn∈Γ=(Θ−Θ)+2πZ\chi\delta_n\in\Gamma=(\Theta-\Theta)+2\pi\mathbb Z, a set that meets every bounded interval in finitely many points. Now choose kk not parallel to e1e_1 and set ϕ(t)=Q11t2+2Q1kt+Qkk+M2\phi(t)=\sqrt{Q_{11}t^2+2Q_{1k}t+Q_{kk}+M^2}, with Q11=Q(e1)Q_{11}=Q(e_1), Q1k=Q(e1,k)Q_{1k}=Q(e_1,k) the bilinear form and Qkk=Q(k)Q_{kk}=Q(k), so that ϕn=ϕ(2πn)\phi_n=\phi(2\pi n). Then

ϕ′′(t)=Q11(Qkk+M2)−Q1k2(Q11t2+2Q1kt+Qkk+M2)3/2>0,\phi''(t)=\frac{Q_{11}(Q_{kk}+M^2)-Q_{1k}^2}{\big(Q_{11}t^2+2Q_{1k}t+Q_{kk}+M^2\big)^{3/2}}>0,

by the strict Cauchy–Schwarz inequality Q1k2<Q11QkkQ_{1k}^2<Q_{11}Q_{kk} for non-parallel vectors. So ϕ\phi is strictly convex and the increments δn=ϕ(2π(n+1))−ϕ(2πn)\delta_n=\phi(2\pi(n+1))-\phi(2\pi n) are strictly increasing, and they are bounded above by 2πQ112\pi\sqrt{Q_{11}}. A strictly increasing bounded sequence takes infinitely many values in a bounded interval, but χΓ\chi\Gamma meets that interval in finitely many points. Contradiction.

Remark (how much analyticity is used). The proof uses real-analyticity of AA only on the complement of a discrete set that is invariant under the translations by 2πZ32\pi\mathbb Z^3, such as 2πZ32\pi\mathbb Z^3 itself, because that complement is connected. Smoothness off such a set is not enough. Let gg be smooth with g=1g=1 for ∣k∣≤1\lvert k\rvert\le1 and g=0g=0 for ∣k∣≥2\lvert k\rvert\ge2, and let A(k)A(k) be the 2π2\pi-periodic extension of ∣k∣ g(∣k∣)\lvert k\rvert\,g(\lvert k\rvert) from the cube [−π,π]3[-\pi,\pi]^3; this is a scalar (N=1N=1) symbol, smooth away from 2πZ32\pi\mathbb Z^3, equal to the exact cone ∣k∣\lvert k\rvert near the origin. The conical point at k=0k=0 makes the hopping amplitudes, the Fourier coefficients of AA, decay as the power law −1/(π2∣x∣4)-1/(\pi^2\lvert x\rvert^4) in the lattice distance ∣x∣\lvert x\rvert; the transition region of gg is where analyticity fails. This is a toy band, not a graviton. For C∞C^\infty symbols with non-vanishing hard scattering, the question is whether crystal-momentum periodicity alone excludes an exact cone.

Remark (periodicity at complex momentum; conditional). In a formulation with a Euclidean time direction the on-shell variety is the complex cone ∑μpμ2=0\sum_\mu p_\mu^2=0, and the kinetic symbol is not a finite Hermitian or unitary matrix, so Theorem 2 does not apply. Let P(p)P(p), the symbol of a conserved lattice momentum as in A conserved lattice momentum and graviton scattering, be 2π2\pi-periodic in p1p_1 and real-analytic in the closed strip ∣Im pμ∣≤π/3\lvert\mathrm{Im}\,p_\mu\rvert\le\pi/\sqrt3, and assume P(p)=p1P(p)=p_1 on the complex cone inside the strip. This assumption is the Euclidean analogue of Theorem 1 and is not derived. If pp and p−2πe1p-2\pi e_1 are both on the cone, subtracting the two equations of the cone gives 4πp1−4π2=04\pi p_1-4\pi^2=0, so p1=πp_1=\pi and p22+p32+p42=−π2p_2^2+p_3^2+p_4^2=-\pi^2. Writing pj=xj+iyjp_j=x_j+iy_j gives ∑jyj2=π2+∑jxj2≥π2\sum_jy_j^2=\pi^2+\sum_jx_j^2\ge\pi^2, hence max⁡j∣yj∣≥π/3\max_j\lvert y_j\rvert\ge\pi/\sqrt3, with equality at p=(π,iπ/3,iπ/3,iπ/3)p=(\pi,i\pi/\sqrt3,i\pi/\sqrt3,i\pi/\sqrt3). At this point P(p)=πP(p)=\pi and P(p−2πe1)=−πP(p-2\pi e_1)=-\pi, which contradicts P(p)=P(p−2πe1)P(p)=P(p-2\pi e_1). A strip of smaller height does not contain such a pair of points, and the argument gives nothing there.

Corollary 3 (no analytic Bloch symbol satisfies the hypotheses of Theorem 1). There is no lattice Hamiltonian or Floquet model with a finite-dimensional, real-analytic, periodic Bloch symbol whose graviton band satisfies (B) and (S1)–(S4).

Proof. Theorem 1 makes the band an exact cone on a punctured ball, which is a non-empty open set not containing 00; in the original axes it is Ω2=Q\Omega^2=Q with QQ positive definite. Theorem 2 forbids this for any such symbol at any lattice spacing.

The corollary concerns exact soft universality in the sense of (S1)–(S4). Section 5 measures how close a band can come to satisfying the necessary condition (3).

5. The order of the defect

The condition (3) is necessary for a gauge-invariant soft factor of the form (S1). For a band satisfying (B) and a coupling that is admissible in the following sense, we measure how far a band is from satisfying it. A coupling is admissible if β=∑m≥1βm\beta=\sum_{m\ge1}\beta_m with βm\beta_m real-analytic on R3∖{0}\mathbb R^3\setminus\{0\} and homogeneous of degree mm, and β1=γ∣q∣\beta_1=\gamma\lvert q\rvert with γ≠0\gamma\ne0, which is the lowest-order form that Step 2 of Theorem 1 forces. The defect of a hard quartet q∈Mρq\in\mathcal M_\rho is

Dβ(ρ,q)=∑iσi β(ρqi) (1,∇Ω(ρqi)).(5)\mathcal D_\beta(\rho,q)=\sum_i\sigma_i\,\beta(\rho q_i)\,\big(1,\nabla\Omega(\rho q_i)\big). \tag{5}

It vanishes identically for the exact cone with β=γ∣k∣\beta=\gamma\lvert k\rvert, and it is the failure of the pole term of the soft factor to be gauge invariant.

Proposition 4 (order of the defect). Assume (B) and let Ω2=∣k∣2+Rd+O(∣k∣d+1)\Omega^2=\lvert k\rvert^2+R_d+O(\lvert k\rvert^{d+1}) with Rd≢0R_d\not\equiv0 homogeneous of degree d≥3d\ge3.

(i) For β=γΩ\beta=\gamma\Omega the energy component of Dβ\mathcal D_\beta vanishes and its momentum components are γ2ρd−1∑iσi∂jRd(qi)+O(ρd)\tfrac\gamma2\rho^{d-1}\sum_i\sigma_i\partial_jR_d(q_i)+O(\rho^d).

(ii) Let O⊂M0∗\mathcal O\subset\mathcal M_0^\ast be non-empty and relatively open. For no admissible β\beta is Dβ(ρ,qρ)=o(ρd−1)\mathcal D_\beta(\rho,q_\rho)=o(\rho^{d-1}) for all families qρ∈Mρq_\rho\in\mathcal M_\rho with qρ→q∈Oq_\rho\to q\in\mathcal O.

Relative to ∑i∣ρqi∣∝ρ\sum_i\lvert\rho q_i\rvert\propto\rho, the defect is therefore of order ρd−2\rho^{d-2}, and no admissible β\beta gives a smaller one. We call r=d−2r=d-2 the order of the defect. For a separable band Ω2=∑j=13fw(kj)\Omega^2=\sum_{j=1}^3f_w(k_j) built from the order-2w2w central-difference function fwf_w for k2k^2, d=2w+2d=2w+2 and r=2wr=2w.

Proof. (i) The energy component is ∑iσiγΩ(ρqi)=0\sum_i\sigma_i\gamma\Omega(\rho q_i)=0 on the energy shell. The momentum components are γ2∑iσi∂jΩ2(ρqi)=γ2∑iσi∂jR(ρqi)\tfrac\gamma2\sum_i\sigma_i\partial_j\Omega^2(\rho q_i)=\tfrac\gamma2\sum_i\sigma_i\partial_jR(\rho q_i) by ∑iσiki=0\sum_i\sigma_ik_i=0, and the leading term of ∂jR\partial_jR is ∂jRd\partial_jR_d.

(ii) Suppose Dβ=o(ρd−1)\mathcal D_\beta=o(\rho^{d-1}) along all such families. By Lemma B(iii) every q∈Oq\in\mathcal O is the limit of such a family. The energy component of Dβ\mathcal D_\beta equals ∑iσig(ρqi)\sum_i\sigma_ig(\rho q_i) for g=β−γΩg=\beta-\gamma\Omega, as in Step 3 of Theorem 1. If gmg_m is the lowest non-zero term of gg and m≤d−1m\le d-1, then ∑iσigm=0\sum_i\sigma_ig_m=0 on O\mathcal O, and Lemma B(ii) gives gm=0g_m=0. Hence βm=γΩm\beta_m=\gamma\Omega_m for m≤d−1m\le d-1. Here Ωm=0\Omega_m=0 for 2≤m≤d−22\le m\le d-2 and Ωd−1=Rd/(2∣k∣)\Omega_{d-1}=R_d/(2\lvert k\rvert). The degree-(d−1)(d-1) part of β ∂jΩ\beta\,\partial_j\Omega is then γ∣q∣∂jΩd−1+γΩd−1q^j=γ2∂jRd\gamma\lvert q\rvert\partial_j\Omega_{d-1}+\gamma\Omega_{d-1}\hat q_j=\tfrac\gamma2\partial_jR_d, and the lower parts are either conserved exactly (γqj\gamma q_j in degree one) or zero. So the momentum components of Dβ\mathcal D_\beta are γ2ρd−1∑iσi∂jRd(qi)+O(ρd)\tfrac\gamma2\rho^{d-1}\sum_i\sigma_i\partial_jR_d(q_i)+O(\rho^d), and o(ρd−1)o(\rho^{d-1}) would make ∂jRd\partial_jR_d a collision invariant of degree d−1≥2d-1\ge2 on O\mathcal O. By Lemma B(ii) it vanishes, hence Rd=0R_d=0, a contradiction.

For the separable bands we use fw(k)=2∑ι=1w(2sin⁡(k/2))2ιι2(2ιι)f_w(k)=2\sum_{\iota=1}^{w}\dfrac{(2\sin(k/2))^{2\iota}}{\iota^2\binom{2\iota}{\iota}}, the truncation at ι=w\iota=w of the series k2=2∑ι≥1(2sin⁡(k/2))2ι/(ι2(2ιι))k^2=2\sum_{\iota\ge1}(2\sin(k/2))^{2\iota}/(\iota^2\binom{2\iota}{\iota}) (the series of (2arcsin⁡y)2(2\arcsin y)^2 with y=sin⁡(k/2)y=\sin(k/2)). It is a trigonometric polynomial of degree ww and for w=1,2,3w=1,2,3 it equals 2−2cos⁡k2-2\cos k, 52−83cos⁡k+16cos⁡2k\tfrac52-\tfrac83\cos k+\tfrac16\cos2k and 4918−3cos⁡k+310cos⁡2k−145cos⁡3k\tfrac{49}{18}-3\cos k+\tfrac3{10}\cos2k-\tfrac1{45}\cos3k. The omitted terms (ι>w\iota>w) give fw(k)−k2=cwk2w+2+O(k2w+4)f_w(k)-k^2=c_wk^{2w+2}+O(k^{2w+4}) with cw=−2(w!)2/(2w+2)!c_w=-2(w!)^2/(2w+2)!, which is non-zero for every ww; for w=1,2,3w=1,2,3 it is −112-\tfrac1{12}, −190-\tfrac1{90}, −1560-\tfrac1{560}. So R2w+2=cw∑jkj2w+2≢0R_{2w+2}=c_w\sum_jk_j^{2w+2}\not\equiv0.

Proposition 4 concerns the algebraic condition (3). The condition can be met to relative order ρd−2\rho^{d-2} by a band whose Ω2−∣k∣2\Omega^2-\lvert k\rvert^2 starts at degree dd, with β=γΩ\beta=\gamma\Omega, and to no higher order by any admissible β\beta. No amplitude is constructed whose soft factor is gauge invariant to that order, since (S1) and (S2) are hypotheses of exact invariance.

Example 1 (a spin-ss walk with d=4d=4). Let Σj\Sigma_j be the Pauli matrices, Rj(t)=cos⁡t−isin⁡t Σj\mathsf R_j(t)=\cos t-i\sin t\,\Sigma_j and

V(k)=Rx(kx/2) Ry(ky/2) Rz(kz) Ry(ky/2) Rx(kx/2)∈SU(2),V(k)=\mathsf R_x(k_x/2)\,\mathsf R_y(k_y/2)\,\mathsf R_z(k_z)\,\mathsf R_y(k_y/2)\,\mathsf R_x(k_x/2)\in SU(2),

which is 2π2\pi-periodic in each kjk_j. Write V=cos⁡ϑ−isin⁡ϑ n^⋅ΣV=\cos\vartheta-i\sin\vartheta\,\hat n\cdot\Sigma. A direct computation gives cos⁡ϑ=cos⁡kxcos⁡kycos⁡kz\cos\vartheta=\cos k_x\cos k_y\cos k_z. Let D(s)D^{(s)} be the (2s+1)(2s+1)-dimensional irreducible representation of SU(2)SU(2) and take the one-step symbol W(k)=D(s)(V(k/2s))W(k)=D^{(s)}(V(k/2s)), which acts on 2s+12s+1 internal states per cell. Since D(s)(V)=e−2iϑ n^⋅SD^{(s)}(V)=e^{-2i\vartheta\,\hat n\cdot S}, its eigenphases are 2ςϑ(k/2s)2\varsigma\vartheta(k/2s) for ς=−s,…,s\varsigma=-s,\dots,s, and the helicity-ss band is

Ωs(k)=2s arccos⁡ ⁣(∏j=13cos⁡kj2s),Ωs2=∣k∣2−112s2∑i<jki2kj2+O(k6).(6)\Omega_s(k)=2s\,\arccos\!\Big(\prod_{j=1}^{3}\cos\frac{k_j}{2s}\Big),\qquad \Omega_s^2=\lvert k\rvert^2-\frac1{12s^2}\sum_{i<j}k_i^2k_j^2+O(k^6). \tag{6}

The symbol WW has period 4πs4\pi s in each kjk_j, so in these units the lattice spacing is a=1/(2s)a=1/(2s) and the quartic coefficient is −a2/3-a^2/3; for s=2s=2 it is −1/48-1/48. The band is not a cone for any ss (Theorem 2), R4=−112s2∑i<jki2kj2≢0R_4=-\frac1{12s^2}\sum_{i<j}k_i^2k_j^2\not\equiv0, and Proposition 4 gives d=4d=4, r=2r=2.

Remark (equivalence-principle estimate; not a theorem). Weinberg’s theorem and its S-matrix sharpenings assume an exactly Poincaré-covariant S-matrix, and no result is known that forces exact universality from an ultraviolet theory with a preferred frame. Suppose that the S-matrix is covariant up to relative corrections of order (Ea)r(Ea)^r and that universality holds up to the same corrections. Suppose further that the relative violation of the equivalence principle for nucleon-bound matter is δEP∼(pFa)r\delta_{\rm EP}\sim(p_Fa)^r, with the nucleon Fermi momentum pF≈0.25p_F\approx0.25 GeV and a composition-dependent coefficient of order one. The MICROSCOPE experiment found for the Eötvös parameter ηE(Ti,Pt)=[−1.5±2.3 (stat)±1.5 (syst)]×10−15\eta_{\rm E}(\mathrm{Ti},\mathrm{Pt})=[-1.5\pm2.3\,(\mathrm{stat})\pm1.5\,(\mathrm{syst})]\times10^{-15}, the statistical errors at 1σ1\sigma [Touboul et al.]. Adding the two uncertainties in quadrature and taking the central value plus two standard deviations gives ∣ηE∣<7×10−15\lvert\eta_{\rm E}\rvert<7\times10^{-15}. Setting δEP\delta_{\rm EP} equal to this bound gives 1/a≳pF/(7×10−15)1/r=3×1061/a\gtrsim p_F/(7\times10^{-15})^{1/r}=3\times10^6 GeV for r=2r=2, 9×1029\times10^2 GeV for r=4r=4 and 5757 GeV for r=6r=6. The scaling δEP∼(pFa)r\delta_{\rm EP}\sim(p_Fa)^r is dimensional. The one-loop shifts of the speed of light in the band-limited model of A conserved lattice momentum and graviton scattering are a separate cost of a preferred frame and do not enter this estimate.

6. The Weinberg–Witten theorem and a local lattice energy density

The Weinberg–Witten theorem could give a second obstruction. Its hypothesis is a Lorentz-covariant, conserved, local energy-momentum tensor Θμν\Theta^{\mu\nu} with the Poincaré charges as its integrals. The standard argument for a massless particle of helicity ss is the following sketch. Take the momenta pp and −p-p along an axis, where the one-particle states have angular momentum ss and −s-s about it. A rank-2 tensor operator changes the angular momentum about that axis by at most 2 in absolute value, so ⟨−p,s∣Θμν∣p,s⟩=0\langle -p,s\vert \Theta^{\mu\nu}\vert p,s\rangle=0 if ∣2s∣>2\lvert2s\rvert>2. Boosts map any pair of non-parallel massless momenta to a back-to-back pair, so covariance gives ⟨p′,s∣Θμν∣p,s⟩=0\langle p',s\vert \Theta^{\mu\nu}\vert p,s\rangle=0 for all non-parallel p,p′p,p' when ∣s∣>1\lvert s\rvert>1. If the matrix element is continuous at p′=pp'=p, this contradicts the normalisation ⟨p,s∣Θ00∣p,s⟩=Ω(p)\langle p,s\vert \Theta^{00}\vert p,s\rangle=\Omega(p). We call the vanishing for non-parallel momenta the covariant benchmark.

On a lattice the energy density is a local operator that need not be covariant, so the benchmark need not hold. We measure the departure from it in a model. The one-particle space at small wave vector is the spin-ss irreducible representation of SU(2)SU(2), and the helicity-ss band has the eigenvector u(k)=∣n^(k)⟩u(k)=\lvert\hat n(k)\rangle, the spin coherent state along n^(k)=k^+O(∣k∣)\hat n(k)=\hat k+O(\lvert k\rvert). A graviton has only the helicities ±2\pm2. The model embeds the helicity-ss band in a spin-ss multiplet, as the automaton of Example 1 does, and a statement about it does not cover every possible graviton. The matrix element T(k′,k)T(k',k) below is the one-particle matrix element of the lattice energy density τ\tau, the counterpart of ⟨p′,s∣Θ00∣p,s⟩\langle p',s\vert\Theta^{00}\vert p,s\rangle in the covariant case.

Proposition 5 (a local energy density is not covariant on helicity-ss states). In the model above, let the energy density be a bilinear form with kernel τ(k′,k)\tau(k',k), real-analytic at (0,0)(0,0) (as for an exponentially local lattice operator). Let T(k′,k)=u(k′)†τ(k′,k)u(k)T(k',k)=u(k')^\dagger\tau(k',k)u(k) and assume that on the diagonal it is the band energy up to a relative correction that vanishes with kk, T(k,k)=Ω(k) (1+o(1))T(k,k)=\Omega(k)\,(1+o(1)) with Ω(k)=∣k∣(1+o(1))\Omega(k)=\lvert k\rvert(1+o(1)). Then for k=ρn^k=\rho\hat n, k′=ρn^′k'=\rho\hat n' the limit of ∣T(k′,k)∣/Ω(k)\lvert T(k',k)\rvert/\Omega(k) as ρ→0\rho\to0 exists and equals ∣G(n^′,n^)∣\lvert G(\hat n',\hat n)\rvert, where

G(n^′,n^)=⟨n^′∣L(n^′,n^)∣n^⟩.G(\hat n',\hat n)=\langle\hat n'\vert L(\hat n',\hat n)\vert\hat n\rangle .

Here L(k′,k)L(k',k), the part of τ\tau at the origin that is linear in the pair (k′,k)(k',k), is an operator on the spin-ss space, and GG is its value at k′=n^′k'=\hat n', k=n^k=\hat n taken between the coherent states. Moreover ∣G∣\lvert G\rvert is continuous on S2×S2S^2\times S^2, and ∣G(n^,n^)∣=1\lvert G(\hat n,\hat n)\rvert=1. In particular there is θ0>0\theta_0>0, independent of ρ\rho, with ∣T(k′,k)∣≥12Ω(k)\lvert T(k',k)\rvert\ge\tfrac12\Omega(k) for all small ρ\rho and all directions with ∣n^′−n^∣<θ0\lvert\hat n'-\hat n\rvert<\theta_0.

Proof. Expand τ(k′,k)=τ0+L(k′,k)+O(ρ2)\tau(k',k)=\tau_0+L(k',k)+O(\rho^2) with LL linear. On the diagonal, T(k,k)=⟨n^∣τ0∣n^⟩+O(ρ)T(k,k)=\langle\hat n\vert\tau_0\vert\hat n\rangle+O(\rho) must be O(ρ)O(\rho), so ⟨n^∣τ0∣n^⟩=0\langle\hat n\vert\tau_0\vert\hat n\rangle=0 for every unit n^\hat n. The diagonal coherent-state symbol is injective on operators on the spin-ss representation (it is a polynomial in the components of n^\hat n that determines the matrix), so τ0=0\tau_0=0. Since L(ρn^′,ρn^)=ρL(n^′,n^)L(\rho\hat n',\rho\hat n)=\rho L(\hat n',\hat n), T(k′,k)=ρ ⟨n^′∣L(n^′,n^)∣n^⟩+O(ρ2)T(k',k)=\rho\,\langle\hat n'\vert L(\hat n',\hat n)\vert\hat n\rangle+O(\rho^2) up to the arbitrary phases of u(k),u(k′)u(k),u(k'). The modulus of GG does not depend on these phases and is continuous, and ∣G(n^,n^)∣=lim⁡T(k,k)/Ω(k)=1\lvert G(\hat n,\hat n)\rvert=\lim T(k,k)/\Omega(k)=1. The last statement follows from the uniform continuity of ∣G∣\lvert G\rvert on S2×S2S^2\times S^2 and the uniformity of the O(ρ2)O(\rho^2) remainder.

For s≥32s\ge\frac32 and non-parallel momenta k,k′k,k' of equal length with a small angle, the covariant benchmark is T=0T=0 and the proposition gives ∣T∣≥12Ω\lvert T\rvert\ge\tfrac12\Omega. The hypothesis s≥32s\ge\tfrac32 enters only in this comparison and not in the proof. Within the model, a local analytic energy density is therefore not covariant on helicity-ss states, not even approximately: the violation is of relative order one and does not decrease with the scale, so the Weinberg–Witten hypothesis is not met. In the continuum a free massless spin-2 field likewise has no gauge-invariant, covariant, conserved rank-2 energy-momentum tensor; its covariant conserved bilinear is the Bel–Robinson tensor, of rank four.

An explicit family has the one-particle generator Ω^(k)=Ω(k)s n^(k)⋅S\hat\Omega(k)=\frac{\Omega(k)}{s}\,\hat n(k)\cdot S in the spin-ss representation, whose helicity-ss eigenvector is u(k)=∣n^(k)⟩u(k)=\lvert\hat n(k)\rangle with eigenvalue Ω(k)\Omega(k). Take the energy density τ(k′,k)=12[ψ(Ω^(k′))+ψ(Ω^(k))]\tau(k',k)=\tfrac12[\psi(\hat\Omega(k'))+\psi(\hat\Omega(k))] for a function ψ\psi. Because u(k)u(k) is an eigenvector of ψ(Ω^(k))\psi(\hat\Omega(k)) with eigenvalue ψ(Ω(k))\psi(\Omega(k)), and the overlap of two spin-ss coherent states is ∣⟨n^′∣n^⟩∣=((1+n^⋅n^′)/2)s=cos⁡2s(θ/2)\lvert\langle\hat n'\vert\hat n\rangle\rvert=((1+\hat n\cdot\hat n')/2)^s=\cos^{2s}(\theta/2) for θ\theta the angle between n^\hat n and n^′\hat n', the matrix element is exactly

∣T(k′,k)∣=12∣E(k)+E(k′)∣(1+n^(k)⋅n^(k′)2)s,E=ψ(Ω).(7)\lvert T(k',k)\rvert=\tfrac12\lvert E(k)+E(k')\rvert\Big(\frac{1+\hat n(k)\cdot\hat n(k')}2\Big)^{s},\qquad E=\psi(\Omega). \tag{7}

If Ω(k)n^(k)=k\Omega(k)\hat n(k)=k and ψ\psi is real-analytic with ψ(0)=0\psi(0)=0, ψ′(0)=1\psi'(0)=1, then Ω^(k)=k⋅S/s\hat\Omega(k)=k\cdot S/s and τ\tau is analytic at the origin, so the hypotheses of Proposition 5 hold. For s=2s=2 and equal energies, (7) gives ∣T∣/E=cos⁡4(θ/2)=0.8705, 0.5625, 0.25, 0.0625\lvert T\rvert/E=\cos^4(\theta/2)=0.8705,\ 0.5625,\ 0.25,\ 0.0625 at θ=π/6,π/3,π/2,2π/3\theta=\pi/6,\pi/3,\pi/2,2\pi/3, independent of ρ\rho. For s=1s=1 it gives cos⁡2(θ/2)\cos^2(\theta/2), the overlap of two helicity-one polarisation vectors. No contradiction arises there, because the photon has a covariant conserved energy-momentum tensor, Maxwell’s, and the Weinberg–Witten argument does not require the matrix element to vanish for ∣s∣≤1\lvert s\rvert\le1. For s=12s=\tfrac12 it gives cos⁡(θ/2)\cos(\theta/2).

7. Lattice amplitudes with the assumed soft expansion

A derivation of the band and the soft expansion assumed in (B) and (S1)–(S4) for an interacting lattice or automaton graviton would settle whether Theorem 1 applies to any such model. The first question is whether the hard legs have simple poles 1/(si⋅K)1/(\mathsf s_i\cdot K) in discrete time.

Cubic vertices built from the linearised curvature are invariant under the linear gauge transformations by themselves and can be finite-range, but they carry no coupling to energy. For them (S4) fails, and there is no Newtonian limit. A coupling JJ that is a matrix in helicity or species space is outside Lemma A, which treats a multiple of the identity. A band that is non-analytic in the direction k^\hat k at the origin, as can occur when helicity bands are degenerate at leading order, violates (B), and the expansion in Theorem 1 is unavailable.

For band-limited vertices with a charge linear on the band, A conserved lattice momentum and graviton scattering shows that a nonlinear gauge transformation which keeps band-limited fields band-limited (requirement (G) there, not derived from gauge invariance) makes the band the whole zone. (S2) is then available for such vertices only through a gauge deformation that does not close on band-limited fields, and whether one exists is open there.

Theorem 2 uses a finite-dimensional fibre. A continuum free massless field regarded as a lattice theory with one unit cell has all bands exactly conical, Ω=∣k+2πn∣\Omega=\lvert k+2\pi n\rvert, with an unbounded fibre. Do bounded infinite-dimensional Hamiltonian fibres, or unitary ones, escape the theorem?

The estimate of section 5 is dimensional, and how the order rr of the defect translates into observable violations of the equivalence principle depends on the dynamics. It is open whether the non-covariant energy density of section 6 can be the source to which a cubic graviton vertex couples with a soft limit of defect order rr, as the gravitational pseudotensor does in general relativity.

References

  1. Weinberg, “Photons and gravitons in S-matrix theory: derivation of charge conservation and equality of gravitational and inertial mass,” Phys. Rev. 135, B1049 (1964).
  2. Weinberg, “Infrared photons and gravitons,” Phys. Rev. 140, B516 (1965).
  3. Weinberg, Witten, “Limits on massless particles,” Phys. Lett. B 96, 59 (1980).
  4. Porrati, “Universal limits on massless high-spin particles,” Phys. Rev. D 78, 065016 (2008), arXiv:0804.4672.
  5. Benincasa, Cachazo, “Consistency conditions on the S-matrix of massless particles,” arXiv:0705.4305.
  6. Kato, Sakamoto, So, “No-go theorem of Leibniz rule and supersymmetry on the lattice,” PoS LATTICE2008, 233, arXiv:0810.2360.
  7. Touboul et al. (MICROSCOPE collaboration), “MICROSCOPE mission: final results of the test of the Equivalence Principle,” Phys. Rev. Lett. 129, 121102 (2022), arXiv:2209.15487.