Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The MERA Scaling Spectrum at χ = 6: A Tensor-Network Realization of the ⟨2, 3⟩ Modular Ratio Set

Daland Montgomery

Published May 16, 2026• 0 citations
Publisher

Abstract

Description We construct the multiscale entanglement renormalization ansatz (MERA) at bond dimension χ = 6 with the natural binary × ternary coarse-graining of the framework's 2 × 3 cell grid and identify its scaling structure at the conserved-current fixed point. The scaling-weight operator W on End(ℂ⁶) — the MERA descending operator at the fixed point — has finite point spectrum {1, 1/2, 1/3, 1/6} on four primary sectors of dimensions (1, 3, 8, 24) summing to 36 = χ². All non-trivial sectors carry primary scaling dimension Δ = 1 — the universal conserved-current dimension forced by the Sugawara construction (Knizhnik–Zamolodchikov 1984; Witten 1984) — and the four scaling weights come from the four branching ratios b ∈ {1, 2, 3, 6} of the cell grid. The corresponding scaling superoperator 𝒮(X) := W X W⁻¹ has finite point spectrum equal to the nine-element ratio set {1, 1/2, 1/3, 1/6, 2/3, 3/2, 2, 3, 6}. The paper establishes the operator-algebraic bridge between the tensor-network and modular-theory readings of the framework's finite-dimensional core. The MERA scaling-weight operator W and the framework's per-cell Bisognano–Wichmann density ρ_ω on A_F = ℂ ⊕ M₂(ℂ) ⊕ M₃(ℂ) (Montgomery 2026a) are the same diagonal positive operator in the Wedderburn–Artin basis; the scaling superoperator 𝒮 and the Tomita–Takesaki modular operator Δ_ω are the same multiplicative conjugation; the nine-element spectra coincide exactly. The Sd-invariant is ⟨2, 3⟩, Pontryagin-dual to the compact 2-torus 𝕋². Three further results sharpen the construction. (i) The fixed-point projection is stated for the physical scaling superoperator 𝒮_phys, the restriction of 𝒮 to the descending sector ℰ_≤ = ℰ_< ⊕ ℰ_= — the rigorous form of the MERA causal-cone / local-operator-support condition — with convergence ‖𝒮_phys^n(X) − P_singlet(X)‖_F ≤ (2/3)^n ‖X − P_singlet(X)‖_F at the slowest contracting eigenvalue 2/3. (ii) A falsifiable prediction: the framework's affine current sector at level 1 is û(1) ⊕ ŝu(2)₁ ⊕ ŝu(3)₁ — the Standard Model gauge group's affine Kac–Moody algebras — with total Sugawara central charge c_total = 1 + 1 + 2 = 4, testable by variational MERA optimization on A_F at χ = 6. (iii) A three-witness convergence: the same four-element spectrum {1, 1/2, 1/3, 1/6} arises from Wedderburn-sector traces (algebraic / NCG), MERA scaling at Δ = 1 (tensor-network), and Sugawara WZW current scaling (CFT) — three independent operator-theoretic calculations on the same algebraic data. Companion materials in the repository: a self-contained Python script verifying seven receipts (Wedderburn closure, both spectra, geometric contraction at rate ≤ 2/3, Sd-invariant generators, Pontryagin pairing, the W = ρ_ω / 𝒮 = Δ_ω bridge, Weyl equidistribution); kernel-checked Lean 4 and Agda proofs of the discrete content (closed by decide and refl with no sorry, no axioms beyond standard, no mathlib); a three.js interactive demo of the MERA tree, the four primary sectors, and the nine-element ratio lattice. Built on the Crystal Topos foundation: load the four foundational papers BEFORE this one, one at a time in separate uploads, to give an LLM full context — single bundled uploads compress each paper into a thin summary and lose the algebra. Recommended order: (1) Discrete Core https://doi.org/10.5281/zenodo.20232800 (2) MERA Scaling Spectrum https://doi.org/10.5281/zenodo.20236511 (3) Type II∞ Core https://doi.org/10.5281/zenodo.20257603 (4) Triangular Duality https://doi.org/10.5281/zenodo.20257953 After all four are ingested, upload this paper as the fifth turn. Keywords MERA, multiscale entanglement renormalization ansatz, tensor networks, scaling superoperator, conserved-current fixed point, Wess–Zumino–Witten model, Sugawara construction, affine Kac–Moody algebra, level-1 representation, central charge, c=4, Standard Model gauge group, û(1), ŝu(2)₁, ŝu(3)₁, modular theory, Tomita–Takesaki, Bisognano–Wichmann, Connes Sd-invariant, almost-periodic state, Pontryagin duality, 2-torus, Weyl equidistribution, Mihăilescu, Catalan conjecture, A_F, Wedderburn–Artin, three-witness convergence, χ=6, bond dimension, ratio set, ⟨2,3⟩, Crystal Topos, WACA, discrete spectral-triple holography, kernel-checked Lean, Agda Copyright © 2026 Daland Montgomery. Creative Commons Attribution Non Commercial No Derivatives 4.0 International For a private, proprietary license, contact: quidbit@icloud.com Software Implementation: The formulas and constants derived in this work are implemented in the CrystalAgent engine, available under the AGPL-3.0 license at: https://github.com/CrystalToe/CrystalAgent.

DOI
10.5281/zenodo.20236511

Related papers

Approximation of the Riemann Zeta Function via MERA Tensor Networks: A Study of Convergence and Error Bounds Based on Mellin Operators

Juan Arroyo

