# Beane-Davoudi-Savage lattice scaling thread

Date: 2026-09-09

Scope: exact quantitative scaling in Silas R. Beane, Zohreh Davoudi, Martin J. Savage, "Constraints on the Universe as a Numerical Simulation," arXiv:1210.1847v2, submitted 2012-10-04, revised 2012-11-09, published Eur. Phys. J. A 50 (2014) 148, DOI 10.1140/epja/i2014-14148-0.

Primary sources read:

- Established: arXiv abstract and metadata: https://arxiv.org/abs/1210.1847
- Established: arXiv PDF: https://arxiv.org/pdf/1210.1847
- Established: arXiv source package: https://arxiv.org/e-print/1210.1847
- Established: journal DOI metadata: https://doi.org/10.1140/epja/i2014-14148-0

Searches performed for follow-up/reanalysis:

- Established: INSPIRE record for arXiv:1210.1847, recid 1189720, citation_count 13: https://inspirehep.net/literature/1189720
- Established: INSPIRE `refersto:recid:1189720` citation list; 13 records returned.
- Established: OpenAlex DOI record W2021145183; 54 citing works returned by `filter=cites:W2021145183`.
- Established: Crossref DOI metadata; `is-referenced-by-count` 22.
- Established: arXiv API title/keyword searches for `1210.1847`, title phrase, `cubic symmetry` + `cosmic rays`, `rotational symmetry breaking` + `cosmic rays` + `lattice`, and `simulation hypothesis` + `cosmic rays`.
- Established: Web searches for exact Beane title plus Auger/cubic anisotropy, `Kubic harmonic`, `Y_4^0`/`Y_4^4`/GZK, `octahedral cosmic ray anisotropy Auger`, and Telescope Array spherical harmonics.
- Established with limitation: Google Scholar query for exact title was fetched, but the automatable page exposed only a normal search result page and not a usable citation list. I did not rely on it for citation enumeration.

## 1. Exact chain to `b^{-1} \gtrsim 10^{11} GeV`

Established: The paper's abstract states the headline result: "Among the observables that are considered are the muon $g-2$ and the current differences between determinations of $\alpha$, but the most stringent bound on the inverse lattice spacing of the universe, $b^{-1}\gsim 10^{11}~{\rm GeV}$, is derived from the high-energy cut off of the cosmic ray spectrum."

Established: The relevant section is Section IV, "ROTATIONAL SYMMETRY BREAKING," subsection 4, "The Energy-Momentum Relation and Cosmic Rays." Section I says: "Section IV considers the simplest effects of Lorentz symmetry breaking operators that first appear at ${\cal O}(\lattspace^2)$, and modifications to the energy-momentum relation. Constraints on the energy-momentum relation due to cosmic ray events are found to provide the most stringent bound on $\lattspace$."

Established: The mathematical input is the lattice dispersion relation. The paper gives boson and Wilson-fermion dispersion relations as equations (16) and (17):

> "The dispersion relations satisfied by bosons and Wilson fermions in a lattice simulation (in Minkowski space) are" followed by equation (16) for bosons and equation (17) for Wilson fermions.

Equation (16), label `dispersionrelationsb`:

```tex
\sinh^2({b E_b\over 2}) - \sum_{j=1,2,3} \sin^2({b k_j\over 2}) - ({b m_b\over 2})^2 = 0 ;
E_b = \sqrt{|{\bf k}|^2+m_b^2} + {\cal O}(b^2) .
```

Equation (17), label `dispersionrelationsf`:

```tex
\sinh^2(b E_f)-\sum_{j=1,2,3} \sin^2(b k_j) - [ b m_f + 2 r ( \sum_{j=1,2,3} \sin^2({b k_j\over 2}) - \sinh^2({b E_f\over 2}) )]^2 = 0 ;
E_f = \sqrt{|{\bf k}|^2+m_f^2} - {r b m_f^3 \over 2 \sqrt{|{\bf k}|^2+m_f^2}} + {\cal O}(b^2) .
```

