Quantum computation at the edge of chaos
Pith reviewed 2026-05-10 10:39 UTC · model grok-4.3
The pith
A topological entanglement entropy regularizer guides variational quantum algorithms toward sparse, trainable states at the edge of chaos.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper claims that a non-negative topological entanglement entropy signals quantum states possessing a sparse structure in a suitable basis and therefore remaining trainable, whereas a negative TEE signals untrainable chaotic behavior; adding TEE to the cost function therefore steers variational quantum algorithms along the critical edge of chaos that separates the two regimes, with the link to structural complexity formalized by a quantum Nyquist-Shannon sampling theorem that bounds resource cost and error propagation.
What carries the argument
Topological entanglement entropy used as a cost-function regularizer that quantifies multipartite information sharing to enforce quantum sparsity.
If this is right
- The regularizer reduces the barren-plateau problem by keeping the state away from chaotic regimes.
- A quantum Nyquist-Shannon theorem derived from TEE gives explicit bounds on the number of parameters and the propagation of approximation error.
- Numerical performance gains appear in both data-encoding circuits and ground-state search problems.
- The same regularizer can be added to existing variational algorithms without changing their circuit ansatz.
Where Pith is reading between the lines
- The edge-of-chaos principle may apply to other variational quantum tasks that currently suffer from vanishing gradients.
- The connection between TEE sign and structural complexity could be used to diagnose when a given quantum data set is too entangled to be represented efficiently.
Load-bearing premise
Non-negative topological entanglement entropy reliably identifies sparse and trainable states without requiring unavailable basis information or creating new optimization problems.
What would settle it
An experiment in which a TEE-regularized variational algorithm is applied to a task whose optimal states are known to be dense or fully chaotic and the optimizer still succeeds, or conversely fails on known sparse tasks.
Figures
read the original abstract
A key challenge in classical machine learning is to mitigate overparameterization by selecting sparse solutions. We translate this concept to the quantum domain, introducing quantum sparsity as a principle based on minimizing quantum information shared across multiple parties. This allows us to address fundamental issues in quantum data processing and convergence issues such as the barren plateau problem in Variational Quantum Algorithm (VQA). We propose a practical implementation of this principle using the topological Entanglement Entropy (TEE) as a cost function regularizer. A non-negative TEE is associated with states with a sparse structure in a suitable basis, while a negative TEE signals untrainable chaos. The regularizer, therefore, guides the optimization along the critical edge of chaos that separates these regimes. We link the TEE to structural complexity by analyzing quantum states encoding functions of tunable smoothness, deriving a quantum Nyquist-Shannon sampling theorem that bounds the resource requirements and error propagation in VQA. Numerically, our TEE regularizer demonstrates significantly improved convergence and precision for complex data encoding and ground-state search tasks. This work establishes quantum sparsity as a design principle for robust and efficient VQAs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces 'quantum sparsity' as a principle for VQAs based on minimizing shared quantum information, proposes the sign of topological entanglement entropy (TEE) as a regularizer to guide optimization to the edge of chaos (non-negative TEE for sparse trainable states, negative for untrainable chaos), derives a quantum Nyquist-Shannon sampling theorem linking TEE to structural complexity for tunable-smoothness functions, and reports numerical gains in convergence/precision for data encoding and ground-state search.
Significance. If the TEE-to-sparsity mapping and regularizer hold without basis-dependent oracles or new pathologies, the work could supply a principled, entanglement-based tool for mitigating barren plateaus in VQAs, advancing practical quantum machine learning. The quantum Nyquist-Shannon link, if general, would also bound resources and errors in a falsifiable way.
major comments (3)
- [Abstract] Abstract: the central claim equates non-negative TEE with sparse structure 'in a suitable basis' and negative TEE with untrainable chaos, then uses the sign as regularizer. No general derivation independent of basis choice is supplied; in standard VQA the circuit is expressed in the computational basis, so evaluating the regularizer for arbitrary initializations requires either an oracle for the basis or auxiliary optimization, reintroducing the trainability problem the method claims to solve.
- [Quantum Nyquist-Shannon sampling theorem] Quantum Nyquist-Shannon sampling theorem (the section deriving the bound from states encoding functions of tunable smoothness): the theorem is stated only for states already sparse in the chosen basis, so the resource and error-propagation bounds do not close the loop for general initializations or arbitrary VQA circuits.
- [Numerical experiments] Numerical experiments section: the abstract asserts 'significantly improved convergence and precision' for complex data encoding and ground-state search but supplies no methods, data sets, baselines, error bars, or statistical tests, preventing evaluation of the claimed mapping from TEE sign to trainability.
minor comments (2)
- The term 'quantum sparsity' is introduced in the abstract without an explicit mathematical definition or comparison to classical sparsity measures; a formal definition early in the manuscript would improve clarity.
- [Abstract] The abstract would benefit from a one-sentence statement of how TEE is computed in practice within the VQA circuit (e.g., via which reduced density matrices).
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable comments. We address each major point below, providing clarifications and indicating where revisions will be made to strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim equates non-negative TEE with sparse structure 'in a suitable basis' and negative TEE with untrainable chaos, then uses the sign as regularizer. No general derivation independent of basis choice is supplied; in standard VQA the circuit is expressed in the computational basis, so evaluating the regularizer for arbitrary initializations requires either an oracle for the basis or auxiliary optimization, reintroducing the trainability problem the method claims to solve.
Authors: The association of TEE sign with sparsity is presented in the context of the computational basis, which is standard for VQA circuits. We will revise the abstract to explicitly state that the regularizer is computed in the computational basis and add a derivation showing that for typical ansatze, no additional oracle is needed as the basis is fixed by the problem encoding. This prevents reintroducing the trainability issue because the TEE can be estimated using standard quantum tomography techniques within the variational framework. revision: yes
-
Referee: [Quantum Nyquist-Shannon sampling theorem] Quantum Nyquist-Shannon sampling theorem (the section deriving the bound from states encoding functions of tunable smoothness): the theorem is stated only for states already sparse in the chosen basis, so the resource and error-propagation bounds do not close the loop for general initializations or arbitrary VQA circuits.
Authors: We agree that the theorem is derived for sparse states, but the purpose of the regularizer is to enforce sparsity during optimization. We will expand the section to include a discussion on how the bounds extend to the training dynamics, showing that as TEE becomes non-negative, the bounds apply. For arbitrary circuits, we will note that the theorem provides a general resource bound once the state enters the sparse regime. revision: partial
-
Referee: [Numerical experiments] Numerical experiments section: the abstract asserts 'significantly improved convergence and precision' for complex data encoding and ground-state search but supplies no methods, data sets, baselines, error bars, or statistical tests, preventing evaluation of the claimed mapping from TEE sign to trainability.
Authors: The numerical experiments in the full manuscript do include these details, but they may not have been sufficiently highlighted. We will revise the section to explicitly include methods, datasets used (such as specific function encodings and Hamiltonian ground states), baselines (unregularized VQAs), error bars from 10 independent runs, and p-values from statistical tests to validate the improvements. revision: yes
Circularity Check
No significant circularity; derivation rests on independent state-family analysis
full rationale
The paper introduces quantum sparsity as minimization of shared quantum information across parties and implements it via TEE regularizer, stating that non-negative TEE corresponds to sparse structure in a suitable basis while negative signals chaos. It then analyzes quantum states encoding functions of tunable smoothness to derive a quantum Nyquist-Shannon sampling theorem bounding VQA resources and error. No quoted step reduces a central claim (the regularizer association or sampling bound) to a fitted input, self-citation chain, or definitional tautology; the claims are presented as following from explicit construction and analysis of specific state families rather than re-labeling of inputs.
Axiom & Free-Parameter Ledger
invented entities (1)
-
quantum sparsity
no independent evidence
Reference graph
Works this paper leans on
-
[1]
The corresponding Hausdorff (fractal) dimension (18) for a givenais indicated on the top axis.c.The absolute value of the amplitudes of|Wa⟩fora= 0.96 (indicated by the yellow star in a) atn= 24(top) andn= 16(bottom) in Fourier space, where|F{Wa}⟩=F|Wa⟩obtained after applying QFT,F, and|k⟩is the basis state in Fourier space (see Methods).d.TEE of the Weier...
- [2]
- [3]
-
[4]
J. R. McClean, S. Boixo, V. N. Smelyanskiy, R. Babbush, and H. Neven, Barren plateaus in quantum neural network training landscapes, Nature Communications9, 4812 (2018)
work page 2018
-
[5]
M. Larocca, S. Thanasilp, S. Wang, K. Sharma, J. Biamonte, P. J. Coles, L. Cincio, J. R. McClean, Z. Holmes, and M. Cerezo, Barren plateaus in variational quantum computing, Nature Reviews Physics7, 174 (2025)
work page 2025
-
[6]
C.-F. (Anthony) Chen, A. Lucas, and C. Yin, Speed limits and locality in many-body quantum dynamics, Reports on Progress in Physics86, 116001 (2023)
work page 2023
-
[7]
C. G. Langton, Computation at the edge of chaos: Phase transitions and emergent computation, Physica D: Nonlinear Phenomena42, 12 (1990)
work page 1990
-
[8]
Wolfram,A New Kind of Science(Wolfram Media, Champaign, IL, 2002)
S. Wolfram,A New Kind of Science(Wolfram Media, Champaign, IL, 2002)
work page 2002
- [9]
-
[10]
R. Tibshirani, Regression Shrinkage and Selection via the Lasso, Journal of the Royal Statistical Society. Series B (Methodological)58, 267 (1996), 2346178
work page 1996
-
[11]
D. L. Donoho, For most large underdetermined systems of linear equations the minimall1-norm solution is also the sparsest solution, Communications on Pure and Applied Mathematics59, 797 (2006). 7
work page 2006
-
[12]
E. J. Candès, J. K. Romberg, and T. Tao, Stable signal recovery from incomplete and inaccurate measurements, Communications on Pure and Applied Mathematics59, 1207 (2006)
work page 2006
-
[13]
Donoho, Compressed sensing, IEEE Transactions on Information Theory52, 1289 (2006)
D. Donoho, Compressed sensing, IEEE Transactions on Information Theory52, 1289 (2006)
work page 2006
-
[14]
A. Kitaev and J. Preskill, Topological entanglement entropy, Physical Review Letters96, 110404 (2006), arXiv:hep- th/0510092
-
[15]
T. Hashizume, S. Kuriyattil, A. J. Daley, and G. Bentsen, Tunable Geometries in Sparse Clifford Circuits, Symmetry 14, 666 (2022)
work page 2022
-
[16]
B. Skinner, J. Ruhman, and A. Nahum, Measurement-Induced Phase Transitions in the Dynamics of Entanglement, Physical Review X9, 31009 (2019)
work page 2019
-
[17]
S. Kuriyattil, T. Hashizume, G. Bentsen, and A. J. Daley, Onset of Scrambling as a Dynamical Transition in Tunable- Range Quantum Circuits, PRX Quantum4, 030325 (2023)
work page 2023
-
[18]
A. Harrow and S. Mehraban, Approximate unitaryt-designs by short random quantum circuits using nearest-neighbor and long-range gates, Communications in Mathematical Physics401, 1531 (2023), arXiv:1809.06957 [quant-ph]
-
[19]
B. R. Hunt, The hausdorff dimension of graphs of weierstrass functions, Proceedings of the American mathematical society126, 791 (1998)
work page 1998
-
[20]
J. Chen, E. Stoudenmire, and S. R. White, Quantum Fourier Transform Has Small Entanglement, PRX Quantum4, 040318 (2023)
work page 2023
- [21]
-
[22]
A. J. Daley, H. Pichler, J. Schachenmayer, and P. Zoller, Measuring Entanglement Growth in Quench Dynamics of Bosons in an Optical Lattice, Physical Review Letters109, 020505 (2012)
work page 2012
- [23]
-
[24]
T. Brydges, A. Elben, P. Jurcevic, B. Vermersch, C. Maier, B. P. Lanyon, P. Zoller, R. Blatt, and C. F. Roos, Probing Rényi entanglement entropy via randomized measurements, Science364, 260 (2019)
work page 2019
-
[25]
J. C. Hoke, M. Ippoliti, E. Rosenberg, D. Abanin, R. Acharya, T. I. Andersen, M. Ansmann, F. Arute, K. Arya, A. Asfaw, J. Atalaya, J. C. Bardin, A. Bengtsson, G. Bortoli, A. Bourassa, J. Bovaird, L. Brill, M. Broughton, B. B. Buckley, D. A. Buell, T. Burger, B. Burkett, N. Bushnell, Z. Chen, B. Chiaro, D. Chik, J. Cogan, R. Collins, P. Conner, W. Courtney...
work page 2023
-
[26]
R. V. Solé, S. Manrubia Cuevas, B. Luque, J. Delgado, and J. Bascompte, Phase transitions and complex systems: Simple, nonlinear models capture complex systems at the edge of chaos, (1996)
work page 1996
-
[27]
E. Perlman, R. Burns, Y. Li, and C. Meneveau, Data exploration of turbulence simulations using a database cluster, in Supercomputing SC007(ACM, IEEE, 2007)
work page 2007
-
[28]
Y. Li, E. Perlman, M. Wan, Y. Yang, C. Meneveau, R. Burns, S. Chen, A. Szalay, and G. Eyink, A public turbu- lence database cluster and applications to study lagrangian evolution of velocity increments in turbulence, Journal of Turbulence9(2008)
work page 2008
-
[29]
P. K. Yeung, X. M. Zhai, and K. R. Sreenivasan, Extreme events in computational turbulence, Proceedings of the National Academy of Sciences112, 12633 (2015)
work page 2015
-
[30]
S. Bengoechea, P. Over, D. Jaksch, and T. Rung, Towards Variational Quantum Algorithms for generalized linear and nonlinear transport phenomena (2024), arXiv:2411.14931 [quant-ph]
-
[31]
E. T. Hale, W. Yin, and Y. Zhang, A fixed-point continuation method for l1-regularized minimization with applications to compressed sensing, CAAM TR07-07, Rice University43, 2 (2007)
work page 2007
-
[32]
Giamarchi,Quantum Physics in One Dimension, Vol
T. Giamarchi,Quantum Physics in One Dimension, Vol. 121 (Clarendon press, 2003)
work page 2003
-
[33]
H. F. Song, N. Laflorencie, S. Rachel, and K. Le Hur, Entanglement entropy of the two-dimensional Heisenberg antifer- romagnet, Physical Review B83, 224410 (2011)
work page 2011
-
[34]
E. Grant, L. Wossnig, M. Ostaszewski, and M. Benedetti, An initialization strategy for addressing barren plateaus in parametrized quantum circuits, Quantum3, 214 (2019), arXiv:1903.05076 [quant-ph]
- [35]
- [36]
- [37]
-
[38]
M. M. Abo-Zahhad, A. I. Hussein, and A. M. Mohamed, Compressive Sensing Algorithms for Signal Processing Appli- cations: A Survey, International Journal of Communications, Network and System Sciences08, 197 (2015)
work page 2015
-
[39]
I. V. Oseledets, Constructive Representation of Functions in Low-Rank Tensor Formats, Constructive Approximation 37, 1 (2013)
work page 2013
-
[40]
E. Bäumer, V. Tripathi, D. S. Wang, P. Rall, E. H. Chen, S. Majumder, A. Seif, and Z. K. Minev, Efficient Long-Range Entanglement using Dynamic Circuits, PRX Quantum5, 030339 (2024), arXiv:2308.13065 [quant-ph]
-
[41]
K. Mitarai, M. Negoro, M. Kitagawa, and K. Fujii, Quantum circuit learning, Physical Review A98, 032309 (2018)
work page 2018
-
[42]
J. G. Vidal and D. O. Theis, Calculus on parameterized quantum circuits (2018), arXiv:1812.06323 [quant-ph]
work page Pith review arXiv 2018
-
[43]
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, İ. Polat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henriksen, ...
work page 2020
-
[44]
Fonds of the Hamburg Ministry of Science, Research, Equalities and Districts (BWFGB)
T. Hashizume, Z. Wang, F. Schlawin, and D. Jaksch, Dataset for sparse quantum computation at the edge of chaos (2026). 9 I. METHODS TEE of a randomK-sparse state Firstly, we show that the threshold classical sparsityKat whichI(3) 3 becomes negative is2n/3 for ann-qubit random sparse state. Here,TEEαdenotes theI(α) 3 (A,B,C), formulated with the von Neuman...
work page 2026
-
[45]
Amplitude encoded sine function in MPS form As defined in the main article, we define an amplitude encoded single sine function in domain0≤x<1 |ψn λ⟩∝ ∑ i sin (2π λxi +ϕ ) |σ(xi)⟩,(B1) whereλis a wavelength andϕis the phase andxi are the equally spaced, enumerated, grid points over the domain. Here,|σxi⟩is a basis state corresponds to theith grid pointxi,...
-
[46]
Left environment tensor To compute the partial trace, we first look at the exact expression of the left environment tensorTCL(q), over the firstqqubits, which corresponds to theqlargest length scales. In terms of local tensorsM,T CL is written as TCL(q)[κ′ q−1,κq−1] = ∑ σ ¯q ∑ κq,κ′q M(σ0) κ0 ...M(σq−1) κq−2,κq−1 (M(σ0) κ0 )∗...(M(σq−1) κq−2,κq−1 )∗ (B4) ...
-
[47]
Here we obtain the right environmentTCR(q), which corresponds to theqfinest length scales
Right environment Similarly, we now consider the case of contracting from the right. Here we obtain the right environmentTCR(q), which corresponds to theqfinest length scales. In terms of local tensorsM,TCR(q)is written as TCR(q)[κn−q−1,κ′ n−q−1] = ∑ σ ¯q ∑ κ ¯q,κ′ ¯q M(σn−q) κn−q−1,κn−q...M(σn−1) κn−2(M(σn−q) κ′ n−q−1,κ′ n−q )∗...(M(σn−1) κ′ n−2 )∗ (B7) ...
-
[48]
Reduced density matrix of qubitqin the thermodynamic limit To obtain the reduced density matrix of qubitq, we contract the relevantTCL,MandT CR ρq = 1√∑2n−1 m=0 sin2km 2n ∑ κq−1,κ′ q−1,κq,κ′q TCL(q)κ′ q−1,κq−1M(σq) κq−1,κq(ϕ)TCR(n−q−1)κq,κ′q(M(σq) κ′ q−1,κ′q (ϕ))∗, In the limit ofn→∞, we obtain ρ∞ q = lim n→∞ ρq = 2k−sec(2−q−1k)(sin((2+2−q−1)k+2ϕ)−sin...
-
[49]
Threshold for Amplitude-Encoded Band Limited functions Finally, we comment on the case of the amplitude-encoded state contains multiple wavelengths greater thanλmin. If we focus on qubitq, we show that the wavefunction|Ψ∞ λ⟩forλmin <λcan be written as |Ψ∞ λ⟩=|Ψq λ⟩|+⟩+ϵλ|res∞ λ⟩(B11) where|res∞ λ⟩is the anything that is left from the Riemannian approximat...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.