This preprint presents a rigorous, self-contained mathematical framework for numerically approximating the Riemann Zeta function using Multiscale Entanglement Renormalization Ansatz (MERA) tensor networks. By formally establishing a unitary equivalence between the MERA coarse-graining operator and a discretized Mellin integral (Lemma 3.1), we analytically derive the theoretical infrared fixed point (Theorem 3.1) to be γtarget=log⁡(2π)2π2≈0.0931γtarget=2π2log(2π)≈0.0931. This derivation leverages the polygamma function evaluated at the center of the critical strip. A numerical implementation with bond dimension χ=2χ=2 achieves this target with an absolute error of 1.3×10−51.3×10−5 and a relative error of +0.0139%+0.0139%. We also provide strict, exponentially converging error bounds based on the singular value decomposition of Hankel/Mellin-type operators (Theorem 3.2). Notably, this mathematically derived fixed point provides the fundamental scale parameter for the Renormalization Group Scale (RG-S) within the R5+ modified gravity model. Keywords: Riemann Zeta Function, MERA, Tensor Networks, Mellin Transform, Renormalization Group, Functional Equation, Fixed Point, Numerical Analysis, R5+ Gravity.

•Aug 11, 2026•0 citations

Approximation of the Riemann Zeta Function via MERA Tensor Networks: A Study of Convergence and Error Bounds Based on Mellin Operators

Juan Arroyo

This preprint presents a rigorous, self-contained mathematical framework for numerically approximating the Riemann Zeta function using Multiscale Entanglement Renormalization Ansatz (MERA) tensor networks. By formally establishing a unitary equivalence between the MERA coarse-graining operator and a discretized Mellin integral (Lemma 3.1), we analytically derive the theoretical infrared fixed point (Theorem 3.1) to be γtarget=log⁡(2π)2π2≈0.0931γtarget=2π2log(2π)≈0.0931. This derivation leverages the polygamma function evaluated at the center of the critical strip. A numerical implementation with bond dimension χ=2χ=2 achieves this target with an absolute error of 1.3×10−51.3×10−5 and a relative error of +0.0139%+0.0139%. We also provide strict, exponentially converging error bounds based on the singular value decomposition of Hankel/Mellin-type operators (Theorem 3.2). Notably, this mathematically derived fixed point provides the fundamental scale parameter for the Renormalization Group Scale (RG-S) within the R5+ modified gravity model. Keywords: Riemann Zeta Function, MERA, Tensor Networks, Mellin Transform, Renormalization Group, Functional Equation, Fixed Point, Numerical Analysis, R5+ Gravity.

•Aug 11, 2026•0 citations

The Crystal Toolbox: A Complete Dictionary of Crystal, QFT, and MERA Operations

Daland Montgomery

Abstract: The finite algebra A_F = ℂ ⊕ M₂(ℂ) ⊕ M₃(ℂ) at KMS inverse temperature β = 2π provides a complete, finite set of operations for computing physical constants. We classify every operation into five levels and establish a bijective dictionary between spectral operations and standard quantum field theory. The ascending superoperator S on End(ℂ⁶) is the MERA transfer matrix. Its eigenvalues {1, 1/2, 1/3, 1/6} are the anomalous dimensions. Its degeneracies {1, 3, 8, 24} are the channel multiplicities. The first-order channel count Σd = 36 gives tree-level results (33/38 within 1% of measurement). The Schur square Σd² = 650 gives all corrections: the universal MERA retention factor 35/36 = 1 − 1/χ² is the wavefunction renormalisation Z; the Plancherel resolvent Σd²ln(1/(1−λ)) = 137.205 is the spectral zeta function; the Ward anomalies (1−λ_k) = {0, 1/2, 2/3, 5/6} are the Ward–Takahashi identities that ensure one-loop renormalisability (van Nuland & van Suijlekom, JHEP 2022). The spectral truncation error O(1/χ) at finite bond dimension χ = 6 corresponds to UV cutoff artefacts in the Connes–van Suijlekom operator system framework (CMP 2020). Standard model running provides boundary-condition transport from the conformal point (D = 0, ~245 GeV) to the measurement scale. The type-I seesaw is built into the tower formula m_ν = v/2⁴², with the ~10% tree-to-loop gap confirmed as universal finite-χ truncation across five independent hierarchical branching domains. The dictionary contains 19 entries. Each is a bijection: every crystal operation maps to exactly one QFT operation. The five outliers outside 1% each resolve through a specific, standard QFT mechanism — wavefunction renormalisation (α_s/α_w), UV cutoff (horizon entropy), Yukawa running (m_b), and lattice discretisation (neutrinos) — applied to the spectral data of A_F. Three open calculations remain: the van Nuland–van Suijlekom one-loop spectral action evaluated for the specific triple (A_F, H, D); the derivation of Z = 35/36 from first principles for each observable; and the full ΛCDM integration with crystal inputs. The toolbox is finite: 36 channels give tree-level, 650 channels give all corrections, 1296 paths give mixing. No infinite tower of corrections, no ad hoc parameters. The crystal does not replace quantum field theory. It is quantum field theory, written in one algebra. Keywords: spectral action, noncommutative geometry, ascending superoperator, MERA, tensor network, Schur square, wavefunction renormalisation, Ward–Takahashi identities, anomalous dimensions, spectral zeta function, one-loop renormalisability, spectral truncation, operator systems, KMS state, Bisognano–Wichmann, Plancherel measure, finite bond dimension, tree-level closure, seesaw mechanism, crystal toolbox, WACA

•Mar 19, 2026•0 citations

Entanglement compression in scale space: From the multiscale entanglement renormalization ansatz to matrix product operators

Karel Van Acoleyen, Andrew Hallam, Matthias Bal +3 more

