CPT parity check: cubic lattice operators versus SME photon bounds
Thread date: 2026-09-10
Question: does the Beane-Davoudi-Savage cubic-lattice inference evade the strongest vacuum-birefringence bounds because it populates CPT-even, nonbirefringent dimension-6 photon-sector operators rather than CPT-odd dimension-5 birefringent ones?
Sources read
- Established. Silas R. Beane, Zohreh Davoudi, Martin J. Savage, "Constraints on the Universe as a Numerical Simulation," Eur. Phys. J. A 50, 148 (2014), arXiv:1210.1847, https://arxiv.org/abs/1210.1847 and readable HTML https://ar5iv.labs.arxiv.org/html/1210.1847.
- Established. V. Alan Kostelecky and Matthew Mewes, "Electrodynamics with Lorentz-violating operators of arbitrary dimension," Phys. Rev. D 80, 015020 (2009), arXiv:0905.0031, https://arxiv.org/abs/0905.0031 and readable HTML https://ar5iv.labs.arxiv.org/html/0905.0031.
- Established. V. Alan Kostelecky and Neil Russell, "Data Tables for Lorentz and CPT Violation," Rev. Mod. Phys. 83, 11 (2011), arXiv:0801.0287; January 2026 update v19, https://arxiv.org/abs/0801.0287v19 and https://arxiv.org/html/0801.0287v19.
- Established. A. Albert et al. (HAWC Collaboration), "Constraints on Lorentz Invariance Violation from HAWC Observations of Gamma Rays above 100 TeV," Phys. Rev. Lett. 124, 131101 (2020), arXiv:1911.08070, https://arxiv.org/abs/1911.08070.
1. What the cubic lattice produces
Established. Beane, Davoudi, and Savage explicitly frame the regulator as a hypercubic/cubic lattice: "In lattice QCD, space-time is replaced by a finite hyper-cubic grid of points over which the fields are defined" and "The grid breaks Lorentz symmetry (and hence rotational symmetry)" (Beane, Davoudi, Savage 2014, Eur. Phys. J. A 50, 148, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
Established. For a boson on the lattice, their exact dispersion relation is
sinh^2(b E_b / 2) - sum_{j=1,2,3} sin^2(b k_j / 2) - (b m_b / 2)^2 = 0,
with small-momentum expansion
E_b = sqrt(|k|^2 + m_b^2) + O(b^2).
The same section says: "The summations are performed over the components along the lattice Cartesian axes corresponding to the x,y, and z spatial directions" and "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" (Beane, Davoudi, Savage 2014, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
Your own inference, derived from their exact equation. Expanding the boson relation gives
E^2 - |k|^2 - m^2 + (b^2/12)(E^4 + sum_j k_j^4) + O(b^4) = 0,
so, on shell,
E^2 = |k|^2 + m^2 - (b^2/12)[(|k|^2 + m^2)^2 + sum_j k_j^4] + O(b^4).
For a massless photonlike boson this is
E = |k| - (b^2 |k|^3 / 24)(1 + sum_j n_j^4) + O(b^4 |k|^5), n_j = k_j / |k|.
Your own inference. The leading correction is order b^2 |k|^3 in energy, equivalently order b^2 |k|^4 in E^2; in effective-field-theory language this is a dimension-6 correction with coefficient scale ~ b^2. Its angular dependence is the cubic invariant sum_j n_j^4, which decomposes into rotational j=0 and cubic-anisotropic j=4 pieces.
Your own inference. The correction is even under momentum reversal and contains no Levi-Civita/helicity sign, so it is CPT-even rather than CPT-odd. Beane et al. do not state a CPT classification in this section; I looked in their dispersion-relation and Lorentz-violation discussion and found no direct SME/CPT mapping.
Your own inference. The displayed boson dispersion is a scalar relation and contains no +/- helicity branch. It is therefore nonbirefringent at this leading order. Beane et al. do note a finite-lattice polarization artifact: "the polarizations of the massless vector fields are not exactly perpendicular to their direction of propagation for some directions of propagation with respect to the lattice axes, with longitudinal components present for non-zero lattice spacings" (Beane, Davoudi, Savage 2014, arXiv:1210.1847, https://arxiv.org/abs/1210.1847), but that sentence does not by itself establish vacuum birefringence because it does not state a helicity-dependent phase velocity.
2. SME photon-sector classification
Established. Kostelecky and Mewes define the vacuum photon coefficients as four sets: "The net result ... vacuum behavior is controlled by four sets of vacuum coefficients c_(I)jm^(d), k_(E)jm^(d), k_(B)jm^(d), and k_(V)jm^(d)" (Kostelecky and Mewes 2009, Phys. Rev. D 80, 015020, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
Established. Their precise CPT/dimension count is: "For even d, there are (d-1)^2 coefficients c_(I)jm^(d), (d-1)^2 - 4 coefficients k_(E)jm^(d), and (d-1)^2 - 4 coefficients k_(B)jm^(d), giving a total of 3(d-1)^2 - 8 independent components associated with CPT-even Lorentz violation in vacuum propagation. For odd d, the coefficients k_(V)jm^(d) for CPT-odd Lorentz violation have (d-1)^2 independent components" (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
Established. Their birefringence rule is explicit: "The coefficients k_(E)jm^(d) and k_(B)jm^(d) are associated with CPT-even operators that lead to birefringence, with the propagating modes being linearly polarized. The coefficients k_(V)jm^(d) control CPT-odd birefringence, and the corresponding eigenmodes are circularly polarized. In contrast, the CPT-even operators associated with the coefficients c_(I)jm^(d) are nonbirefringent" (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
Established. The same paper says a purely nonbirefringent vacuum model with only c_(I)jm^(d) has "Lorentz-violating operators only in even dimensions d = 4, 6, 8, ..." (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
Established. In minimal notation, k_AF is the CPT-odd photon coefficient and k_F is the CPT-even photon coefficient; at arbitrary dimension the vacuum classes sort into k_(V) for CPT-odd birefringent terms, k_(E) and k_(B) for CPT-even birefringent terms, and c_(I) for CPT-even nonbirefringent terms (Kostelecky and Mewes 2009, Phys. Rev. D 80, 015020, arXiv:0905.0031, https://arxiv.org/abs/0905.0031; Kostelecky and Russell 2011/2026 Data Tables, Rev. Mod. Phys. 83, 11, arXiv:0801.0287v19, https://arxiv.org/abs/0801.0287v19).
3. Where the cubic lattice lands
Your own inference. The Beane-Davoudi-Savage leading photonlike lattice artifact maps to SME c_(I)jm^(6): it is even-dimension, CPT-even, dispersive, direction-dependent, and common to both helicities. The cubic angular factor 1 + sum_j n_j^4 populates an isotropic component plus cubic-anisotropic j=4 components, not the spin-weighted polarization structures k_(E) or k_(B).
Established plus your own inference. Kostelecky-Mewes establish that c_(I)jm^(d) is nonbirefringent; applying that rule to the scalar Beane lattice dispersion makes the leading cubic-lattice photon-sector artifact nonbirefringent (Kostelecky and Mewes 2009, arXiv:0905.0031; Beane, Davoudi, Savage 2014, arXiv:1210.1847).
Not found. I did not find a peer-reviewed paper that explicitly says "the Beane-Davoudi-Savage cubic lattice equals SME c_(I)jm^(6)" with a coefficient-by-coefficient dictionary. The mapping above is a direct classification from the published lattice dispersion and the published SME vacuum taxonomy.
4. Best bounds on nonbirefringent CPT-even dimension-6 photon coefficients
Established. The January 2026 Kostelecky-Russell Data Tables summarize the strongest isotropic sensitivities as k_(V)00^(5): 10^-34 GeV^-1 and c_(I)00^(6): 10^-30 GeV^-2 (Kostelecky and Russell 2011/2026, Rev. Mod. Phys. 83, 11, arXiv:0801.0287v19, Summary Table S3, https://arxiv.org/abs/0801.0287v19). The units differ, so raw exponent comparison is not physically meaningful; as suppression scales they are roughly 1e34 GeV for dimension 5 versus 1e15 GeV for dimension 6.
Established. The strongest listed c_(I)jm^(6) bounds in the 2026 Data Tables come from HAWC high-energy gamma-ray threshold/reaction physics, not a laboratory experiment. Table D22 lists, from A. Albert et al. (HAWC Collaboration), Phys. Rev. Lett. 124, 131101 (2020), arXiv:1911.08070, https://arxiv.org/abs/1911.08070:
| coefficient or direction combination |
bound |
evidence class |
| ` |
c_(I)00^(6) |
` |
| ` |
sum Y_jm(103.45 deg, 276.41 deg)c_(I)jm^(6) |
` |
| ` |
sum Y_jm(83.75 deg, 286.95 deg)c_(I)jm^(6) |
` |
| ` |
sum Y_jm(67.96 deg, 83.6 deg)c_(I)jm^(6) |
` |
| ` |
sum Y_jm(53.26 deg, 304.94 deg)c_(I)jm^(6) |
` |
Established. HAWC's mechanism is explicitly threshold/reaction based: "Superluminal LIV enables the decay of photon at high energy" and "HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources" (Albert et al. 2020, Phys. Rev. Lett. 124, 131101, arXiv:1911.08070, https://arxiv.org/abs/1911.08070). The abstract reports that HAWC excludes the corresponding LIV energy scale to 2.2 x 10^31 eV for the n=1 case, "over 1800 times the Planck energy"; for the dimension-6/nonbirefringent n=2 SME entries, the Data Tables conversion above gives the relevant ~10^-30 GeV^-2 sensitivities (Albert et al. 2020; Kostelecky and Russell 2026 Data Tables).
Established. Other c_(I)jm^(6) bounds listed in Table D22 are much weaker than HAWC but still astrophysical rather than laboratory: Guerrero, Campoy-Ordaz, Potting, and Gaug, Phys. Rev. D 112, 104002 (2025), arXiv:2508.02883, https://arxiv.org/abs/2508.02883, report individual c_(I)jm^(6) sensitivities around 10^-15 GeV^-2; Wei, Liu, Wei, Zhang, and Wu, Universe 8, 519 (2022), arXiv:2210.03897, https://arxiv.org/abs/2210.03897, list |c_(I)00^(6)| = 4.25(-1.63+1.60) x 10^-15 GeV^-2; Agrawal, Singirikonda, and Desai, JCAP 05, 029 (2021), arXiv:2102.11248, https://arxiv.org/abs/2102.11248, list a direction combination at 10^(-14.2 +/- 0.1) GeV^-2; Du et al., Astrophys. J. 906, 8 (2021), arXiv:2010.16029, https://arxiv.org/abs/2010.16029, list |c_(I)00^(6)| < 2.4 x 10^-12 GeV^-2; Rubtsov, Satunin, and Sibiryakov, JCAP 05, 049 (2017), arXiv:1611.10125, https://arxiv.org/abs/1611.10125, list c_(I)00^(6) < 4 x 10^-23 GeV^-2; Vasileiou et al., Phys. Rev. D 87, 122001 (2013), arXiv:1305.3463, https://arxiv.org/abs/1305.3463, list c_(I)00^(6) ranges at 10^-20 GeV^-2 scale; and Fermi-LAT/H.E.S.S./MAGIC time-of-flight bounds are listed in the 10^-19 to 10^-21 GeV^-2 range (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19).
Established. Laboratory/terrestrial bounds on nonminimal c_(I)jm^(6) are far weaker: Kostelecky, Melissinos, and Mewes, Phys. Lett. B 761, 1 (2016), arXiv:1608.02592, https://arxiv.org/abs/1608.02592, are listed with interferometry sensitivities such as |c_(I)21^(6)| < 1.0 x 10^-2 GeV^-2 and |c_(I)22^(6)| < 1.3 x 10^-4 GeV^-2 (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19).
Your own inference. Comparing the best quoted isotropic dimension-5 birefringent bound |k_(V)00^(5)| <~ 6 x 10^-34 GeV^-1 to the HAWC/Data Tables isotropic dimension-6 nonbirefringent bound |c_(I)00^(6)| < 1.24 x 10^-30 GeV^-2, the dimension-5 birefringence bound probes a suppression scale of ~1.7 x 10^33 GeV, while the dimension-6 threshold bound probes ~9 x 10^14 GeV if c ~ 1/M^2. That is about 18 to 19 orders of magnitude difference in mass scale, not 37 orders for the presently listed astrophysical nonbirefringent bound.
Important qualification, established plus your own inference. HAWC's direct paper emphasizes superluminal photon decay/splitting; the simplest Beane boson expansion above has E^2 - p^2 < 0 at leading order, i.e. a subluminal sign in the displayed convention. Therefore the HAWC ~10^-30 GeV^-2 threshold number is the strongest published Data Tables entry for the coefficient class, but applying it to the literal Beane lattice sign requires checking the sign convention and reaction channel. Sign-robust time-of-flight/dispersion bounds in the same tables are around 10^-15 GeV^-2, far weaker.
5. Does anything rescue birefringence for a lattice?
Established. Yes in the SME generally: CPT-even dimension-6 photon operators can be birefringent. Kostelecky-Mewes explicitly say k_(E)jm^(d) and k_(B)jm^(d) are CPT-even and birefringent (Phys. Rev. D 80, 015020 (2009), arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
Established. Existing birefringent dimension-6 constraints can be extremely strong. The 2026 Data Tables list, for Kostelecky and Mewes, Phys. Rev. Lett. 110, 201601 (2013), arXiv:1301.5367, https://arxiv.org/abs/1301.5367, direction-combination bounds such as |sum _2Y_jm(27 deg, 6 deg)(k_(E)jm^(6) + i k_(B)jm^(6))| <=~ 1 x 10^-31 GeV^-2 and |sum _2Y_jm(112 deg, 279 deg)(k_(E)jm^(6) + i k_(B)jm^(6))| <=~ 1 x 10^-32 GeV^-2 (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19). Friedman et al., Phys. Rev. D 102, 043008 (2020), arXiv:2003.00647, https://arxiv.org/abs/2003.00647, list individual dimension-6 birefringent photon-coefficient sensitivities around 6.6 x 10^-18 to 8.8 x 10^-18 GeV^-2.
Your own inference. This does not rescue birefringence for the specific Beane-style leading lattice artifact. A cubic direction dependence in a scalar dispersion relation is not the same as polarization-dependent propagation; it lands in c_(I)jm^(6), not k_(E)jm^(6) or k_(B)jm^(6).
Serious speculation. A more complicated discrete model could generate CPT-even birefringent k_(E)/k_(B) terms if its gauge-field discretization or microscopic update rule treats polarizations differently, violates additional symmetries, or introduces parity/electric-magnetic asymmetric structures. That would be a different model-dependent claim, not the leading artifact established by the Beane-Davoudi-Savage dispersion relation.
6. Lattice-spacing implications
Established. Beane et al.'s own bound is b^-1 ~ 10^11 GeV: "For both the fermions and the bosons, the cut off from the dispersion relation is E^max ~ 1/b. Equating this to the [GZK] cut off corresponds to a lattice spacing of b ~ 10^-12 fm, or a mass scale of b^-1 ~ 10^11 GeV. Therefore, the lattice spacing ... must be b <=~ 10^-12 fm" (Beane, Davoudi, Savage 2014, Eur. Phys. J. A 50, 148, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
Your own inference. If the best HAWC/Data Tables nonbirefringent dimension-6 coefficient bound is applicable with order-one normalization c_(I)^(6) ~ b^2, then |c_(I)00^(6)| < 1.24 x 10^-30 GeV^-2 implies b^-1 >~ 9 x 10^14 GeV; including the explicit lattice expansion factor of order 1/10 to 1/24 gives a few 10^14 GeV. That would be roughly 3 to 4 orders stronger than Beane et al.'s 10^11 GeV bound.
Your own inference. If the HAWC photon-decay/splitting bound is not applicable to the subluminal sign of the simple Beane boson dispersion, then the sign-robust dispersion bounds near 10^-15 GeV^-2 imply only b^-1 ~ 10^7 GeV for an order-one b^2 coefficient, well below Beane et al.'s cosmic-ray cutoff estimate. The sign issue is therefore not cosmetic; it decides whether modern nonbirefringent astrophysics beats the original lattice-spacing number.
Gaps and null results
- Not found. I did not find an explicit published coefficient dictionary mapping Beane-Davoudi-Savage's cubic lattice spacing
b to each SME c_(I)jm^(6) spherical coefficient, including signs and normalization.
- Not found. I did not find, in the sources checked, a statement that the Beane-Davoudi-Savage leading photon lattice artifact produces helicity-dependent propagation. The available evidence points the other way: scalar dispersion, cubic anisotropy, no helicity splitting.
- Not found. I did not find a current subluminal, nonbirefringent dimension-6 threshold bound comparable to the HAWC
10^-30 GeV^-2 superluminal photon-decay/splitting entry. The two-sided dispersion constraints I found in the Data Tables are much weaker, around 10^-15 GeV^-2.
VERDICT
Argus's inference is PARTIAL. It is right that the Beane-Davoudi-Savage cubic lattice's leading photonlike artifact is CPT-even, dimension-6, cubic-anisotropic, and nonbirefringent, so the famous CPT-odd dimension-5 vacuum-birefringence bound does not directly test that artifact. It is wrong or outdated if it implies that the relevant nonbirefringent class is constrained only in the lab or is 37 orders behind: the decisive published number is the 2026 Data Tables/HAWC entry |c_(I)00^(6)| < 12.4 x 10^-31 GeV^-2 (Albert et al. 2020, Phys. Rev. Lett. 124, 131101, arXiv:1911.08070; Kostelecky and Russell 2026, arXiv:0801.0287v19). The key qualification is sign: HAWC's mechanism is superluminal photon decay/splitting, while the simple lattice dispersion derived from Beane et al. is subluminal, so applying the strongest number to the literal lattice requires a sign-convention/reaction-channel check. If applicable it implies b^-1 ~ few x 10^14 to 10^15 GeV; if not, sign-robust dispersion bounds are only around 10^-15 GeV^-2 and Beane et al.'s b^-1 ~ 10^11 GeV remains the stronger simple estimate.
View exactly as delivered (raw text)
# CPT parity check: cubic lattice operators versus SME photon bounds
Thread date: 2026-09-10
Question: does the Beane-Davoudi-Savage cubic-lattice inference evade the strongest vacuum-birefringence bounds because it populates CPT-even, nonbirefringent dimension-6 photon-sector operators rather than CPT-odd dimension-5 birefringent ones?
## Sources read
- **Established.** Silas R. Beane, Zohreh Davoudi, Martin J. Savage, "Constraints on the Universe as a Numerical Simulation," Eur. Phys. J. A 50, 148 (2014), arXiv:1210.1847, https://arxiv.org/abs/1210.1847 and readable HTML https://ar5iv.labs.arxiv.org/html/1210.1847.
- **Established.** V. Alan Kostelecky and Matthew Mewes, "Electrodynamics with Lorentz-violating operators of arbitrary dimension," Phys. Rev. D 80, 015020 (2009), arXiv:0905.0031, https://arxiv.org/abs/0905.0031 and readable HTML https://ar5iv.labs.arxiv.org/html/0905.0031.
- **Established.** V. Alan Kostelecky and Neil Russell, "Data Tables for Lorentz and CPT Violation," Rev. Mod. Phys. 83, 11 (2011), arXiv:0801.0287; January 2026 update v19, https://arxiv.org/abs/0801.0287v19 and https://arxiv.org/html/0801.0287v19.
- **Established.** A. Albert et al. (HAWC Collaboration), "Constraints on Lorentz Invariance Violation from HAWC Observations of Gamma Rays above 100 TeV," Phys. Rev. Lett. 124, 131101 (2020), arXiv:1911.08070, https://arxiv.org/abs/1911.08070.
## 1. What the cubic lattice produces
**Established.** Beane, Davoudi, and Savage explicitly frame the regulator as a hypercubic/cubic lattice: "In lattice QCD, space-time is replaced by a finite hyper-cubic grid of points over which the fields are defined" and "The grid breaks Lorentz symmetry (and hence rotational symmetry)" (Beane, Davoudi, Savage 2014, Eur. Phys. J. A 50, 148, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
**Established.** For a boson on the lattice, their exact dispersion relation is
```text
sinh^2(b E_b / 2) - sum_{j=1,2,3} sin^2(b k_j / 2) - (b m_b / 2)^2 = 0,
```
with small-momentum expansion
```text
E_b = sqrt(|k|^2 + m_b^2) + O(b^2).
```
The same section says: "The summations are performed over the components along the lattice Cartesian axes corresponding to the x,y, and z spatial directions" and "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" (Beane, Davoudi, Savage 2014, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
**Your own inference, derived from their exact equation.** Expanding the boson relation gives
```text
E^2 - |k|^2 - m^2 + (b^2/12)(E^4 + sum_j k_j^4) + O(b^4) = 0,
```
so, on shell,
```text
E^2 = |k|^2 + m^2 - (b^2/12)[(|k|^2 + m^2)^2 + sum_j k_j^4] + O(b^4).
```
For a massless photonlike boson this is
```text
E = |k| - (b^2 |k|^3 / 24)(1 + sum_j n_j^4) + O(b^4 |k|^5), n_j = k_j / |k|.
```
**Your own inference.** The leading correction is order `b^2 |k|^3` in energy, equivalently order `b^2 |k|^4` in `E^2`; in effective-field-theory language this is a dimension-6 correction with coefficient scale `~ b^2`. Its angular dependence is the cubic invariant `sum_j n_j^4`, which decomposes into rotational `j=0` and cubic-anisotropic `j=4` pieces.
**Your own inference.** The correction is even under momentum reversal and contains no Levi-Civita/helicity sign, so it is CPT-even rather than CPT-odd. Beane et al. do not state a CPT classification in this section; I looked in their dispersion-relation and Lorentz-violation discussion and found no direct SME/CPT mapping.
**Your own inference.** The displayed boson dispersion is a scalar relation and contains no `+/-` helicity branch. It is therefore nonbirefringent at this leading order. Beane et al. do note a finite-lattice polarization artifact: "the polarizations of the massless vector fields are not exactly perpendicular to their direction of propagation for some directions of propagation with respect to the lattice axes, with longitudinal components present for non-zero lattice spacings" (Beane, Davoudi, Savage 2014, arXiv:1210.1847, https://arxiv.org/abs/1210.1847), but that sentence does not by itself establish vacuum birefringence because it does not state a helicity-dependent phase velocity.
## 2. SME photon-sector classification
**Established.** Kostelecky and Mewes define the vacuum photon coefficients as four sets: "The net result ... vacuum behavior is controlled by four sets of vacuum coefficients `c_(I)jm^(d)`, `k_(E)jm^(d)`, `k_(B)jm^(d)`, and `k_(V)jm^(d)`" (Kostelecky and Mewes 2009, Phys. Rev. D 80, 015020, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
**Established.** Their precise CPT/dimension count is: "For even d, there are `(d-1)^2` coefficients `c_(I)jm^(d)`, `(d-1)^2 - 4` coefficients `k_(E)jm^(d)`, and `(d-1)^2 - 4` coefficients `k_(B)jm^(d)`, giving a total of `3(d-1)^2 - 8` independent components associated with CPT-even Lorentz violation in vacuum propagation. For odd d, the coefficients `k_(V)jm^(d)` for CPT-odd Lorentz violation have `(d-1)^2` independent components" (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
**Established.** Their birefringence rule is explicit: "The coefficients `k_(E)jm^(d)` and `k_(B)jm^(d)` are associated with CPT-even operators that lead to birefringence, with the propagating modes being linearly polarized. The coefficients `k_(V)jm^(d)` control CPT-odd birefringence, and the corresponding eigenmodes are circularly polarized. In contrast, the CPT-even operators associated with the coefficients `c_(I)jm^(d)` are nonbirefringent" (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
**Established.** The same paper says a purely nonbirefringent vacuum model with only `c_(I)jm^(d)` has "Lorentz-violating operators only in even dimensions `d = 4, 6, 8, ...`" (Kostelecky and Mewes 2009, arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
**Established.** In minimal notation, `k_AF` is the CPT-odd photon coefficient and `k_F` is the CPT-even photon coefficient; at arbitrary dimension the vacuum classes sort into `k_(V)` for CPT-odd birefringent terms, `k_(E)` and `k_(B)` for CPT-even birefringent terms, and `c_(I)` for CPT-even nonbirefringent terms (Kostelecky and Mewes 2009, Phys. Rev. D 80, 015020, arXiv:0905.0031, https://arxiv.org/abs/0905.0031; Kostelecky and Russell 2011/2026 Data Tables, Rev. Mod. Phys. 83, 11, arXiv:0801.0287v19, https://arxiv.org/abs/0801.0287v19).
## 3. Where the cubic lattice lands
**Your own inference.** The Beane-Davoudi-Savage leading photonlike lattice artifact maps to SME `c_(I)jm^(6)`: it is even-dimension, CPT-even, dispersive, direction-dependent, and common to both helicities. The cubic angular factor `1 + sum_j n_j^4` populates an isotropic component plus cubic-anisotropic `j=4` components, not the spin-weighted polarization structures `k_(E)` or `k_(B)`.
**Established plus your own inference.** Kostelecky-Mewes establish that `c_(I)jm^(d)` is nonbirefringent; applying that rule to the scalar Beane lattice dispersion makes the leading cubic-lattice photon-sector artifact nonbirefringent (Kostelecky and Mewes 2009, arXiv:0905.0031; Beane, Davoudi, Savage 2014, arXiv:1210.1847).
**Not found.** I did not find a peer-reviewed paper that explicitly says "the Beane-Davoudi-Savage cubic lattice equals SME `c_(I)jm^(6)`" with a coefficient-by-coefficient dictionary. The mapping above is a direct classification from the published lattice dispersion and the published SME vacuum taxonomy.
## 4. Best bounds on nonbirefringent CPT-even dimension-6 photon coefficients
**Established.** The January 2026 Kostelecky-Russell Data Tables summarize the strongest isotropic sensitivities as `k_(V)00^(5)`: `10^-34 GeV^-1` and `c_(I)00^(6)`: `10^-30 GeV^-2` (Kostelecky and Russell 2011/2026, Rev. Mod. Phys. 83, 11, arXiv:0801.0287v19, Summary Table S3, https://arxiv.org/abs/0801.0287v19). The units differ, so raw exponent comparison is not physically meaningful; as suppression scales they are roughly `1e34 GeV` for dimension 5 versus `1e15 GeV` for dimension 6.
**Established.** The strongest listed `c_(I)jm^(6)` bounds in the 2026 Data Tables come from HAWC high-energy gamma-ray threshold/reaction physics, not a laboratory experiment. Table D22 lists, from A. Albert et al. (HAWC Collaboration), Phys. Rev. Lett. 124, 131101 (2020), arXiv:1911.08070, https://arxiv.org/abs/1911.08070:
| coefficient or direction combination | bound | evidence class |
|---|---:|---|
| `|c_(I)00^(6)|` | `< 12.4 x 10^-31 GeV^-2` | **Established** |
| `|sum Y_jm(103.45 deg, 276.41 deg)c_(I)jm^(6)|` | `< 3.5 x 10^-31 GeV^-2` | **Established** |
| `|sum Y_jm(83.75 deg, 286.95 deg)c_(I)jm^(6)|` | `< 4.93 x 10^-31 GeV^-2` | **Established** |
| `|sum Y_jm(67.96 deg, 83.6 deg)c_(I)jm^(6)|` | `< 20.1 x 10^-31 GeV^-2` | **Established** |
| `|sum Y_jm(53.26 deg, 304.94 deg)c_(I)jm^(6)|` | `< 50.3 x 10^-31 GeV^-2` | **Established** |
**Established.** HAWC's mechanism is explicitly threshold/reaction based: "Superluminal LIV enables the decay of photon at high energy" and "HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources" (Albert et al. 2020, Phys. Rev. Lett. 124, 131101, arXiv:1911.08070, https://arxiv.org/abs/1911.08070). The abstract reports that HAWC excludes the corresponding LIV energy scale to `2.2 x 10^31 eV` for the `n=1` case, "over 1800 times the Planck energy"; for the dimension-6/nonbirefringent `n=2` SME entries, the Data Tables conversion above gives the relevant `~10^-30 GeV^-2` sensitivities (Albert et al. 2020; Kostelecky and Russell 2026 Data Tables).
**Established.** Other `c_(I)jm^(6)` bounds listed in Table D22 are much weaker than HAWC but still astrophysical rather than laboratory: Guerrero, Campoy-Ordaz, Potting, and Gaug, Phys. Rev. D 112, 104002 (2025), arXiv:2508.02883, https://arxiv.org/abs/2508.02883, report individual `c_(I)jm^(6)` sensitivities around `10^-15 GeV^-2`; Wei, Liu, Wei, Zhang, and Wu, Universe 8, 519 (2022), arXiv:2210.03897, https://arxiv.org/abs/2210.03897, list `|c_(I)00^(6)| = 4.25(-1.63+1.60) x 10^-15 GeV^-2`; Agrawal, Singirikonda, and Desai, JCAP 05, 029 (2021), arXiv:2102.11248, https://arxiv.org/abs/2102.11248, list a direction combination at `10^(-14.2 +/- 0.1) GeV^-2`; Du et al., Astrophys. J. 906, 8 (2021), arXiv:2010.16029, https://arxiv.org/abs/2010.16029, list `|c_(I)00^(6)| < 2.4 x 10^-12 GeV^-2`; Rubtsov, Satunin, and Sibiryakov, JCAP 05, 049 (2017), arXiv:1611.10125, https://arxiv.org/abs/1611.10125, list `c_(I)00^(6) < 4 x 10^-23 GeV^-2`; Vasileiou et al., Phys. Rev. D 87, 122001 (2013), arXiv:1305.3463, https://arxiv.org/abs/1305.3463, list `c_(I)00^(6)` ranges at `10^-20 GeV^-2` scale; and Fermi-LAT/H.E.S.S./MAGIC time-of-flight bounds are listed in the `10^-19` to `10^-21 GeV^-2` range (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19).
**Established.** Laboratory/terrestrial bounds on nonminimal `c_(I)jm^(6)` are far weaker: Kostelecky, Melissinos, and Mewes, Phys. Lett. B 761, 1 (2016), arXiv:1608.02592, https://arxiv.org/abs/1608.02592, are listed with interferometry sensitivities such as `|c_(I)21^(6)| < 1.0 x 10^-2 GeV^-2` and `|c_(I)22^(6)| < 1.3 x 10^-4 GeV^-2` (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19).
**Your own inference.** Comparing the best quoted isotropic dimension-5 birefringent bound `|k_(V)00^(5)| <~ 6 x 10^-34 GeV^-1` to the HAWC/Data Tables isotropic dimension-6 nonbirefringent bound `|c_(I)00^(6)| < 1.24 x 10^-30 GeV^-2`, the dimension-5 birefringence bound probes a suppression scale of `~1.7 x 10^33 GeV`, while the dimension-6 threshold bound probes `~9 x 10^14 GeV` if `c ~ 1/M^2`. That is about 18 to 19 orders of magnitude difference in mass scale, not 37 orders for the presently listed astrophysical nonbirefringent bound.
**Important qualification, established plus your own inference.** HAWC's direct paper emphasizes superluminal photon decay/splitting; the simplest Beane boson expansion above has `E^2 - p^2 < 0` at leading order, i.e. a subluminal sign in the displayed convention. Therefore the HAWC `~10^-30 GeV^-2` threshold number is the strongest published Data Tables entry for the coefficient class, but applying it to the literal Beane lattice sign requires checking the sign convention and reaction channel. Sign-robust time-of-flight/dispersion bounds in the same tables are around `10^-15 GeV^-2`, far weaker.
## 5. Does anything rescue birefringence for a lattice?
**Established.** Yes in the SME generally: CPT-even dimension-6 photon operators can be birefringent. Kostelecky-Mewes explicitly say `k_(E)jm^(d)` and `k_(B)jm^(d)` are CPT-even and birefringent (Phys. Rev. D 80, 015020 (2009), arXiv:0905.0031, https://arxiv.org/abs/0905.0031).
**Established.** Existing birefringent dimension-6 constraints can be extremely strong. The 2026 Data Tables list, for Kostelecky and Mewes, Phys. Rev. Lett. 110, 201601 (2013), arXiv:1301.5367, https://arxiv.org/abs/1301.5367, direction-combination bounds such as `|sum _2Y_jm(27 deg, 6 deg)(k_(E)jm^(6) + i k_(B)jm^(6))| <=~ 1 x 10^-31 GeV^-2` and `|sum _2Y_jm(112 deg, 279 deg)(k_(E)jm^(6) + i k_(B)jm^(6))| <=~ 1 x 10^-32 GeV^-2` (Kostelecky and Russell 2026 Data Tables, Table D22, arXiv:0801.0287v19). Friedman et al., Phys. Rev. D 102, 043008 (2020), arXiv:2003.00647, https://arxiv.org/abs/2003.00647, list individual dimension-6 birefringent photon-coefficient sensitivities around `6.6 x 10^-18` to `8.8 x 10^-18 GeV^-2`.
**Your own inference.** This does not rescue birefringence for the specific Beane-style leading lattice artifact. A cubic direction dependence in a scalar dispersion relation is not the same as polarization-dependent propagation; it lands in `c_(I)jm^(6)`, not `k_(E)jm^(6)` or `k_(B)jm^(6)`.
**Serious speculation.** A more complicated discrete model could generate CPT-even birefringent `k_(E)`/`k_(B)` terms if its gauge-field discretization or microscopic update rule treats polarizations differently, violates additional symmetries, or introduces parity/electric-magnetic asymmetric structures. That would be a different model-dependent claim, not the leading artifact established by the Beane-Davoudi-Savage dispersion relation.
## 6. Lattice-spacing implications
**Established.** Beane et al.'s own bound is `b^-1 ~ 10^11 GeV`: "For both the fermions and the bosons, the cut off from the dispersion relation is `E^max ~ 1/b`. Equating this to the [GZK] cut off corresponds to a lattice spacing of `b ~ 10^-12 fm`, or a mass scale of `b^-1 ~ 10^11 GeV`. Therefore, the lattice spacing ... must be `b <=~ 10^-12 fm`" (Beane, Davoudi, Savage 2014, Eur. Phys. J. A 50, 148, arXiv:1210.1847, https://arxiv.org/abs/1210.1847).
**Your own inference.** If the best HAWC/Data Tables nonbirefringent dimension-6 coefficient bound is applicable with order-one normalization `c_(I)^(6) ~ b^2`, then `|c_(I)00^(6)| < 1.24 x 10^-30 GeV^-2` implies `b^-1 >~ 9 x 10^14 GeV`; including the explicit lattice expansion factor of order `1/10` to `1/24` gives a few `10^14 GeV`. That would be roughly 3 to 4 orders stronger than Beane et al.'s `10^11 GeV` bound.
**Your own inference.** If the HAWC photon-decay/splitting bound is not applicable to the subluminal sign of the simple Beane boson dispersion, then the sign-robust dispersion bounds near `10^-15 GeV^-2` imply only `b^-1 ~ 10^7 GeV` for an order-one `b^2` coefficient, well below Beane et al.'s cosmic-ray cutoff estimate. The sign issue is therefore not cosmetic; it decides whether modern nonbirefringent astrophysics beats the original lattice-spacing number.
## Gaps and null results
- **Not found.** I did not find an explicit published coefficient dictionary mapping Beane-Davoudi-Savage's cubic lattice spacing `b` to each SME `c_(I)jm^(6)` spherical coefficient, including signs and normalization.
- **Not found.** I did not find, in the sources checked, a statement that the Beane-Davoudi-Savage leading photon lattice artifact produces helicity-dependent propagation. The available evidence points the other way: scalar dispersion, cubic anisotropy, no helicity splitting.
- **Not found.** I did not find a current subluminal, nonbirefringent dimension-6 threshold bound comparable to the HAWC `10^-30 GeV^-2` superluminal photon-decay/splitting entry. The two-sided dispersion constraints I found in the Data Tables are much weaker, around `10^-15 GeV^-2`.
## VERDICT
Argus's inference is **PARTIAL**. It is right that the Beane-Davoudi-Savage cubic lattice's leading photonlike artifact is CPT-even, dimension-6, cubic-anisotropic, and nonbirefringent, so the famous CPT-odd dimension-5 vacuum-birefringence bound does not directly test that artifact. It is wrong or outdated if it implies that the relevant nonbirefringent class is constrained only in the lab or is 37 orders behind: the decisive published number is the 2026 Data Tables/HAWC entry `|c_(I)00^(6)| < 12.4 x 10^-31 GeV^-2` (Albert et al. 2020, Phys. Rev. Lett. 124, 131101, arXiv:1911.08070; Kostelecky and Russell 2026, arXiv:0801.0287v19). The key qualification is sign: HAWC's mechanism is superluminal photon decay/splitting, while the simple lattice dispersion derived from Beane et al. is subluminal, so applying the strongest number to the literal lattice requires a sign-convention/reaction-channel check. If applicable it implies `b^-1 ~ few x 10^14` to `10^15 GeV`; if not, sign-robust dispersion bounds are only around `10^-15 GeV^-2` and Beane et al.'s `b^-1 ~ 10^11 GeV` remains the stronger simple estimate.