Established: Their bound comes from the maximum energy implied by those lattice dispersion relations, not from fitting a measured anisotropy amplitude. Exact quote:

> "The lattice spacing itself introduces a cut off to the cosmic ray spectrum. For both the fermions and the bosons, the cut off from the dispersion relation is $E^{\rm max} \sim 1/ \lattspace$. Equating this to the GKZ cut off corresponds to a lattice spacing of $b\sim 10^{-12}~{\rm fm}$, or a mass scale of $b^{-1}\sim 10^{11}~{\rm GeV}$. Therefore, the lattice spacing used in the lattice simulation of the universe must be $b \lsim 10^{-12}~{\rm fm}$ in order for the GZK cut off to be present or for the lattice spacing itself to provide the cut off in the cosmic ray spectrum."

Established: Dimensional conversion of their stated value is consistent at order-of-magnitude level: `1 fm = 5.0677 GeV^{-1}`, so `b = 10^{-12} fm = 5.0677e-12 GeV^{-1}`, and `b^{-1} = 1.97e11 GeV`, quoted by them as `~10^11 GeV`.

Established: The experimental input they cite is not a single event but observed suppression/decline in the UHECR flux near the GZK region. They write:

> "Processes such as $\gamma_{\rm CMB}+N\rightarrow\Delta$ give rise to the predicted GKZ-cut off scale [50, 51] of $\sim 6\times 10^{20}~{\rm eV}$ in the spectrum of high energy cosmic rays. Recent experimental observations show a decline in the fluxes starting around this value [52, 53], indicating that the GKZ-cut off (or some other cut off mechanism) is present in the cosmic ray flux."

Established: References [52] and [53] are:

- Pierre Auger Collaboration, J. Abraham et al., "Measurement of the energy spectrum of cosmic rays above 10^18 eV using the Pierre Auger Observatory," Phys. Lett. B 685, 239 (2010), arXiv:1002.1975, DOI 10.1016/j.physletb.2010.02.013. URL: https://arxiv.org/abs/1002.1975
- P. Sokolsky for the HiRes Collaboration, "Final Results from the High Resolution Fly's Eye (HiRes) Experiment," PoS ICHEP2010, 444 (2010), arXiv:1010.2690. URL: https://arxiv.org/abs/1010.2690

Established: The Auger paper's abstract says: "Above the ankle the spectrum is described by a power law with index 2.6 followed by a flux suppression, above about log10(E/eV) = 19.5, detected with high statistical significance." Its summary gives the more quantitative fit: "In comparison to the power law extrapolation, the spectrum is suppressed by a factor two at log10(E1/2/eV) = 19.61 +/- 0.03. The significance of the suppression is larger than 20 sigma."

Established: The HiRes conference paper abstract says: "We observe a cutoff consistent with the GZK predictions with a five sigma significance." In the body it says the GZK prediction is "a cutoff near 6 x 10^19 eV" and "The statistical strength of the monocular observation of the GZK cutoff is 5.3 sigma."

Own inference: There is a source-level numerical looseness here. Beane et al. quote the predicted cutoff scale as `~6 x 10^20 eV`, but the Auger fit they cite gives `log10(E1/2/eV)=19.61`, i.e. `~4.1 x 10^19 eV`, and HiRes text says `~6 x 10^19 eV`. Their final `b^{-1} ~ 10^11 GeV = 10^20 eV` should be read as an order-of-magnitude scale, not a precision bound from the Auger/HiRes measurement.

## 2. Quantitative relationship between `b` and UHECR anisotropy

Established: There is no explicit formula in Beane et al. that maps lattice spacing `b` to a fractional arrival-direction anisotropy amplitude, event-count excess, angular power-spectrum coefficient, or Auger/TA observable. The arrival-direction statement is qualitative.

Established: The only explicit angular formula is equation (18), label `eq:gkzkins`, for the direction-dependent photon energy threshold for head-on `gamma_CMB + N -> Delta` interactions when the lattice rest frame equals the CMB rest frame:

