REVIEW 3 major objections 4 minor 21 references
Optimal packing of attractor states in neural representations
T0 review · 3 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read For uniform and cyclic Markov chains, the optimal packing of neural attractor states is solved under a symmetry ansatz, yielding a simplex for uniform statistics and a generally non-square quadrilateral for four-state cycles.
desk verdict A fresh normative framework for neural geometry with a genuinely non-square M=4 result, but the symmetry-inheritance claim is over-sold and an appendix arithmetic slip needs fixing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
The authors write a cost function that combines three pieces: how well the brain's transitions match the environment's, how accurately the brain can decode which state it is in, and how much activity the states require. They show the cost only depends on the pairwise distances between states, not on where the whole cluster sits. Then they use symmetry: if the environment's transition rules look the same under some rearrangements, they look for arrangements of states that also look the same under those rearrangements. For a fully uniform environment, symmetry forces all states to be equally far apart, forming a simplex, like the corners of a triangle. For a four-state ring, symmetry leaves two distances free: the ring edges and the diagonals. The optimum generally does not form a perfect square, which is a concrete, unexpected prediction.
The results are exact only within the symmetric ansatz, and the paper does not prove this ansatz is the true global optimum. One formula in the uniform case also has a small counting error. Still, the framework gives a clean way to think about why ring attractors and grid-like codes have the geometry they do, and it offers toy predictions that could be tested in simulations or data.
Extended reading notes
Core claim
For Markov chains with uniform transition statistics, the optimal arrangement of M attractor states is an (M-1)-simplex with all pairwise distances equal; for M=4 cyclic statistics, the optimal arrangement is a D4-symmetric quadrilateral whose edge length d and diagonal L generally do not satisfy L/d = sqrt(2) (Sections 4.1, 4.2). If the paper is correct, symmetries of environmental transition statistics are reflected in the geometry of optimal neural representations.
Load-bearing premise
The global minimum of the objective J (Eq. 8) is assumed to inherit the symmetries of the Markov chain, via an unproved application of the Purkiss principle. The paper states the problem is 'probably not convex', that it must 'settle for something weaker', and that a rigorous analysis is not pursued; if the true optimum is not symmetric, the reported packings are only symmetry-constrained stationary points, not the optimal packing claimed in the title and abstract. Location: Section 3, paragraph beginning 'For convex optimization problems' and the sentence 'A more rigorous analysis of Eq. 11 may be able to show that Waterhouse's precise formulation of the Purkiss principle applies, although we do not pursue such an analysis here.'
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript formulates a normative toy model for arranging M internal attractor states in a neural representation. Environment statistics are modeled by a Markov chain, internal transitions by distance-dependent rates, and encoding/decoding by a noisy symmetric channel; the objective J (Eq. 8) trades off the KL divergence between the desired and actual internal transition distributions against a firing-rate regularization. After a change of variables, J depends only on the mean location and pairwise distances (Eqs. 10–11), and the mean is set to zero. Relying on a Purkiss-style symmetry argument, the paper restricts to symmetric configurations and claims to solve two cases: uniform statistics, for which all distances are equal (an (M−1)-simplex), and M=4 cyclic statistics, for which the optimum is a D4-symmetric quadrilateral with two distances d and L. A large-b approximation for general cyclic topologies is also sketched. Code is provided.
Significance. The paper's modeling pipeline is clear and the reduction to pairwise distances is a useful simplification; the uniform-prior decoding step and the symmetry of the objective under Markov-chain permutations are correctly derived. The two-parameter model yields concrete, falsifiable geometric predictions, and the availability of code supports reproducibility. However, the central optimality claim is not established: the symmetry restriction is an assumption about the global minimizer, not a proof, and there are concrete arithmetic and geometric errors in key results. In its current form the contribution is best read as a symmetry-constrained analysis of a plausible objective rather than as a proof of optimal packing.
major comments (3)
- [Section 3, Sections 4.1–4.2] The paper's central claim that it has found the optimal packing is not supported. The text explicitly says Eq. (8) is 'probably not convex', that the authors 'settle for something weaker' in the spirit of the Purkiss principle, and that a rigorous analysis of Eq. (11) 'is not pursued'. Sections 4.1 and 4.2 then optimize only the symmetry-constrained ansatz (all distances equal, or two distances d and L under D4). Because a symmetric stationary point of a nonconvex function need not be a global minimum, and because the objective is defined over all configurations in R^D, the results as stated ('optimal packing', 'optimal solution') go beyond what is proven. This is load-bearing for the abstract and title. The authors should either prove global optimality for the two solved cases (e.g., by a convexity or monotonicity argument on the reduced objectives, combined with explicit feasible regions) or explicitly reframe the contribution as the optimal symmetric packing.
- [Appendix B, Section 4.1] The coefficient of e^{-d^2} in the uniform-case p_int and in the displayed J is wrong. For x≠y, the number of ordered pairs (a,b) with a≠x, b≠y, and a≠b is (M−1)^2 − (M−2) = M^2 − 3M + 3, not M(M−1) − 2(M−2) + 1 = M^2 − 3M + 5. Direct counting for M=2 gives one cross term, whereas the printed formula gives three. Consequently the objective plotted in Fig. 2 and the associated phase-transition statements must be recomputed; this is an arithmetic error in a central equation, not a typo.
- [Section 4.2] The claimed large-b degenerate solution 'd = 0 and L ≠ 0' is not realizable as a Euclidean configuration. If the four side distances d12=d23=d34=d41 are zero, then all four vertices coincide and the diagonals L=d13=d24 must also be zero. The feasible line-degenerate configuration for the four-cycle is the opposite limiting case, L=0 with d>0 (two pairs of states coalesce). The text and Fig. 3 panels that describe this phase therefore need correction, and the optimization should explicitly enforce the Euclidean embedding constraints on (d,L), at minimum the triangle inequalities for the face triangles.
minor comments (4)
- [Section 4.3] The section is titled 'Optimal packing for cyclic topology' but, for M>4, it only derives an approximate large-b system of equations and states that it cannot be solved analytically. The section should be retitled or clearly scoped as an approximate analysis, since the abstract's 'cyclic transition statistics' case is only exactly solved for M=4.
- [Section 2] The statement that setting p(x|x)=0 is 'without loss of generality' needs a sentence of justification, since deleting self-loops generally changes transition probabilities on a fixed time scale.
- [Appendix D, Eqs. (33)–(35)] The exponents in the sums after the derivative appear to mix d_k^2 and d_k^2/2; please recheck these equations against the bracket in Eq. (32) and the derivative of Z.
- [Figure 3b] The caption says the objective is plotted versus d and L, but the main text describes a two-variable landscape; please specify whether the figure shows a contour or a surface and state the fixed parameter values.
Assumptions & free parameters
free parameters (2)
- b (encoding bias)
- alpha (regularization strength)
assumptions (4)
- domain assumption Internal transition probabilities depend only on pairwise distance via q(z_y|z_x) = (1-delta_xy) exp(-D^2/2)/Z(z_x) (Eq. 2)
- domain assumption Encoding noise is isotropic in the noise-whitened metric and controlled by a single bias b (Eq. 3)
- domain assumption Decoding uses a uniform prior over states (Eq. 4)
- domain assumption The global optimum of J inherits the Markov chain's symmetries (Purkiss principle, Waterhouse 1983)
Cite this review
Pith. "Pith review of Optimal packing of attractor states in neural representations." pith.science (2026). https://pith.science/paper/FQIL7U3Q
@misc{pith2026250412429,
author = {Pith},
title = {Pith review of: Optimal packing of attractor states in neural representations},
year = {2026},
howpublished = {\url{https://pith.science/paper/FQIL7U3Q}},
note = {Machine review of arXiv:2504.12429}
}
read the original abstract
Animals' internal states reflect variables like their position in space, orientation, decisions, and motor actions -- but how should these internal states be arranged? Internal states which frequently transition between one another should be close enough that transitions can happen quickly, but not so close that neural noise significantly impacts the stability of those states, and how reliably they can be encoded and decoded. In this paper, we study the problem of striking a balance between these two concerns, which we call an `optimal packing' problem since it resembles mathematical problems like sphere packing. While this problem is generally extremely difficult, we show that symmetries in environmental transition statistics imply certain symmetries of the optimal neural representations, which allows us in some cases to exactly solve for the optimal state arrangement. We focus on two toy cases: uniform transition statistics, and cyclic transition statistics. Code is available at https://github.com/john-vastola/optimal-packing-neurreps23 .
Figures
Reference graph
Works this paper leans on
-
[1]
Keinath, Xue-Xin Wei, and Mark P
Zaki Ajabi, Alexandra T. Keinath, Xue-Xin Wei, and Mark P. Brandon. Population dynamics of head-direction neurons during drift and reorientation. Nature, 615 0 (7954): 0 892--899, Mar 2023. ISSN 1476-4687. doi:10.1038/s41586-023-05813-2. URL https://doi.org/10.1038/s41586-023-05813-2
-
[2]
C. Barry and N. Burgess. Neural mechanisms of self-location. Current Biology, 24 0 (8): 0 R330--R339, 2014. ISSN 0960-9822. doi:https://doi.org/10.1016/j.cub.2014.02.049. URL https://www.sciencedirect.com/science/article/pii/S0960982214002176
-
[3]
Yoram Burak and Ila R. Fiete. Fundamental limits on persistent activity in networks of noisy neurons. Proceedings of the National Academy of Sciences, 109 0 (43): 0 17645--17650, 2012. doi:10.1073/pnas.1117386109. URL https://www.pnas.org/doi/abs/10.1073/pnas.1117386109
-
[4]
Zheng, and Kevin Chen-Chuan Chang
HongYun Cai, Vincent W. Zheng, and Kevin Chen-Chuan Chang. A comprehensive survey of graph embedding: Problems, techniques, and applications. IEEE Transactions on Knowledge and Data Engineering, 30 0 (9): 0 1616--1637, 2018. doi:10.1109/TKDE.2018.2807452
arXiv 2018
-
[5]
Neural representations of motor plans, desired trajectories, and controlled objects
Paul Cisek. Neural representations of motor plans, desired trajectories, and controlled objects. Cognitive Processing, 6 0 (1): 0 15--24, Mar 2005. ISSN 1612-4790. doi:10.1007/s10339-004-0046-7. URL https://doi.org/10.1007/s10339-004-0046-7
-
[6]
The sphere packing problem in dimension 24
Henry Cohn, Abhinav Kumar, Stephen Miller, Danylo Radchenko, and Maryna Viazovska. The sphere packing problem in dimension 24 . Annals of Mathematics, 185 0 (3): 0 1017 -- 1033, 2017. doi:10.4007/annals.2017.185.3.8. URL https://doi.org/10.4007/annals.2017.185.3.8
-
[7]
Driscoll, Lea Duncker, and Christopher D
Laura N. Driscoll, Lea Duncker, and Christopher D. Harvey. Representational drift: Emerging theories for continual learning and experimental future directions. Current Opinion in Neurobiology, 76: 0 102609, 2022. ISSN 0959-4388. doi:https://doi.org/10.1016/j.conb.2022.102609. URL https://www.sciencedirect.com/science/article/pii/S0959438822001039
arXiv 2022
-
[8]
A. Aldo Faisal, Luc P. J. Selen, and Daniel M. Wolpert. Noise in the nervous system. Nature Reviews Neuroscience, 9 0 (4): 0 292--303, Apr 2008. ISSN 1471-0048. doi:10.1038/nrn2258. URL https://doi.org/10.1038/nrn2258
doi:10.1038/nrn2258 2008
Show all 21 references
-
[9]
Gold and Michael N
Joshua I. Gold and Michael N. Shadlen. The neural basis of decision making. Annual Review of Neuroscience, 30 0 (1): 0 535--574, 2007. doi:10.1146/annurev.neuro.29.051605.113038. URL https://doi.org/10.1146/annurev.neuro.29.051605.113038. PMID: 17600525
2007
-
[10]
Reaction-rate theory: fifty years after K ramers
Peter H\"anggi, Peter Talkner, and Michal Borkovec. Reaction-rate theory: fifty years after K ramers. Rev. Mod. Phys., 62: 0 251--341, Apr 1990. doi:10.1103/RevModPhys.62.251. URL https://link.aps.org/doi/10.1103/RevModPhys.62.251
1990 doi
-
[11]
Hulse and Vivek Jayaraman
Brad K. Hulse and Vivek Jayaraman. Mechanisms underlying the neural computation of head direction. Annual Review of Neuroscience, 43 0 (1): 0 31--54, 2020. doi:10.1146/annurev-neuro-072116-031516. URL https://doi.org/10.1146/annurev-neuro-072116-031516. PMID: 31874068
2020 doi
-
[12]
H.A. Kramers. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica, 7 0 (4): 0 284--304, 1940. ISSN 0031-8914. doi:https://doi.org/10.1016/S0031-8914(40)90098-2. URL https://www.sciencedirect.com/science/article/pii/S0031891440900982
1940 doi
-
[13]
Basnak, Rachel I
Anna Kutschireiter, Melanie A. Basnak, Rachel I. Wilson, and Jan Drugowitsch. Bayesian inference in ring attractor networks. Proceedings of the National Academy of Sciences, 120 0 (9): 0 e2210622120, 2023. doi:10.1073/pnas.2210622120. URL https://www.pnas.org/doi/abs/10.1073/p...
2023 doi
-
[14]
Beck, Peter E
Wei Ji Ma, Jeffrey M. Beck, Peter E. Latham, and Alexandre Pouget. Bayesian inference with probabilistic population codes. Nature Neuroscience, 9 0 (11): 0 1432--1438, Nov 2006. ISSN 1546-1726. doi:10.1038/nn1790. URL https://doi.org/10.1038/nn1790
2006 doi
-
[15]
Zavatone-Veth
Paul Masset, Shanshan Qin, and Jacob A. Zavatone-Veth. Drifting neuronal representations: Bug or feature? Biological Cybernetics, 116 0 (3): 0 253--266, Jun 2022. ISSN 1432-0770. doi:10.1007/s00422-021-00916-3. URL https://doi.org/10.1007/s00422-021-00916-3
2022 doi
-
[16]
Carsen Stringer, Marius Pachitariu, Nicholas Steinmetz, Matteo Carandini, and Kenneth D. Harris. High-dimensional geometry of population responses in visual cortex. Nature, 571 0 (7765): 0 361--365, Jul 2019. ISSN 1476-4687. doi:10.1038/s41586-019-1346-5. URL https://doi.org/1...
2019 doi
-
[17]
van Vreeswijk and H
C. van Vreeswijk and H. Sompolinsky. Chaos in neuronal networks with balanced excitatory and inhibitory activity. Science, 274 0 (5293): 0 1724--1726, 1996. doi:10.1126/science.274.5293.1724. URL https://www.science.org/doi/abs/10.1126/science.274.5293.1724
1996
-
[18]
Vastola, Zach Cohen, and Jan Drugowitsch
John J. Vastola, Zach Cohen, and Jan Drugowitsch. Is the information geometry of probabilistic population codes learnable? In Sophia Sanborn, Christian Shewmake, Simone Azeglio, Arianna Di Bernardo, and Nina Miolane, editors, Proceedings of the 1st NeurIPS Workshop on Symmetry...
2023
-
[19]
Viazovska
Maryna S. Viazovska. The sphere packing problem in dimension 8. Annals of Mathematics, 185 0 (3): 0 991--1015, 2017. ISSN 0003486X. URL http://www.jstor.org/stable/26395747
2017
-
[20]
Waterhouse
William C. Waterhouse. Do symmetric problems have symmetric solutions? The American Mathematical Monthly, 90 0 (6): 0 378--387, 1983. ISSN 00029890, 19300972. URL http://www.jstor.org/stable/2975573
1983
-
[21]
Sphere packings
Chuanming Zong. Sphere packings. Springer Science & Business Media, 2008
2008
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.