{"id":"a812bda5-ebf8-4399-a9c5-2b2c2af3909b","arxiv_id":"2507.23030","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Apparent power-law scaling under the PRG arises from heterogeneity of independent place fields, not from critical dynamics.","lead":"This paper finds that a standard coarse-graining analysis, the phenomenological renormalization group (PRG), can produce power-law scaling from independent place cells with heterogeneous spatial fields, with no critical brain dynamics required. The result challenges earlier claims that such scaling is evidence for criticality in the hippocampus.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's explanatory claim requires that all structured correlations in the data are spatial; residual non-spatial shared variability is untested and could make the exponent match coincidental.","rationale":"The reader's CONDITIONAL verdict already identifies the same weakest spot, and I agree. The paper's strongest evidence is the match between independent place-field simulations and data across four PRG exponents, plus the subsampling dependencies in Fig. 4. For that match to support the causal/explanatory claim, the model must be a faithful null for the joint activity, not merely for single-cell field statistics. The proposed residual-correlation check would settle this directly: it tests conditional independence given position, the exact assumption that separates the model's 'non-interacting' explanation from alternatives. If the check fails, the title's claim is not established, although the broader anti-criticality point would remain. I therefore do not change the verdict; the paper should be CONDITIONAL pending this test, along with error bars and code/data availability, which are secondary but relevant.","tokens_in":9781,"tokens_out":24436,"duration_ms":308830,"concrete_test":"Using the recorded position trace, estimate each neuron's position-conditional activity probability from its binarized trace and form residuals R_i(t) = σ_i(t) − P(active | x(t)). Compute the cross-neuron correlation of these residuals (for example, the largest eigenvalue of the residual correlation matrix, or the mean absolute off-diagonal correlation). Perform the same computation on simulations from the fitted Gamma-Poisson model with matched session length and neuron count. If the data residuals show significantly larger shared variance than the model residuals, the data contain non-spatial shared variability that the model does not capture, and the exponent match in Fig. 2F cannot be attributed to place-field heterogeneity. If the residual correlation statistics are statistically indistinguishable, the conditional-independence assumption is supported and the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central explanatory claim — that place-cell heterogeneity, not critical interactions, underlies the PRG power laws measured in CA1 — rests on the Gamma-Poisson model of Fig. 2 being the correct statistical null for the recorded population. That model is a position-conditioned independent Poisson process with a per-cell constant background rate: given the animal's position, every neuron is independent. Real hippocampal recordings during navigation are subject to global shared modulations (running speed, arousal, theta, sharp-wave ripples, reward-related transients) that are not explained by position. If such non-spatial correlations are present in the binarized traces, the greedy correlation-based clustering of the PRG could be driven by them, and the agreement between model and data exponents in Fig. 2F could occur for reasons other than place-field overlap. The random-PRG control in Fig. 1J only establishes that correlation-based pairing matters; it does not identify which correlations are responsible. Thus the load-bearing assumption is untested: the data may violate conditional independence given position. Note that the weaker logical point — that a non-interacting model can generate PRG power laws, so PRG scaling alone is not evidence for criticality — survives even if this assumption fails; the concern specifically targets the stronger 'underlies' claim in the title and abstract.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the phenomenological renormalization group (PRG) of Meshulam et al. to calcium imaging recordings of CA1 place cells from mice navigating a 40 m virtual linear track, reproduces the previously reported power-law scaling, and shows that a random-pair control yields trivial exponents. It then builds a statistical model of place cells as independent Poisson processes with Gamma-distributed field counts and position-conditioned firing rates, and shows that PRG exponents obtained from simulated independent neurons match the empirical exponents. A simplified one-dimensional tile model is solved analytically in the homogeneous case, demonstrating that no exact power law emerges and that two scaling regimes meet at a cluster size equal to the number of simultaneously active units; adding heterogeneity in field sizes produces approximate power laws for all four PRG observables. The model predicts a systematic dependence of the exponents on population size and environment size, and subsampling the data confirms these trends. The authors conclude that PRG power-law scaling in hippocampal data does not constitute evidence for critical dynamics, because population heterogeneity alone can generate it.","tokens_in":10101,"tokens_out":5740,"duration_ms":71119,"significance":"If the conclusions hold, the paper has a clear conceptual payoff: it identifies a concrete non-critical mechanism—heterogeneous spatial tuning with independent neurons—for the PRG power laws observed in hippocampus, and it offers testable predictions (dependence on population size and environment size) that distinguish this account from scale-invariant critical dynamics. The homogeneous tile model is derived analytically rather than fitted, and the exponents in the detailed model are not free parameters fitted to the PRG outputs; they emerge from forward simulations with parameters calibrated to field counts, transient rates, and transient durations. The random-PRG control and the subsampling analyses strengthen the empirical side of the paper. The main risk is the untested conditional-independence assumption, which is load-bearing for the stronger \"underlies\" interpretation in the title and abstract, though not for the weaker claim that PRG scaling alone is not evidence for criticality.","major_comments":[{"comment":"The explanatory claim that place-cell heterogeneity, rather than critical interactions, underlies the measured PRG power laws rests on the Gamma-Poisson independent place-field model being an accurate statistical null for the recorded population. The model assumes that, given position, every neuron is independent; real CA1 recordings are subject to global non-spatial modulations (running speed, theta, sharp-wave ripples, arousal, reward-related transients) that are not explained by position. If such non-spatial correlations are present in the binarized traces, the greedy correlation-based clustering of the PRG could be driven by them, and the agreement between model and data exponents in Fig. 2F could arise for reasons other than place-field overlap. The random-PRG control in Fig. 1J shows that correlation-based pairing matters, but it does not identify which correlations are responsible. Please add a control that incorporates a common latent gain or arousal process in the model, or an analysis showing that PRG pairing in the data is explained by spatial tuning rather than by non-spatial co-fluctuations; otherwise the title's \"underlies\" claim is not supported.","section":"Results, Fig. 2"},{"comment":"The PRG exponents are reported as point estimates without error bars or confidence intervals, and the text does not state the number of K points used for each fit or the fit quality (e.g., R² or residuals). This makes the comparisons across mice, between data and model, and across subsampled conditions difficult to weigh. For example, the mouse-to-mouse spread in Fig. 1I and the model-versus-data agreement in Fig. 2F are presented visually, but there is no statistical measure of whether the differences are significant. Please provide per-session bootstrap errors for alpha, beta, z, and mu, and report fit quality for each power-law fit, at least in the supplement.","section":"Results, Figs. 1I, 2F, and 4"},{"comment":"The heterogeneous tile model's smoothing effect is asserted as \"likely due to the heterogeneous length scales across neurons\" but is only demonstrated by simulation, with no quantitative criterion for \"more closely matched to a power law.\" The central explanatory role of the simplified model requires a quantitative comparison: please report log-log residuals, R² values, or an information-theoretic comparison between the homogeneous and heterogeneous fits, and, if possible, provide an analytic argument for why heterogeneity removes the Nactive transition. Without this, the simplified model remains suggestive rather than established.","section":"Results, Fig. 3, and End Matter"},{"comment":"The subsampling validation is not fully specified. It is unclear how \"environment size fraction\" is implemented in the data: is the track truncated to smaller spatial segments, or is the environment rescaled? Truncation may introduce boundary effects and alter running speed, reward structure, and lap structure, which the constant-speed model does not capture. The same ambiguity applies to the population-size subsampling, where the model curves and data curves may not use the same subsampling procedure. Please clarify the exact protocol and, if truncation is used, apply the same truncation protocol to the model before comparing with the data.","section":"Results, Fig. 4C-D"}],"minor_comments":[{"comment":"There is a typo in the abstract: \"heteregeneous\" should be \"heterogeneous.\"","section":"Abstract"},{"comment":"The Table I caption contains a typo: \"repectively\" should be \"respectively.\"","section":"End Matter, Table I"},{"comment":"The condition \"For K ≥ ρ\" uses an undefined symbol; it should be \"K ≥ Nactive\" consistently with the preceding text. Please define Nactive before its first use in Eq. (1).","section":"End Matter, Eqs. (1)-(2)"},{"comment":"The Fig. 3A caption refers to \"K = ρ,\" but ρ is not defined in the main text or figure; please use Nactive for consistency with the End Matter.","section":"Fig. 3 caption"},{"comment":"The statement that the random-PRG control gives \"trivial scaling exponents for α ≈ 1 and β ≈ 1, as expected from combining uncorrelated Gaussian processes\" is not derived; please provide a short justification or a citation.","section":"Results, Fig. 1J"},{"comment":"Reference [25] is an unpublished conference abstract; if it is used to support a claim, please provide more details or remove it.","section":"References"},{"comment":"The Discussion notes that the tile model uses heterogeneity in field size while the biological model uses heterogeneity in field number; a formal connection between these two forms of heterogeneity would strengthen the simplified model's relevance to the detailed simulations.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and makes a genuinely useful contribution by showing that a non-interacting, spatially heterogeneous model reproduces PRG power laws and that the exponents depend on recording size and environment size. The main issue is that the title and abstract claim that heterogeneity 'underlies' the observed scaling, but the key assumption of conditional independence given position is untested. If the authors add the recommended control for non-spatial shared variability, or soften the explanatory claim accordingly, I would be happy to see the paper published. The lack of error bars and the purely simulated heterogeneous tile model are also fixable and should be addressed in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Useful paper, and one I expect will change how people read PRG results in hippocampus. The core result is genuinely new: a non-interacting model with heterogeneous place fields reproduces the PRG exponents from CA1 data, and an analytically solvable homogeneous tile model shows why homogeneity does not produce a true power law. That analytic derivation is the strongest part; it gives an explicit mechanism—two length scales producing a kink in scaling at K = N_active—and the heterogeneous tile model shows smoothing of that kink. I think the paper is right on the main logical point: PRG scaling alone is not evidence for critical dynamics. The subsampling analysis in Fig. 4, where model and data agree on how exponents shift with population and environment size, is persuasive as far as it goes.\n\nThe soft spots are real but mostly addressable. The title says heterogeneity 'underlies' the power laws, but the model only demonstrates sufficiency, not necessity. The model builds in conditional independence given position; real CA1 activity during navigation is surely affected by running speed, arousal, theta, reward, and other shared modulations not explained by position. If those non-spatial correlations are present, the greedy correlation-based clustering could pair units for reasons unrelated to place-field overlap, and the exponent match could be partly coincidental. The random-PRG control shows correlation-based pairing matters, but it does not identify which correlations. I would want the paper to test or at least discuss this, e.g., by including a shared speed or arousal signal in the model, or by conditioning on behavior. That said, the weaker claim—that PRG scaling is not evidence for criticality—survives even if the model is not the exact null.\n\nMinor issues: exponents are reported without error bars (only per-mouse markers), and no code or data are provided, which makes it hard to check the fits and the calibration. The heterogeneous smoothing mechanism is described as 'likely' and is only simulated; an analytic handle would strengthen it, but it is not necessary. The comparison with Meshulam 2019 is speculative, as the authors acknowledge.\n\nWho is this for? Anyone using PRG on neural data, and anyone working on criticality in biological systems. It deserves a serious referee. My recommendation: engage with it, but push for a test of conditional independence and for release of the code. With that, the 'underlies' claim can be evaluated properly.","headline":"Confirms that PRG power laws do not imply criticality, but the stronger 'underlies' claim needs a test of conditional independence before it fully lands.","tokens_in":10555,"tokens_out":2571,"would_cite":true,"duration_ms":33038,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C20","82B28"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that power-law scaling under a phenomenological renormalization group applied to hippocampal data arises from heterogeneous, non-interacting place cells, not from critical dynamics.","keywords":["place cells","power laws","criticality","phenomenological renormalization group","heterogeneity","hippocampus","coarse-graining","spatial coding"],"falsifier":"A concrete test would be to compute PRG exponents from hippocampal data and from surrogate spike trains in which each neuron's spikes are shuffled within small spatial bins, preserving the spatial rate map, field counts, and field sizes while destroying any non-spatial correlations; if the surrogate exponents deviate substantially from the real ones, the model misses an essential interaction, whereas if they match, spatial heterogeneity alone suffices.","tokens_in":9608,"feed_emoji":"🧠","tokens_out":3121,"duration_ms":42062,"temperature":0.7,"pith_summary":"The paper argues that the power-law scaling seen when hippocampal recordings are coarse-grained with a phenomenological renormalization group (PRG) is a byproduct of heterogeneity in place fields, not a sign of critical dynamics. It reproduces the experimental PRG exponents with a model of independent place cells whose field counts and rates come from a Gamma-Poisson distribution, and it shows analytically that a homogeneous tile model does not produce exact power laws because a transition occurs at the number of active fields per location. Adding heterogeneity in field sizes smooths that transition, producing apparent scale invariance. The authors conclude that power-law scaling of hippocampal data under PRG does not constitute evidence for critical dynamics.","feed_headline":"Place-cell variety, not criticality, yields brain scaling laws","feed_subtitle":"Independent, heterogeneous place cells reproduce the hippocampal power laws previously tied to critical dynamics.","key_machinery":"The load-bearing object is the Gamma-Poisson place-field model: each neuron forms fields as a Poisson process with a Gamma-distributed rate, giving heterogeneous field numbers and sizes fitted to recordings, and neurons are independent conditional on spatial location. The argument also relies on a simplified binary tile model whose coarse-graining can be analyzed exactly: a homogeneous tile model has a scaling transition at $K^* = N_{\\rm active} = l/\\delta x$, the number of active fields at a location, and heterogeneity in field sizes smooths this transition into apparent power laws. This transition structure is what turns unit-level variability into apparent scale invariance.","core_discovery":"The central claim is that population heterogeneity alone can produce power-law scaling under the PRG when units are independent, so the PRG power laws observed in hippocampal CA1 activity do not indicate critical interactions. The experimental exponents are matched by a statistical model in which each neuron has independent, spatially structured place fields drawn from a Gamma-Poisson distribution, with correlations arising only from shared spatial location. A simplified binary tile model shows that a homogeneous population of place fields has non-power-law scaling with two regimes separated by $N_{\\rm active}$, the number of units active at a location, while heterogeneity in field size removes the sharp transition and yields apparent power laws. The resulting exponents depend systematically on population size and environment size, and the paper confirms these dependencies by subsampling the recorded data.","pith_inferences":["The same heterogeneity mechanism may apply to PRG analyses of other brain areas where neurons are tuned to continuous variables, such as orientation or head direction, even though the paper does not test those systems. ","A practical testable extension would be to run the PRG on surrogate spike trains that preserve each neuron's spatial rate map but destroy all non-spatial correlations; the paper's account predicts unchanged exponents. ","The paper implicitly suggests that null models preserving empirically measured heterogeneity should become standard controls before any biological coarse-graining result is attributed to criticality. "],"forward_implications":["PRG power-law exponents previously read as signatures of criticality in hippocampal data should be reinterpreted as consequences of heterogeneous spatial tuning. ","PRG exponents depend on population size and environment size, so comparisons between studies must control for these experimental parameters. ","Subsampling predictions from the model match the recorded data, supporting the claim that the scaling is not scale-invariant in origin. ","The distribution of tuning curves is a more direct and intuitive quantity to analyse than PRG exponents when interpreting neural population activity. "],"supporting_citations":[{"why":"Reports coarse-graining and hints of scaling in a large neuronal population, providing the earlier evidence of power-law scaling under PRG.","marker":"[5]"},{"why":"Reports coarse graining, fixed points, and scaling in a large population of neurons, the main previous claim linking PRG exponents to criticality.","marker":"[6]"},{"why":"Supplies a previous latent-variable model that failed to reproduce the PRG scaling, motivating the need for the heterogeneity mechanism.","marker":"[16]"},{"why":"Provides the large-environment data and the Gamma-Poisson statistical structure that the place-field model is built on.","marker":"[21]"},{"why":"Supplies the recorded CA1 dataset, the fitted model parameters, and the experimental details used for the comparisons.","marker":"[22]"},{"why":"Provides the scale-invariant regimes for one-dimensional droplet growth and coalescence used to interpret the tile-model PRG behavior.","marker":"[24]"}],"fun_headline_variants":["Power laws in hippocampus come from cell variety, not criticality","Heterogeneous place cells mimic critical brain scaling","Brain power laws explained by place-cell diversity, not criticality","Place-cell heterogeneity reproduces hippocampal power laws","No criticality needed: place-cell variety yields power laws"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fitted Gamma-Poisson place-field model captures the correlation structure of the recorded CA1 population, so that matching PRG exponents with independent neurons is not a coincidence; if unmodeled non-spatial correlations or interactions are present, the exponents could match for other reasons.","fun_headline_variants_meta":{"raw":{"variants":["Power laws in hippocampus come from cell variety, not criticality","Heterogeneous place cells mimic critical brain scaling","Brain power laws explained by place-cell diversity, not criticality","Place-cell heterogeneity reproduces hippocampal power laws","No criticality needed: place-cell variety yields power laws"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000148,"raw_usage":{"total_tokens":1128,"prompt_tokens":825,"completion_tokens":303,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":441,"completion_tokens_details":{"reasoning_tokens":225}},"tokens_in":441,"tokens_out":303,"duration_ms":3707,"temperature":1.0,"reasoning_tokens":225,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:06:28.980511+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test would be to compute PRG exponents from hippocampal data and from surrogate spike trains in which each neuron's spikes are shuffled within small spatial bins, preserving the spatial rate map, field counts, and field sizes while destroying any non-spatial correlations; if the surrogate exponents deviate substantially from the real ones, the model misses an essential interaction, whereas if they match, spatial heterogeneity alone suffices.","supporting_citations":[{"cited_title":"Meshulam, J","cited_arxiv_id":null,"evidence_quote":"Reports coarse graining, fixed points, and scaling in a large population of neurons, the main previous claim linking PRG exponents to criticality."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies a previous latent-variable model that failed to reproduce the PRG scaling, motivating the need for the heterogeneity mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the large-environment data and the Gamma-Poisson statistical structure that the place-field model is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the recorded CA1 dataset, the fitted model parameters, and the experimental details used for the comparisons."},{"cited_title":"Derrida, C","cited_arxiv_id":null,"evidence_quote":"Provides the scale-invariant regimes for one-dimensional droplet growth and coalescence used to interpret the tile-model PRG behavior."}],"review_version":1}