```tex
\omega = {m_\Delta^2-m_N^2\over 4|{\bf p}|}
\left[ 1 + { \sqrt{\pi}\lattspace^2 |{\bf p}|^2\over 9}
\left( Y_4^0(\theta,\phi) +
\sqrt{5\over 14}\left( Y_4^{+4}(\theta,\phi) + Y_4^{-4}(\theta,\phi)\right)\right)
\right]
- {m_\Delta^3-m_N^3\over 4|{\bf p}|} b r + ...
```

Established: Immediately after equation (18), they write:

> "for $|{\bf p}|\ll 1/\lattspace$, where $\theta$ and $\phi$ are the polar and azimuthal angles of the particle momenta in the rest frame of the lattice, respectively. This represents a lower bound for the energy of photons participating in such a process with arbitrary collision angles."

Established: The fractional direction-dependent shift in this threshold, as displayed in equation (18), is:

```tex
delta omega / omega_0 = (sqrt(pi)/9) b^2 |p|^2 [Y_4^0(theta,phi) + sqrt(5/14)(Y_4^{+4}(theta,phi)+Y_4^{-4}(theta,phi))]
```

where `omega_0 = (m_Delta^2 - m_N^2)/(4|p|)`. This is a threshold anisotropy, not an arrival-direction flux-amplitude formula.

Established: Their qualitative arrival-direction language is:

> "For lattice spacings corresponding to an energy scale comparable to the GKZ cut off, the cosmic ray spectrum will exhibit significant deviations from isotropy, revealing the cubic structure of the lattice. However, for lattice spacings much smaller than the GKZ cut off scale, the GKZ mechanism cuts off the spectrum, effectively hiding the underlying lattice structure."

and:

> "The most striking feature of the scenario in which the lattice provides the cut off to the cosmic ray spectrum is that the angular distribution of the highest energy components would exhibit cubic symmetry in the rest frame of the lattice, deviating significantly from isotropy. For smaller lattice spacings, the cubic distribution would be less significant, and the GKZ mechanism would increasingly dominate the high energy structure."

Own inference: Equation (18) implies the threshold direction-dependence scales as `(b |p|)^2` times a pure `l=4` cubic harmonic, plus an isotropic Wilson-fermion term linear in `b r`. But Beane et al. do not propagate that threshold shift through a source distribution, propagation model, composition model, detector exposure, or magnetic deflection model to get a measurable fractional anisotropy amplitude.

## 3. Multipole/cubic/octahedral structure

Established: The paper says "cubic" and "hyper-cubic" repeatedly, but it does not use `O_h`, `octahedral`, or `octahedral symmetry` anywhere in the LaTeX source. It does not explicitly say which spherical-harmonic multipoles survive under `O_h`.

Established: The exact symmetry statements are:

> "While the former type of modifications could arise from many different BSM scenarios, the latter, particularly modifications that exhibit cubic symmetry, would be suggestive of a structure consistent with an underlying discretization of space-time."

> "At ${\cal O}(\lattspace^2)$ in the lattice spacing expansion of the Wilson action, that is relevant to describing low-energy processes, there is a rotational-symmetry breaking operator that is consistent with the lattice hyper-cubic symmetry," followed by equation (14).

> "The violation of Lorentz invariance resulting from these dispersion relations is due to the fact that they have only cubic symmetry and not full rotational symmetry, as shown in fig. 2."

Established: The only explicit spherical-harmonic structure is in equation (18):

```tex
Y_4^0(theta,phi) + sqrt(5/14)(Y_4^{+4}(theta,phi) + Y_4^{-4}(theta,phi))
```

Own inference: This is the standard real cubic harmonic at `l=4`, usually associated with cubic/octahedral symmetry. But that representation-theory identification is not stated in Beane et al.; it is an inference from their equation (18).

Own inference: Because the first displayed angular correction is an `l=4` combination, the paper gives no displayed `l=1`, `l=2`, or `l=3` anisotropy term for the GZK threshold in this approximation. That is not the same as the authors explicitly proving or listing all allowed/suppressed multipoles.

## 4. GZK cutoff direction-dependence: angular scale and amplitude

