Higher-order interactions renormalize the thermodynamic signature of neural population codes
Abstract
Maximum-entropy models of neural recordings can be treated as Boltzmann distributions at a fictitious temperature T, with the fitted model at T=1. Their specific-heat peak grows with population size and approaches that temperature, which has been read as a sign of criticality. The case rests on pairwise models, yet retinal codes require higher-order interactions. The thermodynamics of a model with such interactions fitted cell by cell has not been computed. Here I compute it exactly for pairwise and third-order models of the same populations of up to $20$ retinal ganglion cells, in salamander and marmoset. The triplet terms, required by the held-out likelihood, shift the peak by more than the finite-sample bias, with opposite sign in the two retinas. The peak is lower and grows more slowly in salamander, higher and grows faster in marmoset, and lies closer to T=1 in both. An exact decomposition traces most of the shift to the energy variance across population spike counts. A constraint on the spike-count distribution alone reproduces all of it in salamander and over half in marmoset. A second recording of each species, re-binning, and a common-input surrogate show that the sign is set by the spike-count structure of a recording at its temporal resolution, not by species. The thermodynamic signature of a neural code is therefore set not by the data alone but by the interaction order of the model used to read it.
Related
- Anisotropic magneto-elastic coupling in the quasi-2D Ising magnet Li2FeSiO4 and thermodynamic signature of low-energy orbital excitations
- Beyond the Binary Hull: Higher-Order Thermodynamic Stabilization in Inorganic Ternary Compounds
- High-throughput thermodynamic fingerprinting of protein–ligand interactions by DNA-directed focal molography
- Generalized Dunkl Quantum Systems with Energy-Dependent Interactions: Exact Solvability and Thermodynamic Properties
- Perturbations of Charged Black Holes with Higher-Order Interactions
BibTeX
@misc{hegade2026higherorder,
title = {{Higher-order interactions renormalize the thermodynamic signature of neural population codes}},
author = {Narendra N. Hegade},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.05336},
eprint = {2611.05336},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.05336},
note = {CC BY-SA 4.0},
}
RIS
TY - UNPB
TI - Higher-order interactions renormalize the thermodynamic signature of neural population codes
AU - Narendra N. Hegade
PY - 2026
DA - 2026/11/03
PB - arXiv
DO - 10.48550/arxiv.2611.05336
UR - https://arxiv.org/abs/2611.05336
AB - Maximum-entropy models of neural recordings can be treated as Boltzmann distributions at a fictitious temperature T, with the fitted model at T=1. Their specific-heat peak grows with population size and approaches that temperature, which has been read as a sign of criticality. The case rests on pairwise models, yet retinal codes require higher-order interactions. The thermodynamics of a model with such interactions fitted cell by cell has not been computed. Here I compute it exactly for pairwise and third-order models of the same populations of up to $20$ retinal ganglion cells, in salamander and marmoset. The triplet terms, required by the held-out likelihood, shift the peak by more than the finite-sample bias, with opposite sign in the two retinas. The peak is lower and grows more slowly in salamander, higher and grows faster in marmoset, and lies closer to T=1 in both. An exact decomposition traces most of the shift to the energy variance across population spike counts. A constraint on the spike-count distribution alone reproduces all of it in salamander and over half in marmoset. A second recording of each species, re-binning, and a common-input surrogate show that the sign is set by the spike-count structure of a recording at its temporal resolution, not by species. The thermodynamic signature of a neural code is therefore set not by the data alone but by the interaction order of the model used to read it.
ER -