Bell Algebra without Bell Violation: CHSH Operators for Ising Anyons in the Kitaev Honeycomb Model

Abstract

We construct four spatially separated Majorana zero modes (MZMs) per party on a finite Kitaev honeycomb lattice in its non-Abelian phase, bound to well-separated Ising vortices, and use them to realize a local, noncommuting CHSH measurement algebra. Two static fusion-parity observables per party, Z_A=iγ_{A,1}γ_{A,2} and X_A=iγ_{A,2}γ_{A,3} (and analogously for Bob), commute exactly between parties (P₂=0) and fail to commute within a party (P₁=‖[Z_A,X_A]‖=2.000000); on the maximally entangled fusion-parity Bell state the CHSH expectation reaches the algebraic Tsirelson value S=2√2=2.828427. This value survives fermion-parity superselection exactly, because the Bell state is a definite-parity eigenstate with unit weight in its own sector — a statement specific to the present four-MZM, measurement-only construction, not to braid-then-fuse protocols, which are a different experiment. A commuting-only control setting returns K=√2 ≤ 2, and the finite-size Majorana footprint leaking across the bipartition decays monotonically with system size (isotropic scaling 17.0% → 10.2% → 2.3%). We stress the critical scope-limiting fact, following Howard and Vala [Phys. Rev. A 85, 022304 (2012)]: a physical Bell-inequality violation using only topologically protected Ising-anyon operations does not follow from this construction, since Ising braiding generates only the Clifford group; that step requires a non-Clifford ("magic") resource, the subject of a companion paper. What we report here is the local realization of the noncommuting Bell-operator algebra on a genuine microscopic lattice substrate, not a claim of Bell nonlocality. A loop-threading relative phase δ reproduces the analytic corollary S(δ)=2√2·cos²(δ/2) to a residual of 4.4×10⁻¹⁶, and an independent gauge-invariant flux-holonomy diagnostic confirms no adiabatic channel mixing at the quoted scale (‖[W,Z_A]‖_F=0.0132), with the sole open discrepancy being a 0.020-rad gap between the accumulated Berry phase at one full flux quantum and the target π, traced to genuine near-degeneracies along the loop rather than to a construction error. About this series: This record is part of a series of related works from my independent research on Fibonacci anyons, with Ising anyons as their natural counterpart. I started in April 2026, and it has been a long and insightful journey in which I learned a lot; the work uses different methods and stays within verifiable, nonspeculative physics. The common thread of the series is a split: Ising anyons are limited to Clifford operations, while Fibonacci anyons are computationally universal, and across the series I map what standard witnesses of nonclassicality can and cannot certify on such systems. I consider Fibonacci anyons a serious candidate for topological quantum computing, given their universality and their topological protection against local noise. A hybrid approach with Ising is conceivable, but problems such as instability and certification would have to be solved first, and each needs research of its own. This paper asks whether a real microscopic lattice can carry the Bell structure, and constructs four Majorana modes per party on a Kitaev honeycomb lattice in its non-Abelian phase, realizing the full CHSH measurement algebra locally, with observables that commute between the parties and fail to commute within a party; the algebra holds, reaching the algebraic Tsirelson value and surviving the fermion-parity superselection, but a physical Bell violation does not follow, because topologically protected Ising operations generate only Clifford gates and the missing step is a non-Clifford resource, which is the subject of a companion paper; within the series it is the substrate test, where algebra on real material is not yet a violation. Use of AI tools: In the research, processing, and writing of this paper and its results I worked together with generative AI tools, in practice a system of multiple coordinated AI instances that I set up and orchestrate (large language models, mainly Claude, by Anthropic, inside Claude Code). At their current context sizes I found it far more effective to work with several specialized instances, each with its own role and its own harness of rules and parameters that I designed and refined through feedback, than to load a single instance with all of the material; for my workflow that would have been inefficient, though this depends on the individual implementation. I lead this collaboration: I choose the research directions, set the goals, and make the final decisions in open exchange with the AI, learning actively as the work proceeds. The AI carries out the drafting, including the mathematical and technical parts, the numerical computation, and the literature search, under my direction. The AI works autonomously only task by task, within the structure I develop through feedback: it completes a task, and at open questions that need me it stops until the point is settled before the next step. Along the way I witness and take many of the decisions that shape the path, and it is common for me to spot things that need improvement. The work spans many separate runs, and a single simulation or build task alone can take up to an hour, so it could not happen all together in one autonomous run; and had I let the AI do all of it together alone, even if it is possible, it would no longer be my work but the AI's. I run multiple verifications at the different stages of the work and one before release, including cross-checks with an unrelated AI model from a different company, and all references are checked against the original sources. In the end what matters are human eyes, a principle that is itself written into the parameters of my system: I reach out to experts after publishing for review and feedback, so I learn what is solid and what must be corrected or falsified. My scripts for reproduction and review are released with this record. These tools are not authors; I am the author, and I take full responsibility for all scientific content and decisions leading to these results and their publication.---------------- Version notes (v1.0 → v1.1, following a comprehensive internal review of the full series): • No previously reported value changed. What is new are descriptive quantities: the lattice sizes and separations of the finite-size scan, and the accuracy to which one construction identity holds. • Correction (printed equation): Eq. (4) corrected — the printed form did not reflect the CHSH functional actually evaluated with Bob's rotated settings. All reported values are unchanged. The equation was written in terms of Bob's unrotated parities while his rotated settings were defined inline immediately below it; the functional is now written in terms of those settings, and their definition is a numbered equation that it references. • Correction (terminology in a statement about prior work): the resource that Brennen et al. require in addition to Ising braiding is now called a non-stabilizer resource rather than a non-Clifford resource, in the three places where the requirement is stated. This is what their argument actually requires. • The finite-size quantity is renamed and its scan is made explicit: what was reported as the "point-value cross-party commutator footprint" is the mean cross-party support weight of Alice's modes on Bob's region. The scan is now given with its lattice sizes (24x12, 32x16 and 44x22) and inter-cluster separations (d = 18, 24 and 33), and a sentence states that both finite-size quantities vanish as L grows and that neither enters a Bell value. • Scope of the Howard-Vala re-scoping — sharpened: the construction is now described as using only Clifford combinations of static fusion-parity observables, with Alice measuring the parities directly and Bob's CHSH settings being the 45-degree-rotated combinations; and the Howard-Vala condition is stated to require additionally that the measurements be stabilizer operations. • Scope limit added: B and B' are treated at the level of the ideal logical Majorana algebra; their realization as individual topologically protected measurements is not constructed here. • Gauge invariance — now discussed rather than assumed: the parity observables are bilinears in the matter Majoranas at spatially separated vortex cores and are therefore not by themselves invariant under the Z2 gauge generators; the gauge-invariant observable carries a Wilson line between the two cores. Within a fixed flux sector — the setting used throughout — that Wilson line is a number +-1 fixed by the enclosed flux, so the dressed and undressed operators act identically and every value reported here is unchanged. A reference is added for this (Petrova et al.). • Terminology disambiguated: the basis invariance of W is the non-Abelian Wilczek-Zee gauge freedom of the degenerate subspace, not the Z2 gauge redundancy generated by the D_i. A reference is added (Wilczek and Zee). • Orthonormality restated: the lattice modes enter the algebra already orthonormal — they are real Schur vectors of the real antisymmetric A-hat, so G^T G = I is an identity of the construction rather than a numerical finding, holding to the accuracy of the Schur decomposition (2.3e-14) before any orthonormalization step. • The delta-corollary is qualified where it is first stated: the qualifier that the body already carried — that S(delta) is a consistency check on the algebra rather than a free-standing result — is now also carried in the summary passage, together with the lattice size on which the corollary is reproduced (28x16). A qualifier that only appears deep in the body does not reach a reader of the summary. • Newly deposited: the numerical driver for the four locked targets (p5b_m0_numeric.py) and the finite-size scan it produces (p5b_m0_finite_size.json), so that the reported targets and the scan can be recomputed from deposited data. • Licensing — now stated in the paper, and