Established: Their treatment is kinematic and threshold-level. They consider head-on high-energy proton/CMB photon interactions through the `Delta` resonance in the lattice/CMB rest frame and write the direction-dependent threshold equation (18).

Established: Their explicit assumptions in the equation are:

> "When the lattice rest frame coincides with the CMB rest frame, head-on interactions between a high energy proton with momentum $|{\bf p}|$ and a photon of (very-low) energy $\omega$ can proceed through the $\Delta$ resonance when" equation (18) holds.

Established: They do not state an angular scale in degrees. They do not state an amplitude such as `x% anisotropy`. They use qualitative phrases: "significant deviations from isotropy," "deviating significantly from isotropy," and "less significant" for smaller `b`.

Own inference: The angular pattern in equation (18) is `l=4`; a spherical harmonic of degree `l=4` corresponds to a characteristic angular scale of order `180 deg / 4 ~ 45 deg`, but this angular-scale translation is not stated in the paper.

Own inference: From equation (18), the threshold modulation amplitude is order `(b|p|)^2` times the `l=4` cubic harmonic. At `|p| ~ 1/b` that would be order unity, but equation (18) is explicitly written for `|p| << 1/b`, so using it at the cutoff is not controlled.

## 5. Attackable assumptions

Established: Classical-computer simulation assumption. Section I says:

> "In this work, we take a pedestrian approach to the possibility that our universe is a simulation, by assuming that a classical computer (i.e. the classical limit of a quantum computer) is used to simulate the quantum universe (and its classical limit), as is done today on a very small scale, and ask if there are any signatures of this scenario that might be experimentally detectable."

Established: The authors explicitly set aside information-movement/computational-resource constraints:

> "Further, we do not consider the implications of, and constraints upon, the underlying information, and its movement, that are required to perform such extensive simulations."

Established: They focus on a cubic lattice despite unknown future algorithms/hardware:

> "It is the case that the method of simulation, the algorithms, and the hardware that are used in future simulations are unknown, but it is conceivable that some of the ingredients used in present day simulations of quantum fields remain in use, or are used in other universes, and so we focus on one aspect only: the possibility that the simulations of the future employ an underlying cubic lattice structure."

Established: They assume the historical development of universe simulations parallels lattice QCD and starts with cheap unimproved discretizations:

> "Given this low-energy description, we would like to investigate the hypothesis that we are a simulation with the assumption that the development of simulations of the universe in some sense parallels the development of lattice QCD calculations. That is, early simulations use the computationally 'cheapest' discretizations with no improvement."

Established: They assume a hyper-cubic grid and unimproved Wilson action:

> "In particular, we will assume that the simulation of our universe is done on a hyper-cubic grid ... and, as a starting point, we will assume that the simulator is using an unimproved Wilson action, that produces ${\cal O}(b)$ artifacts of the form of the Sheikholeslami-Wohlert operator in the low-energy theory."

Established: They acknowledge improved/chiral-preserving discretizations can evade the `O(b)` precision constraints:

> "For more sophisticated simulations in which chiral symmetry is preserved by the lattice discretization, the coefficient $ {\cal C}_p$ will vanish or will be exponentially small. As a result, the bound on the lattice spacing derived from the muon $g-2$ and from the differences between determinations of $\alpha$ will be significantly weaker."

Established: They assume unresolved future ingredients exist:

> "A number of elements required for a simulation of our universe directly from the fundamental laws of physics have not yet been established, and we have assumed that they will, in fact, be developed at some point in the future; two important elements being an algorithm for simulating chiral gauge theories, and quantum gravity."

Established: They assume composite-particle dispersion follows elementary lattice forms:

> "While for the fundamental particles, the dispersion relations in eq. (16) and eq. (17) are valid, for composite particles, such as the proton or pion, the dispersion relations will be dynamically generated. In the present analysis we assume that the dispersion relations for all particles take the form of those in eq. (16) and eq. (17)."

Established: They assume lattice rest frame equals CMB rest frame for equation (18):

> "When the lattice rest frame coincides with the CMB rest frame..."

