Adversarial review: improved-action result
FATAL: the
T(v) = -S(iv)rule is not a property of improvement; it is a spacetime-lattice assumption.The result states that Beane's eqs. (16) and (17) can be put in the form
T(bE) = sum_j S(b k_j)withT(v) = -S(iv), then applies the same rule to Symanzik/Naik improvement and calls the angular function universal for any tree-level improved action. That is only legitimate for an isotropic spacetime discretisation whose temporal operator is the analytic continuation of the spatial operator. It is not legitimate for a Hamiltonian lattice with continuous time and improved space only, which Beane explicitly notes is possible: "Hamiltonian lattice formulations... are also possible" (arXiv:1210.1847, Introduction, footnote 3 in the fetched HTML).Counter-calculation: for a space-only lattice, the massless dispersion is
E^2 = sum_j S(b k_j) / b^2, notT(bE) = sum_j S(bk_j). IfS(u) = u^2 + s6 u^6 + ...withs4 = 0, thenE = k sqrt(1 + s6 P6 (bk)^4 + ...) = k[1 + (s6/2) P6 (bk)^4 + ...]not
k[1 + (s6/2)(P6 - 1)(bk)^4]. Therefore, for space-only Symanzik/Naik:D_gamma = -P6/180 D_e = -3 P6/40 D_e/D_gamma = 27/2The ratio survives, but the sign and angular factor do not. Both photon and electron remain subluminal. The on-axis
O(b^4)term does not vanish; it is largest on-axis. Numerically from the same threshold formula using the LHAASO low edge, the space-only improved body-diagonal bound is1/b > 1.09e10 GeV, not9.08e9 GeV; on-axis it is1.89e10 GeV, not "no O(b^4) effect." The corrected claim should be: for a spacetime-isotropic improved lattice with the same analytically continued temporal stencil,D_gamma = (1-P6)/180andD_e = 3(1-P6)/40; for a continuous-time space-only lattice,D_gamma = -P6/180andD_e = -3P6/40.SERIOUS: the sign flip is real only under the spacetime analytic-continuation convention.
I independently expanded the stated stencils. For the spacetime rule, the code's sign is correct:
S_Symanzik(u) = u^2 - u^6/90 + u^8/1008 + ... T(v) = -S(iv) = v^2 - v^6/90 + ... d2 = (s6/2)(P6 - 1) = (1-P6)/180 >= 0The Naik case similarly gives
d2 = 3(1-P6)/40 >= 0. So there is no arithmetic sign error indisp.pyfor that model.But the prose says "Improvement flips the sign" as if that follows from tree-level improvement itself. It does not. In the space-only case above,
d2 = s6 P6/2 < 0; improvement cancels theO(b^2)term but the leadingO(b^4)term is still subluminal. Corrected claim: sign flip is a feature of the specific spacetime-improved dispersion relation, not of Symanzik/Naik improvement in general.MINOR: the Symanzik and Naik coefficients check out, but the notation needs to stay explicit.
I attacked the two coefficients and did not break them. The 5-point positive Laplacian is
S(u) = 5/2 - (8/3) cos u + (1/6) cos 2u = u^2 - u^6/90 + u^8/1008 + ...The Naik momentum-space derivative is
F(u) = (9/8) sin u - (1/24) sin 3u = u - (3/40)u^5 + (1/56)u^7 + ... S(u) = F(u)^2 = u^2 - (3/20)u^6 + (1/28)u^8 + ...Thus
s6_Symanzik = -1/90,s6_Naik = -3/20, and the spacetime-improved ratio is(3/40)/(1/180) = 27/2. The only caution is convention: older prose in the previous result mentioned a Naik-1/8; the momentum-space coefficient indisp.pyis-1/24multiplyingsin(3u). The corrected claim should specify momentum-space derivative coefficient-1/24, because real-space three-link normalisations are easy to confuse.SERIOUS: the "universal angular function for any tree-level improved action" is overclaimed even apart from the time-continuation issue.
The derivation proves a narrower theorem: if the free operator is separable,
sum_j S(bk_j), and the temporal operator is fixed byT(v) = -S(iv), thens4 = 0givesd2 = (s6/2)(P6 - 1). It does not prove universality for arbitrary tree-level improved gauge/fermion actions. More general hypercubic quadratic kernels can contain mixed momentum structures at sixth order; the cubic-invariant content is not exhausted by a one-dimensional stencil coefficient unless separability is assumed.Corrected claim: the angular function
1-P6is universal for this separable Symanzik/Naik stencil family under spacetime analytic continuation. Do not say "any tree-level improved action" unless the general free quadratic kernel has been derived.MINOR: the photon-decay threshold formula holds.
I re-derived it from scratch and it matches the code. With
E_gamma(k) = k + D_g b^n k^(n+1) E_e(xk) = xk + m^2/(2xk) + D_e b^n (xk)^(n+1)energy conservation gives
b^n k^(n+1)[D_g - D_e{x^(n+1)+(1-x)^(n+1)}] >= m^2/[2k x(1-x)] k^(n+2) b^n W(x) >= m^2 W(x) = 2x(1-x)[D_g - D_e S_n(x)] 1/b > k (k/m_e)^(2/n) W_max^(1/n)The positive control also holds: for the naive body-diagonal case, the low-edge LHAASO photon gives
4.700e14 GeV, central gives5.696e14 GeV.The Rubtsov et al. variable change is also right. Their
x'is the asymmetry in[-1,1]; if the momentum fraction isx, thenx' = 2x - 1and1 - x'^2 = 4x(1-x). Their conditionomega_LV(x') <= 2m^2/[k(1-x'^2)]becomes2k x(1-x) omega_LV(x) <= m^2, which is the same condition.SERIOUS: the Bethe-Heitler "max_x omega_LV < 0" criterion is not established by the cited literature as the operative shower bound.
The result says the literature's
omega_LV(1)prescription silently fails because it came from a model whereomega_LVis x-independent. That is false for arXiv:1312.4368. Their eq. (6) explicitly hasomega_LV(x) = -kappa k - g k^3/(4M^2)(1+3x^2) + xi k^3/(2M^2)so it is x-dependent when electron LV is present. They nevertheless use
omega_LV(1)because, in their words, "the cross section is peaked at the maximal asymmetry between the momenta of the pair, x = +/-1, hence the appearance of omega_LV at x=1" (arXiv:1312.4368, text below eq. (7)). The older explicit calculation says the same: in the suppressed regime the cross section is dominated by configurations where one fermion carries most of the energy (|x| approx 1) and then writes the total cross section withomega_LVtaken atx=1(arXiv:1204.5782, around eqs. (35), (40)).The report's
max_x omega_LV < 0is a sufficient condition for every splitting configuration to be in the negative-omega regime. It is not shown to be the necessary condition for an order-of-magnitude shower-suppression bound. Corrected claim: because the lattice omega changes sign, the published endpoint formula cannot be used directly; one must integrate the sign-changing two-sector cross section or derive a new mixed-regime approximation.SERIOUS: "the suppression factor is b-independent, therefore no b information" is wrong as stated.
The sign window is b-independent; the magnitude is not. The condition for the RSS suppressed formula is
|k omega_LV| >> m^2, andomega_LV = b^n k^(n+1) F(x). That condition turns on as b changes. In the spacetime-improved body-diagonal case:D_g = 2/405 rho = 13.5 omega/D_g = 1 - rho[x^5 + (1-x)^5] sign roots: x = 0.4322686, 0.5677314 positive window width = 0.1354627 W_max = 1/2592 at x=1/2 max |omega_negative| / (W_max/k factor) = 160At the photon-decay threshold (
1/b = 9.08e9 GeV), the endpoint negative wings havek|omega|/m^2 approx 160, deep in the RSS suppression regime. The onset scale for endpoint suppression is larger by160^(1/4) = 3.56, i.e. roughly3.2e10 GeV. So there is a b-dependent interval in which the asymmetric wings are progressively suppressed while photon decay is still forbidden.A rough smooth Bethe-Heitler weight
1 - (4/3)x(1-x)assigns only11.6%of the standard x-weight to the central positive window[0.4323,0.5677]. That is not zero, but it is right on top of the "within an order of magnitude" criterion used in arXiv:1312.4368. The corrected claim should be: Bethe-Heitler probably does not give a clean stronger bound without a new differential calculation, but the report has not proven structural incapability or b-independence.SERIOUS: "Bethe-Heitler suppression on a lattice is structurally impossible, naive or improved" is broader than the calculation.
The code tests self-consistent naive and self-consistent Symanzik+Naik cases. It does not test mixed improvement. The result itself admits this at lines 240-244, but the headline table says "Bethe-Heitler suppression on a lattice: structurally impossible, naive or improved." That headline is too broad.
Example: a naive photon sector with an improved or near-luminal electron sector is essentially the RSS subluminal-photon setup, and BH suppression should turn back on. Conversely, an improved photon sector with a naive electron sector is dominated by the electron's
O(b^2)term and belongs to the photon-decay side. Corrected claim: for the two self-consistent cases studied, naive and jointly Symanzik+Naik-improved, photon decay remains the relevant sign channel; mixed actions were not analysed and can change the conclusion.SERIOUS: the comparison to Beane's
1e11 GeVis directionally fair but rhetorically overclean.The fetched Beane text supports the core correction: the abstract says the bound
b^-1 >= 10^11 GeVis "derived from the high-energy cut off of the cosmic ray spectrum," and the body says, "For both the fermions and the bosons, the cut off from the dispersion relation is E^max ~ b^-1. Equating this to the GKZ cut off ... corresponds to a mass scale of b^-1 ~ 10^11 GeV" (arXiv:1210.1847, Section IV). So yes, it is basically a cutoff argument, not a Lorentz-violation threshold measurement.But comparing it to the improved PeV photon-decay number is not perfectly apples-to-apples. Beane's number uses UHE cosmic rays near
10^11 GeV; the improved photon-decay number uses a1.29e6 GeVgamma ray plus an assumed electromagnetic dispersion model. Saying the improved photon result is numerically below Beane is fine. Saying the dispersion route "buys nothing" is stronger than shown unless the same improved-action assumptions are propagated into whatever particle species produced the UHECR cutoff and unless one acceptsE_max ~ 1/bas action-independent enough. Corrected claim: one level of improvement drops the PeV photon-decay bound below Beane's cutoff estimate; the two are different observables and should be compared as constraints on the same hypotheticalb, not as the same kind of measurement.MINOR: the "perfect action" ladder is heuristic, not a derived limit.
The algebra 1/b > k (k/m)^(2/n) W^(1/n) -> k as n -> infinity is correct if W stays order unity and if a finite-order dispersion correction remains the right description. A true perfect action is not just the n -> infinity member of the same coefficient sequence; it can reproduce continuum dispersion throughout the Brillouin zone until the cutoff/nonlocality enters. The corrected claim should call this a scaling intuition for successive irrelevant-operator improvement, not a derived statement about perfect actions.
Attacks tried that failed
Coefficient attack failed. The Symanzik 5-point Laplacian gives
s6 = -1/90; the Naik momentum-space derivative(9/8)sin u - (1/24)sin 3ugivess6 = -3/20; the ratio27/2follows.Spacetime sign-arithmetic attack failed. Under the specific
T(v) = -S(iv)spacetime rule, both improved species are superluminal and the body-diagonal bound1/b = 9.08e9 GeVat the LHAASO low edge is reproduced.Threshold-formula attack failed. The general
1/b > k (k/m_e)^(2/n) W_max^(1/n)formula is correct, and the naive positive control reproduces4.700e14 GeVlow edge /5.696e14 GeVcentral.Rubtsov variable-change attack failed.
x' = 2x - 1and1 - x'^2 = 4x(1-x)are correct, and the photon-decay inequality maps cleanly to the threshold formula.Beane-provenance attack failed in the narrow sense. The report is right that Beane's
1e11 GeVis a cutoff estimate from the high-energy cosmic-ray spectrum, not the same kind of Lorentz-violation threshold bound as the photon-decay calculation. The objection is only to how cleanly the two are compared in prose.
Argus