The multiscale entanglement renormalization ansatz (MERA) provides a constructive algorithm for realizing wave functions that are inherently scale invariant. Unlike conformally invariant partition functions, however, the finite bond dimension $\ensuremath{\chi}$ of the MERA provides a cutoff in the fields that can be realized. In this paper, we demonstrate that this cutoff is equivalent to the one obtained when approximating a thermal state of a critical Hamiltonian with a matrix product operator (MPO) of finite bond dimension $\ensuremath{\chi}$. This is achieved by constructing an explicit mapping between the isometries of the MERA and the local tensors of the MPO. In terms of energy scales, our results show that a finite bond dimension MERA is equivalent to introducing both an infrared and an ultraviolet scale, characterizing relevant and irrelevant perturbations on the underlying conformal field theory.

arXiv:1912.10572•Physical review. B./Physical review. B•Oct 15, 2020•3 citations

Power-Law Holographic Forces and Causal Transport in MERA Quantum Circuits

Jon Cackler

CORRECTION NOTICE (v5, 2026-09-06). The central result of v1-v4 is WITHDRAWN. This record has not been peer reviewed. Issued 2026-09-06 by the author (v5 of this record). This notice was added to the archived record; the manuscript itself is unaltered. This record has not been peer reviewed. What is wrong The force exponent reported here as alpha = 3.92 +/- 0.16 is an artifact of the fitting procedure, not a property of the data. The fit was unweighted least squares in linear space applied to a steeply falling profile, which places 95.8 percent of its objective on the single point at the shortest retained separation. Re-fitting the same data with nothing discarded, changing only the error model, gives 2.82 (relative-error weighting) or 2.89 (log-space regression). Across fifteen combinations of standard estimator and fit window the exponent spans 2.53 to 3.92, and the value published here is the largest of them. What the data show instead, and how firmly The retraction is certain. The replacement is not. On the open 24-site chain the same data give alpha = 2.6 +/- 0.2, where the interval is the spread of a plateau in the local logarithmic slope over six consecutive separations, not a statistical error (the bootstrap error on the published fit to the same data is 0.63); the plateau value was reproduced at 2.5921 by a separately written implementation run on the same recorded state (SIM-2026-073). This is a single geometry at a single size. On a 24-site periodic ring, analysed with chord distance as pre-registered, two probes that are both energy-density-like (Delta = 1) give 1.83 and 2.37 against a CFT two-point value of 2 (SIM-2026-074): one size, not converged, disagreeing with each other by 29 percent, and not attributable to a single operator channel. No corrected exponent is established. The exchange channel is an open question. Treat 2.6 +/- 0.2 as the provisional open-chain value, not as a result. Also withdrawn or voided alpha = 3.92 +/- 0.16 (24 sites) and alpha = 3.97 +/- 0.15 (16 sites, exact) from stress-tensor exchange. Was: headline result and its first corroboration. Now: withdrawn. Why: the same estimator artifact: both move in lockstep under it (3.92 to 2.56, 3.97 to 2.75 on dropping the contact point), because they share one observable, one fit window and one fitting routine. scalar exchange (alpha = 2) excluded. Was: stated as excluded at about 12 sigma. Now: withdrawn. Why: the 12 sigma used the covariance error 0.16; at the bootstrap error 0.63 computed before release the exclusion is about 2.2 sigma, and under any window starting at separation 3 the point estimate itself falls to about 2.56. alpha = 4 is universal; alpha = 4.00 at Ising, tricritical Ising and Potts. Was: confirmed at three classes. Now: withdrawn; the operator towers at all three classes stand. Why: the exponent was measured only at Ising and inherited by the other two; the towers are spectra, not forces. third independent evaluation alpha = 4.07 (chi = 3 MERA), and with it the statement that this was the first holographic force measurement from an optimized tensor network state. Was: supporting evidence and a priority claim. Now: both withdrawn. Why: the profile changes sign inside its own fit window and cannot be a fixed-sign power law; the priority claim rests on that evaluation alone and was never gated against prior work. the discriminating test: stress tensor (Z_j Z_j+1) versus spin primary (X_j). Was: registered killer test. Now: void. Why: at the critical point these are one operator, related by exact Kramers-Wannier self-duality. Measured on the periodic ring: identical expectation values to 1e-15 and identical force profiles to 1e-13 at 14 sites, 2.4e-9 at 24 sites (SIM-2026-074). Open boundaries break the duality, which is why the open-chain run of this test (SIM-2026-073) showed a spurious difference between the two arms. The test compared the energy density to itself. Neither label was correct: Z_j Z_j+1 is the energy density, not the stress tensor; X_j is its dual, not the spin primary. transport time T = (pi/2) l_AdS to 1.2 percent as a geometric prediction. Was: confirmed. Now: measurement stands (T = 2.239 at N = 16); geometric reading not established. Why: l_AdS cancels in the cited derivation, which also gives pi rather than the pi/2 the manuscript compares against; the chi = 3 estimate of l_AdS gives a 57 percent miss; the transport driver never applies a MERA. The true spin primary Z_j was never in the design, and the direct replacement test is blocked: Z_j is odd under the Z2 symmetry, so the energy-difference observable used throughout this work subtracts across symmetry sectors and is undefined for it. A replacement test using two Z2-even probes has been pre-registered (PREREGISTRATION-T.md, committed before its driver) and run; its result is deliberately not part of this correction and will appear in a later version. Not previously reported The published fit has a reduced chi-squared of 6.008 on 9 degrees of freedom, and a bootstrap standard deviation of 0.63 rather than the quoted 0.16. Both were computed before release and omitted. What is unaffected Of the claims in this record: the two-phase optimization method, the chi = 3 criticality threshold, the exact-diagonalization gates, the measured OPE coefficient C_sigma-sigma-epsilon = 0.50 +/- 0.03 against the exact 1/2, and the conformal towers at all three classes. Not part of this record, mentioned only because they are sometimes read together with it: the later C_TTT stress-tensor self-coupling measurements (SIM-2026-067/068/069, dated after this record) and the matched off-critical ring control (SIM-2026-074) are unaffected by this correction and make no claim about the force exponent. How the error was found and what this record contains The error was found by the author's own internal adversarial audit, not by a reviewer; the journal submission it was found after is closed, and no review of this work has taken place. Files in this record: the full re-analysis (SIM-2026-072.json, driver, figure), the pre-registered chain test (SIM-2026-073.json) and ring test (SIM-2026-074.json) with both pre-registrations as committed before the runs, the runnable drivers with their shared module and the negative test suite for their validation gates, and the error register (ERRATA.md). Concept DOI 10.5281/zenodo.19660302 always resolves to the current version. The abstract of the original manuscript follows unchanged, for the record. Every claim in it about the force exponent should be read against the correction above. We report results from classical simulation of critical spin chains that establish quantitative physical content in the MERA-AdS correspondence. First, the interaction energy between boundary-primary excitations follows a power law F(dr) ~ dr^(-alpha) with three consistent evaluations yielding alpha = 3.92 +/- 0.16 (24q), 3.97 +/- 0.15 (16q), and 4.07 +/- 0.18 (16q chi=3 MERA), all consistent with stress-tensor exchange (alpha = 2*Delta_T = 4) and inconsistent with scalar exchange (alpha = 2); the exponent is explained by the Ising CFT fusion rule eps x eps = [1]. Because Delta_T = 2 is fixed by conformal symmetry in every unitary 2D CFT, alpha = 4 is universal: we confirm this at two further fixed points - the tricritical Ising model (c = 7/10), where a scale-invariant MERA recovers the full M(4,5) primary spectrum (Delta_sigma = 0.092, Delta_eps = 0.204, Delta_T = 1.999), and the three-state Potts model (c = 4/5), where a Z3-symmetric MERA at chi = 12 recovers Delta_sigma = 0.134, Delta_eps = 0.805, Delta_T = 2.001 - giving alpha = 4.00 at all three universality classes. Second, the holographic transport time matches the AdS_3 null geodesic prediction T = (pi/2) l_AdS to 1.2% at N = 16 sites on a time-grid-converged measurement, with the residual shrinking from 2.8% at N = 12. Third, we introduce a two-phase optimization method (ALM + Z2 projection) for Z2-symmetric MERA ground states at 24 qubits. Fourth, extending the MERA bond dimension from chi=2 to chi=3 overcomes the expressibility barrier, achieving 0.15% energy error at 16 qubits and enabling the first holographic force measurement from an optimized tensor network state. The OPE coefficient is also measured directly from the fixed-point tensors: C_sigma-sigma-eps = 0.50 +/- 0.03 against the exact Ising value 1/2. A hardware path is designed for Quantinuum Helios (57 CNOTs, 0.57% logical error). Version 4 changes (relative to v3): Transport: N=16 residual restated as the time-grid-converged value of -1.21% (dt = 1e-3), with -2.78% at N=12 showing the finite-size trend; the 0.72% figure in v3 was an artifact of the dt = 0.05 integration grid and is now labeled as such. The N=24 transport point carries a new reliability caveat (double-lobe trace structure). The l_AdS comparison table is recomputed against the converged transport time (best fit 1.425, deviation 1.21%). Effective-mass derivation rewritten in a consistent AdS3 convention m^2 l^2 = Delta(Delta-2); m_eff = 0.693 = 1/l_AdS follows from Delta = 1 exactly. NEW: universality section extended to a third fixed point - a Z3-symmetric MERA at chi = 12 recovers the three-state Potts (c = 4/5) tower, giving alpha = 4.00 at all three universality classes (c = 1/2, 7/10, 4/5). NEW: the OPE coefficient C_sigma-sigma-eps = 0.50 +/- 0.03 measured directly from the fixed-point tensors, with C_sigma-sigma-sigma = 0 at machine precision as the Z2 structural check. Supplemental table row SIM-007 annotated (record lost; script preserved); data availability expanded; Zenodo DOI added to the front matter.