Established: They note a more complete proton treatment is required. Footnote 10 says:

> "A more complete treatment of this process involves using the parton distributions of the proton to relate its energy to its momentum... More refined explorations of this and other processes are required."

Serious speculation: A critic can attack each of these: rigid preferred-frame lattice; unbroken cubic symmetry over cosmological propagation distances; no adaptive/dynamical/random lattice; no Lorentz-covariant discretization; no improved action; uncertain UHECR mass composition; proton/Delta-resonance simplification; magnetic deflections; unknown source distribution; using a threshold formula valid for `|p| << 1/b` to motivate behavior near `E ~ 1/b`; and the absence of an explicit detector-level anisotropy prediction.

## 6. Published follow-up, critique, reanalysis, and actual Auger/TA tests

### Direct citation trail

Established: INSPIRE reports 13 citing records for the Beane paper as of this search. Most are broad digital-physics, simulation-hypothesis, or quantum-computation citations, not reanalyses of Beane's cosmic-ray prediction.

Established: OpenAlex reports 54 citing works; Crossref reports 22. The broader OpenAlex list is dominated by philosophy of simulation, ethics/AI, popular-physics, and general computation papers, plus a few computational-physics papers. I found no OpenAlex/Crossref title indicating an Auger/TA cubic anisotropy test of Beane et al.

### Follow-ups/criticisms worth keeping

Serious speculation: Timothy D. Andersen, "Lorentz Covariant Lattice Gauge Theory," arXiv:1210.8348 (2012), directly responds to the Beane-style premise that a lattice creates observable preferred directions. Its abstract says: "Lattice gauge theory's discretization of spacetime suffers from a drawback in that Lorentz covariance is lost because the axes of the lattice create preferred directions in spacetime. Smaller and smaller lattice spacings decrease the effect but fail to eliminate it completely. It has been argued recently that detecting such a set of preferred directions or similar constraints would indicate whether the universe itself has an underlying lattice, i.e. the digital universe hypothesis." It then proposes replacing the lattice by a lattice graph and concludes this "suggests that, even in a digital universe, Lorentz covariance can still hold." URL: https://arxiv.org/abs/1210.8348

Serious speculation: Tom Campbell, Houman Owhadi, Joe Sauvageau, David Watkinson, "On testing the simulation theory," arXiv:1703.00058 (2017), proposes conceptual wave/particle duality tests under a finite-resource/render-on-observation assumption. This is a simulation-test follow-up in the broad literature, not a Beane cosmic-ray reanalysis. Abstract: "Guided by this principle we describe conceptual wave/particle duality experiments aimed at testing the simulation theory." URL: https://arxiv.org/abs/1703.00058

Serious speculation: Zura Kakushadze, "Does the Universe have a Hard Drive?" arXiv:1701.07161 (2017), cites Beane in the broader computational-universe trail but does not test the cosmic-ray anisotropy. Abstract: "We discuss an apparent information paradox that arises in a materialist's description of the Universe if we assume that the Universe is 100% quantum." URL: https://arxiv.org/abs/1701.07161

Serious speculation: Martin Leckey, "Quantum Measurement, Complexity and Discrete Physics," arXiv:quant-ph/0310033, v2 2016, is discrete-physics/modified-QM rather than a Beane reanalysis. Abstract: "This paper presents a new modified quantum mechanics, Critical Complexity Quantum Mechanics, which includes a new account of wavefunction collapse." URL: https://arxiv.org/abs/quant-ph/0310033

Established/Serious speculation: Franco Vazza, "Astrophysical constraints on the simulation hypothesis for this Universe: why it is (nearly) impossible that we live in a simulation," arXiv:2504.08461; Frontiers in Physics DOI 10.3389/fphy.2025.1561873 (2025), is the most direct modern astrophysical critique found. It cites Beane explicitly: "A remarkable exception is the work by \citet{2014EPJA...50..148B}, who investigated the potentially observable consequences of the SH, by exploring the particular case of a cubic space-time lattice." It says Vazza's work uses UHECRs/neutrinos "in a totally different way." Its abstract conclusion is: "In all cases, the amounts of energy or power required by any version of the simulation hypothesis are entirely incompatible with physics, or (literally) astronomically large, even in the lowest resolution case." URL: https://arxiv.org/abs/2504.08461