•Apr 20, 2026•1 citation

Multiscale Entanglement Renormalisation Ansatz

Markus Hauru

This thesis reviews the multiscale entanglement renormalisation ansatz or MERA, a numerical tool for the study of quantum many-body systems and a discrete realisation of the AdS/CFT duality. The thesis covers an introduction to the necessary background concepts of entanglement, entanglement entropy and tensor network states, the structure and main features of MERA and its applications in condensed matter theory and holography. Also covered are details on the algorithmic implementation of MERA and some of its generalisations and extensions. MERA belongs to a class of variational ansatze for quantum many-body states known as tensor network states. It is especially well-suited for the study of scale invariant critical points. MERA is based on a real-space renormalisation group procedure called entanglement renormalisation, designed to systematically handle entanglement at different length scales along the coarse-graining flow. Entanglement renormalisation has be used for example to efficiently describe Kitaev states of the toric code, the prime example of topological order, and numerically study the ground state of the highly frustrated spin-1/2 Heisenberg model on a kagome lattice and various other one- and two-dimensional lattice models. The geometric and causal structure of MERA, which underlies its effectiveness as a numerical tool, also makes it a discrete version of the AdS/CFT duality. This duality describes a conformal field theory by a gravity theory in a higher dimensional space, and vice versa. The duality is manifest in the scaling of entanglement entropy in MERA, which is governed by a law highly analogous to the Ryu-Takayanagi formula for holographic entanglement entropy, in the connection between thermal states and a black-hole-like MERA and in the connection between correlation functions and holographic geodesics in a scale invariant MERA. The aim of this thesis is to lead the reader to an understanding of what MERA is, how it works and how it can be used. MERA's core features and uses are presented in a comprehensive and explicit way, and a broad view of possible applications and further directions is given. Plenty of references are also offered to direct the reader to further research on how MERA may relate to his/her interests.

Työväentutkimus Vuosikirja•Nov 19, 2013•6 citations
Description We construct the multiscale entanglement renormalization ansatz (MERA) at bond dimension χ = 6 with the natural binary × ternary coarse-graining of the framework's 2 × 3 cell grid and identify its scaling structure at the conserved-current fixed point. The scaling-weight operator W on End(ℂ⁶) — the MERA descending operator at the fixed point — has finite point spectrum {1, 1/2, 1/3, 1/6} on four primary sectors of dimensions (1, 3, 8, 24) summing to 36 = χ². All non-trivial sectors carry primary scaling dimension Δ = 1 — the universal conserved-current dimension forced by the Sugawara construction (Knizhnik–Zamolodchikov 1984; Witten 1984) — and the four scaling weights come from the four branching ratios b ∈ {1, 2, 3, 6} of the cell grid. The corresponding scaling superoperator 𝒮(X) := W X W⁻¹ has finite point spectrum equal to the nine-element ratio set {1, 1/2, 1/3, 1/6, 2/3, 3/2, 2, 3, 6}. The paper establishes the operator-algebraic bridge between the tensor-network and modular-theory readings of the framework's finite-dimensional core. The MERA scaling-weight operator W and the framework's per-cell Bisognano–Wichmann density ρ_ω on A_F = ℂ ⊕ M₂(ℂ) ⊕ M₃(ℂ) (Montgomery 2026a) are the same diagonal positive operator in the Wedderburn–Artin basis; the scaling superoperator 𝒮 and the Tomita–Takesaki modular operator Δ_ω are the same multiplicative conjugation; the nine-element spectra coincide exactly. The Sd-invariant is ⟨2, 3⟩, Pontryagin-dual to the compact 2-torus 𝕋². Three further results sharpen the construction. (i) The fixed-point projection is stated for the physical scaling superoperator 𝒮_phys, the restriction of 𝒮 to the descending sector ℰ_≤ = ℰ_< ⊕ ℰ_= — the rigorous form of the MERA causal-cone / local-operator-support condition — with convergence ‖𝒮_phys^n(X) − P_singlet(X)‖_F ≤ (2/3)^n ‖X − P_singlet(X)‖_F at the slowest contracting eigenvalue 2/3. (ii) A falsifiable prediction: the framework's affine current sector at level 1 is û(1) ⊕ ŝu(2)₁ ⊕ ŝu(3)₁ — the Standard Model gauge group's affine Kac–Moody algebras — with total Sugawara central charge c_total = 1 + 1 + 2 = 4, testable by variational MERA optimization on A_F at χ = 6. (iii) A three-witness convergence: the same four-element spectrum {1, 1/2, 1/3, 1/6} arises from Wedderburn-sector traces (algebraic / NCG), MERA scaling at Δ = 1 (tensor-network), and Sugawara WZW current scaling (CFT) — three independent operator-theoretic calculations on the same algebraic data. Companion materials in the repository: a self-contained Python script verifying seven receipts (Wedderburn closure, both spectra, geometric contraction at rate ≤ 2/3, Sd-invariant generators, Pontryagin pairing, the W = ρ_ω / 𝒮 = Δ_ω bridge, Weyl equidistribution); kernel-checked Lean 4 and Agda proofs of the discrete content (closed by decide and refl with no sorry, no axioms beyond standard, no mathlib); a three.js interactive demo of the MERA tree, the four primary sectors, and the nine-element ratio lattice. Built on the Crystal Topos foundation: load the four foundational papers BEFORE this one, one at a time in separate uploads, to give an LLM full context — single bundled uploads compress each paper into a thin summary and lose the algebra. Recommended order: (1) Discrete Core https://doi.org/10.5281/zenodo.20232800 (2) MERA Scaling Spectrum https://doi.org/10.5281/zenodo.20236511 (3) Type II∞ Core https://doi.org/10.5281/zenodo.20257603 (4) Triangular Duality https://doi.org/10.5281/zenodo.20257953 After all four are ingested, upload this paper as the fifth turn. Keywords MERA, multiscale entanglement renormalization ansatz, tensor networks, scaling superoperator, conserved-current fixed point, Wess–Zumino–Witten model, Sugawara construction, affine Kac–Moody algebra, level-1 representation, central charge, c=4, Standard Model gauge group, û(1), ŝu(2)₁, ŝu(3)₁, modular theory, Tomita–Takesaki, Bisognano–Wichmann, Connes Sd-invariant, almost-periodic state, Pontryagin duality, 2-torus, Weyl equidistribution, Mihăilescu, Catalan conjecture, A_F, Wedderburn–Artin, three-witness convergence, χ=6, bond dimension, ratio set, ⟨2,3⟩, Crystal Topos, WACA, discrete spectral-triple holography, kernel-checked Lean, Agda Copyright © 2026 Daland Montgomery. Creative Commons Attribution Non Commercial No Derivatives 4.0 International For a private, proprietary license, contact: quidbit@icloud.com Software Implementation: The formulas and constants derived in this work are implemented in the CrystalAgent engine, available under the AGPL-3.0 license at: https://github.com/CrystalToe/CrystalAgent.
This preprint presents a rigorous, self-contained mathematical framework for numerically approximating the Riemann Zeta function using Multiscale Entanglement Renormalization Ansatz (MERA) tensor networks. By formally establishing a unitary equivalence between the MERA coarse-graining operator and a discretized Mellin integral (Lemma 3.1), we analytically derive the theoretical infrared fixed point (Theorem 3.1) to be γtarget=log⁡(2π)2π2≈0.0931γtarget=2π2log(2π)≈0.0931. This derivation leverages the polygamma function evaluated at the center of the critical strip. A numerical implementation with bond dimension χ=2χ=2 achieves this target with an absolute error of 1.3×10−51.3×10−5 and a relative error of +0.0139%+0.0139%. We also provide strict, exponentially converging error bounds based on the singular value decomposition of Hankel/Mellin-type operators (Theorem 3.2). Notably, this mathematically derived fixed point provides the fundamental scale parameter for the Renormalization Group Scale (RG-S) within the R5+ modified gravity model. Keywords: Riemann Zeta Function, MERA, Tensor Networks, Mellin Transform, Renormalization Group, Functional Equation, Fixed Point, Numerical Analysis, R5+ Gravity.
This preprint presents a rigorous, self-contained mathematical framework for numerically approximating the Riemann Zeta function using Multiscale Entanglement Renormalization Ansatz (MERA) tensor networks. By formally establishing a unitary equivalence between the MERA coarse-graining operator and a discretized Mellin integral (Lemma 3.1), we analytically derive the theoretical infrared fixed point (Theorem 3.1) to be γtarget=log⁡(2π)2π2≈0.0931γtarget=2π2log(2π)≈0.0931. This derivation leverages the polygamma function evaluated at the center of the critical strip. A numerical implementation with bond dimension χ=2χ=2 achieves this target with an absolute error of 1.3×10−51.3×10−5 and a relative error of +0.0139%+0.0139%. We also provide strict, exponentially converging error bounds based on the singular value decomposition of Hankel/Mellin-type operators (Theorem 3.2). Notably, this mathematically derived fixed point provides the fundamental scale parameter for the Renormalization Group Scale (RG-S) within the R5+ modified gravity model. Keywords: Riemann Zeta Function, MERA, Tensor Networks, Mellin Transform, Renormalization Group, Functional Equation, Fixed Point, Numerical Analysis, R5+ Gravity.
Abstract: The finite algebra A_F = ℂ ⊕ M₂(ℂ) ⊕ M₃(ℂ) at KMS inverse temperature β = 2π provides a complete, finite set of operations for computing physical constants. We classify every operation into five levels and establish a bijective dictionary between spectral operations and standard quantum field theory. The ascending superoperator S on End(ℂ⁶) is the MERA transfer matrix. Its eigenvalues {1, 1/2, 1/3, 1/6} are the anomalous dimensions. Its degeneracies {1, 3, 8, 24} are the channel multiplicities. The first-order channel count Σd = 36 gives tree-level results (33/38 within 1% of measurement). The Schur square Σd² = 650 gives all corrections: the universal MERA retention factor 35/36 = 1 − 1/χ² is the wavefunction renormalisation Z; the Plancherel resolvent Σd²ln(1/(1−λ)) = 137.205 is the spectral zeta function; the Ward anomalies (1−λ_k) = {0, 1/2, 2/3, 5/6} are the Ward–Takahashi identities that ensure one-loop renormalisability (van Nuland & van Suijlekom, JHEP 2022). The spectral truncation error O(1/χ) at finite bond dimension χ = 6 corresponds to UV cutoff artefacts in the Connes–van Suijlekom operator system framework (CMP 2020). Standard model running provides boundary-condition transport from the conformal point (D = 0, ~245 GeV) to the measurement scale. The type-I seesaw is built into the tower formula m_ν = v/2⁴², with the ~10% tree-to-loop gap confirmed as universal finite-χ truncation across five independent hierarchical branching domains. The dictionary contains 19 entries. Each is a bijection: every crystal operation maps to exactly one QFT operation. The five outliers outside 1% each resolve through a specific, standard QFT mechanism — wavefunction renormalisation (α_s/α_w), UV cutoff (horizon entropy), Yukawa running (m_b), and lattice discretisation (neutrinos) — applied to the spectral data of A_F. Three open calculations remain: the van Nuland–van Suijlekom one-loop spectral action evaluated for the specific triple (A_F, H, D); the derivation of Z = 35/36 from first principles for each observable; and the full ΛCDM integration with crystal inputs. The toolbox is finite: 36 channels give tree-level, 650 channels give all corrections, 1296 paths give mixing. No infinite tower of corrections, no ad hoc parameters. The crystal does not replace quantum field theory. It is quantum field theory, written in one algebra. Keywords: spectral action, noncommutative geometry, ascending superoperator, MERA, tensor network, Schur square, wavefunction renormalisation, Ward–Takahashi identities, anomalous dimensions, spectral zeta function, one-loop renormalisability, spectral truncation, operator systems, KMS state, Bisognano–Wichmann, Plancherel measure, finite bond dimension, tree-level closure, seesaw mechanism, crystal toolbox, WACA
The multiscale entanglement renormalization ansatz (MERA) provides a constructive algorithm for realizing wave functions that are inherently scale invariant. Unlike conformally invariant partition functions, however, the finite bond dimension $\ensuremath{\chi}$ of the MERA provides a cutoff in the fields that can be realized. In this paper, we demonstrate that this cutoff is equivalent to the one obtained when approximating a thermal state of a critical Hamiltonian with a matrix product operator (MPO) of finite bond dimension $\ensuremath{\chi}$. This is achieved by constructing an explicit mapping between the isometries of the MERA and the local tensors of the MPO. In terms of energy scales, our results show that a finite bond dimension MERA is equivalent to introducing both an infrared and an ultraviolet scale, characterizing relevant and irrelevant perturbations on the underlying conformal field theory.
CORRECTION NOTICE (v5, 2026-09-06). The central result of v1-v4 is WITHDRAWN. This record has not been peer reviewed. Issued 2026-09-06 by the author (v5 of this record). This notice was added to the archived record; the manuscript itself is unaltered. This record has not been peer reviewed. What is wrong The force exponent reported here as alpha = 3.92 +/- 0.16 is an artifact of the fitting procedure, not a property of the data. The fit was unweighted least squares in linear space applied to a steeply falling profile, which places 95.8 percent of its objective on the single point at the shortest retained separation. Re-fitting the same data with nothing discarded, changing only the error model, gives 2.82 (relative-error weighting) or 2.89 (log-space regression). Across fifteen combinations of standard estimator and fit window the exponent spans 2.53 to 3.92, and the value published here is the largest of them. What the data show instead, and how firmly The retraction is certain. The replacement is not. On the open 24-site chain the same data give alpha = 2.6 +/- 0.2, where the interval is the spread of a plateau in the local logarithmic slope over six consecutive separations, not a statistical error (the bootstrap error on the published fit to the same data is 0.63); the plateau value was reproduced at 2.5921 by a separately written implementation run on the same recorded state (SIM-2026-073). This is a single geometry at a single size. On a 24-site periodic ring, analysed with chord distance as pre-registered, two probes that are both energy-density-like (Delta = 1) give 1.83 and 2.37 against a CFT two-point value of 2 (SIM-2026-074): one size, not converged, disagreeing with each other by 29 percent, and not attributable to a single operator channel. No corrected exponent is established. The exchange channel is an open question. Treat 2.6 +/- 0.2 as the provisional open-chain value, not as a result. Also withdrawn or voided alpha = 3.92 +/- 0.16 (24 sites) and alpha = 3.97 +/- 0.15 (16 sites, exact) from stress-tensor exchange. Was: headline result and its first corroboration. Now: withdrawn. Why: the same estimator artifact: both move in lockstep under it (3.92 to 2.56, 3.97 to 2.75 on dropping the contact point), because they share one observable, one fit window and one fitting routine. scalar exchange (alpha = 2) excluded. Was: stated as excluded at about 12 sigma. Now: withdrawn. Why: the 12 sigma used the covariance error 0.16; at the bootstrap error 0.63 computed before release the exclusion is about 2.2 sigma, and under any window starting at separation 3 the point estimate itself falls to about 2.56. alpha = 4 is universal; alpha = 4.00 at Ising, tricritical Ising and Potts. Was: confirmed at three classes. Now: withdrawn; the operator towers at all three classes stand. Why: the exponent was measured only at Ising and inherited by the other two; the towers are spectra, not forces. third independent evaluation alpha = 4.07 (chi = 3 MERA), and with it the statement that this was the first holographic force measurement from an optimized tensor network state. Was: supporting evidence and a priority claim. Now: both withdrawn. Why: the profile changes sign inside its own fit window and cannot be a fixed-sign power law; the priority claim rests on that evaluation alone and was never gated against prior work. the discriminating test: stress tensor (Z_j Z_j+1) versus spin primary (X_j). Was: registered killer test. Now: void. Why: at the critical point these are one operator, related by exact Kramers-Wannier self-duality. Measured on the periodic ring: identical expectation values to 1e-15 and identical force profiles to 1e-13 at 14 sites, 2.4e-9 at 24 sites (SIM-2026-074). Open boundaries break the duality, which is why the open-chain run of this test (SIM-2026-073) showed a spurious difference between the two arms. The test compared the energy density to itself. Neither label was correct: Z_j Z_j+1 is the energy density, not the stress tensor; X_j is its dual, not the spin primary. transport time T = (pi/2) l_AdS to 1.2 percent as a geometric prediction. Was: confirmed. Now: measurement stands (T = 2.239 at N = 16); geometric reading not established. Why: l_AdS cancels in the cited derivation, which also gives pi rather than the pi/2 the manuscript compares against; the chi = 3 estimate of l_AdS gives a 57 percent miss; the transport driver never applies a MERA. The true spin primary Z_j was never in the design, and the direct replacement test is blocked: Z_j is odd under the Z2 symmetry, so the energy-difference observable used throughout this work subtracts across symmetry sectors and is undefined for it. A replacement test using two Z2-even probes has been pre-registered (PREREGISTRATION-T.md, committed before its driver) and run; its result is deliberately not part of this correction and will appear in a later version. Not previously reported The published fit has a reduced chi-squared of 6.008 on 9 degrees of freedom, and a bootstrap standard deviation of 0.63 rather than the quoted 0.16. Both were computed before release and omitted. What is unaffected Of the claims in this record: the two-phase optimization method, the chi = 3 criticality threshold, the exact-diagonalization gates, the measured OPE coefficient C_sigma-sigma-epsilon = 0.50 +/- 0.03 against the exact 1/2, and the conformal towers at all three classes. Not part of this record, mentioned only because they are sometimes read together with it: the later C_TTT stress-tensor self-coupling measurements (SIM-2026-067/068/069, dated after this record) and the matched off-critical ring control (SIM-2026-074) are unaffected by this correction and make no claim about the force exponent. How the error was found and what this record contains The error was found by the author's own internal adversarial audit, not by a reviewer; the journal submission it was found after is closed, and no review of this work has taken place. Files in this record: the full re-analysis (SIM-2026-072.json, driver, figure), the pre-registered chain test (SIM-2026-073.json) and ring test (SIM-2026-074.json) with both pre-registrations as committed before the runs, the runnable drivers with their shared module and the negative test suite for their validation gates, and the error register (ERRATA.md). Concept DOI 10.5281/zenodo.19660302 always resolves to the current version. The abstract of the original manuscript follows unchanged, for the record. Every claim in it about the force exponent should be read against the correction above. We report results from classical simulation of critical spin chains that establish quantitative physical content in the MERA-AdS correspondence. First, the interaction energy between boundary-primary excitations follows a power law F(dr) ~ dr^(-alpha) with three consistent evaluations yielding alpha = 3.92 +/- 0.16 (24q), 3.97 +/- 0.15 (16q), and 4.07 +/- 0.18 (16q chi=3 MERA), all consistent with stress-tensor exchange (alpha = 2*Delta_T = 4) and inconsistent with scalar exchange (alpha = 2); the exponent is explained by the Ising CFT fusion rule eps x eps = [1]. Because Delta_T = 2 is fixed by conformal symmetry in every unitary 2D CFT, alpha = 4 is universal: we confirm this at two further fixed points - the tricritical Ising model (c = 7/10), where a scale-invariant MERA recovers the full M(4,5) primary spectrum (Delta_sigma = 0.092, Delta_eps = 0.204, Delta_T = 1.999), and the three-state Potts model (c = 4/5), where a Z3-symmetric MERA at chi = 12 recovers Delta_sigma = 0.134, Delta_eps = 0.805, Delta_T = 2.001 - giving alpha = 4.00 at all three universality classes. Second, the holographic transport time matches the AdS_3 null geodesic prediction T = (pi/2) l_AdS to 1.2% at N = 16 sites on a time-grid-converged measurement, with the residual shrinking from 2.8% at N = 12. Third, we introduce a two-phase optimization method (ALM + Z2 projection) for Z2-symmetric MERA ground states at 24 qubits. Fourth, extending the MERA bond dimension from chi=2 to chi=3 overcomes the expressibility barrier, achieving 0.15% energy error at 16 qubits and enabling the first holographic force measurement from an optimized tensor network state. The OPE coefficient is also measured directly from the fixed-point tensors: C_sigma-sigma-eps = 0.50 +/- 0.03 against the exact Ising value 1/2. A hardware path is designed for Quantinuum Helios (57 CNOTs, 0.57% logical error). Version 4 changes (relative to v3): Transport: N=16 residual restated as the time-grid-converged value of -1.21% (dt = 1e-3), with -2.78% at N=12 showing the finite-size trend; the 0.72% figure in v3 was an artifact of the dt = 0.05 integration grid and is now labeled as such. The N=24 transport point carries a new reliability caveat (double-lobe trace structure). The l_AdS comparison table is recomputed against the converged transport time (best fit 1.425, deviation 1.21%). Effective-mass derivation rewritten in a consistent AdS3 convention m^2 l^2 = Delta(Delta-2); m_eff = 0.693 = 1/l_AdS follows from Delta = 1 exactly. NEW: universality section extended to a third fixed point - a Z3-symmetric MERA at chi = 12 recovers the three-state Potts (c = 4/5) tower, giving alpha = 4.00 at all three universality classes (c = 1/2, 7/10, 4/5). NEW: the OPE coefficient C_sigma-sigma-eps = 0.50 +/- 0.03 measured directly from the fixed-point tensors, with C_sigma-sigma-sigma = 0 at machine precision as the Z2 structural check. Supplemental table row SIM-007 annotated (record lost; script preserved); data availability expanded; Zenodo DOI added to the front matter.
This thesis reviews the multiscale entanglement renormalisation ansatz or MERA, a numerical tool for the study of quantum many-body systems and a discrete realisation of the AdS/CFT duality. The thesis covers an introduction to the necessary background concepts of entanglement, entanglement entropy and tensor network states, the structure and main features of MERA and its applications in condensed matter theory and holography. Also covered are details on the algorithmic implementation of MERA and some of its generalisations and extensions. MERA belongs to a class of variational ansatze for quantum many-body states known as tensor network states. It is especially well-suited for the study of scale invariant critical points. MERA is based on a real-space renormalisation group procedure called entanglement renormalisation, designed to systematically handle entanglement at different length scales along the coarse-graining flow. Entanglement renormalisation has be used for example to efficiently describe Kitaev states of the toric code, the prime example of topological order, and numerically study the ground state of the highly frustrated spin-1/2 Heisenberg model on a kagome lattice and various other one- and two-dimensional lattice models. The geometric and causal structure of MERA, which underlies its effectiveness as a numerical tool, also makes it a discrete version of the AdS/CFT duality. This duality describes a conformal field theory by a gravity theory in a higher dimensional space, and vice versa. The duality is manifest in the scaling of entanglement entropy in MERA, which is governed by a law highly analogous to the Ryu-Takayanagi formula for holographic entanglement entropy, in the connection between thermal states and a black-hole-like MERA and in the connection between correlation functions and holographic geodesics in a scale invariant MERA. The aim of this thesis is to lead the reader to an understanding of what MERA is, how it works and how it can be used. MERA's core features and uses are presented in a comprehensive and explicit way, and a broad view of possible applications and further directions is given. Plenty of references are also offered to direct the reader to further research on how MERA may relate to his/her interests.