Anomaly/Anecdote: J. P. Rachen and Ute G. Gahlings, "Conspiratorial cosmology - the case against the Universe," arXiv:1303.7476 (2013), appears in the broader citation trail but is explicitly satirical/popular-physics; not a serious Beane reanalysis. URL: https://arxiv.org/abs/1303.7476

### Has anyone run the Beane cubic-anisotropy test on Auger or Telescope Array data?

Established null result from search: I found no published paper that explicitly runs the Beane-Davoudi-Savage cubic-symmetry (`l=4` cubic harmonic / lattice rest-frame) test on Auger or Telescope Array data. I searched arXiv, INSPIRE, OpenAlex/Crossref citation trails, Google Scholar access, and targeted web queries for combinations of Beane/Davoudi/Savage, `1210.1847`, Auger, Telescope Array, cubic symmetry, octahedral, `O_h`, `Kubic harmonic`, `Y_4^0`, `Y_4^4`, GZK, and UHECR anisotropy.

Established adjacent result: The Pierre Auger and Telescope Array Collaborations did publish a full-sky harmonic analysis above `10^19 eV`: "Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of 10^19 eV at the Pierre Auger Observatory and the Telescope Array," ApJ 794, 172 (2014), arXiv:1409.3128. Its abstract says: "Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage." It reports: "No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole." URL: https://arxiv.org/abs/1409.3128

Established: The 2014 Auger/TA paper is not a Beane test. Its full text contains no `Beane`, `simulation`, `lattice`, `cubic`, `octahedral`, or `Y_4` hit in the searched HTML text. It does show a generic multipole map truncated at `l=4` and says: "Overall, no significant deviation from isotropy is found from this study." But its explicit upper limits are for dipole/quadrupole, not for an oriented cubic `l=4` template.

Established adjacent result: Pierre Auger Collaboration, "Large-scale cosmic ray anisotropies with 19 years of data from the Pierre Auger Observatory," ApJ 976 (2024) 48, arXiv:2408.05292. Its abstract says: "Additionally, the results for the angular power spectrum are shown, demonstrating no other statistically significant multipoles." In Section III.3 it states that Auger's incomplete sky coverage means "the estimation of the individual a_lm coefficients cannot be carried out with relevant resolution as soon as l_max > 2 and the same is true for the power spectrum (full-sky analyses are carried out by the Pierre Auger and Telescope Array Collaborations together...)." URL: https://arxiv.org/abs/2408.05292

Established: The 2024 Auger analysis is also not a direct Beane test. It contains no searched Beane/simulation/lattice/cubic/octahedral/Y4 language. It reports angular-power-spectrum results over multipoles and energy bins, with significant dipolar structure and no significant post-trials higher multipoles: "All other C_l values in different energy bins are not significant." But it does not fit an oriented cubic harmonic or translate a null result into a bound on `b`.

Established adjacent result: Pierre Auger Collaboration, "Observation of a Large-scale Anisotropy in the Arrival Directions of Cosmic Rays above 8 x 10^18 eV," Science 357 (2017) 1266, arXiv:1709.07321. Abstract: "Using 3 x 10^4 cosmic rays above 8 x 10^18 electron volts... The anisotropy, detected at more than the 5.2 sigma level of significance, can be described by a dipole with an amplitude of 6.5_{-0.9}^{+1.3}%..." URL: https://arxiv.org/abs/1709.07321

Established: The 2017 Auger result is astrophysical dipole evidence, not a cubic-lattice signature. It works at lower energies than the Beane cutoff-scale test and does not discuss cubic symmetry.

Bottom line: Beane et al. give a testable qualitative signature and an explicit `l=4` threshold correction, but I found no published direct implementation of the oriented cubic-anisotropy test on Auger/TA data, and no published reanalysis turning Auger/TA null higher-multipole results into a new bound on `b`.