The Roads and the Traffic: A Neural Test of the Domain Boundary

What happens when you measure the wrong kind of correlation

The Trust Attractor carries a quantitative prediction. The supporting research names it the Dissipative Coordination Principle (DCP): coordination range may grow with energy throughput in systems that coordinate by invitation. We went looking for that signature in the living brain. The surviving spatial null clarified the limits of the measurement; the proposed dynamical signature still requires a validated estimator.

The DCP proposes a power-law relation ξ ~ Φ^ν between a defined coordination length and sustaining energy flux. What we can measure here is how quickly the similarity between two patches of cortex falls off with distance. That is not the length of the wiring, and it is not yet the dynamical reach the DCP is about; measuring that would need a method of its own.

The exponent ν (nu) is the proposed diagnostic. The DCP predicts positive ν where added energy throughput sustains adaptive, longer-range coordination. A positive ν does not settle the matter on its own: fixed wiring can show one too, and an adaptive system can spend extra energy on noise or local activity. Reading the sign needs like-for-like measures and a tested mechanism.

The brain looked like the ideal testing ground. Neural circuits have Hebbian plasticity (named for psychologist Donald Hebb), summarized as “cells that fire together wire together.” Connections strengthen when neurons coordinate and weaken when they do not. The network rewires itself in response to its own activity. Anesthetics provide a clean experimental knob: different drugs at different doses shift cortical metabolism by amounts the literature has estimated, while leaving tissue physically intact.

The neuroscientist Davor Curic and colleagues at the University of Calgary had the dataset we needed. Their 2024 Nature Communications paper asked whether anesthesia pushes the brain away from its critical state, the boundary between order and chaos where information processing is richest. They documented multiple transitions in mice using widefield calcium imaging that captures a wide expanse of dorsal cortex (a 9.5 × 9.5 mm field of view) at fifty frames per second. They shared the data: fifty-four recordings across nine conditions, three mechanistically distinct anesthetics (isoflurane, ketamine, pentobarbital) at multiple doses.

What We Measured

Calcium imaging relies on a fluorescent indicator inside the neurons that glows when the cell fires. The glow is slower than the firing that causes it, smearing each spike into a lingering flare. Deconvolving the raw fluorescence signals works backward from the smear to the spike, recovering the timing of the underlying neural activity and producing a grid of 4,666 cortical pixels, each pooling many neurons and reporting when that patch fired.

From these traces we generated binary spike matrices (a fired-or-not record for every pixel at every frame) using Curic’s recommended method: a thresholded derivative that marks the onset of each calcium transient. We computed pairwise Pearson correlations (a standard measure of how similarly two signals behave) between all pixel pairs and binned them by physical distance. We then fitted exponential decay curves to extract the spatial correlation length ξ for each recording.

Pick any two points on the cortical surface and ask how correlated their activity is. Nearby points will be more correlated than distant ones. The distance over which that correlation decays is the spatial reach of coordinated activity.

We also computed the susceptibility χ (chi), which integrates total excess correlation above the background floor. Think of ξ as how far a rumor can travel, and χ as how many people end up repeating it. Susceptibility can diverge at a continuous critical point in an appropriate thermodynamic limit; finite recordings retain finite values. It captures both the reach and the amplitude of correlation.

Two variants of χ appear below, distinguished by the geometry of the sum. χ_1D adds up the excess correlation along the distance axis alone, one contribution per separation. χ_2D adds it up over the cortical sheet itself, weighting each separation by how many pixel pairs actually sit that far apart, which gives the distant pairs, vastly more numerous on a two-dimensional surface, proportionally more say. The integration supplies distance units: χ_1D has pixels and χ_2D has square pixels here, or length and area after physical calibration. Rescaling distance by a multiplies them by a and a² respectively. Within-series ratios remain meaningful under the same spatial calibration; absolute values also depend on integration limits and field of view. One note on provenance: susceptibility is our addition rather than Curic’s. The published paper reports the correlation length, the correlation floor, the amplitude, and the spike rate, and does not report a susceptibility at all, so both variants below come from our own integration of their correlation functions.¹¹⁹¹

Figure N1: Spatial correlation C(r) as a function of inter-pixel distance for each condition. Nearby pixels are more correlated than distant ones; ξ is the characteristic distance over which that correlation decays. The curves shift vertically (different correlation floors and amplitudes) while maintaining similar decay lengths.

The Null Result

Across fifty recordings (four of the fifty-four were flagged as outliers in Curic’s metadata and excluded), spanning metabolic rates from 30% to 110% of awake baseline, the spatial correlation length did not vary. The metabolic rate Φ assigned to each condition is a calibrated estimate drawn from the anesthesia literature. For isoflurane, the values come from cerebral-metabolism measurements in mice (Wei et al., 2024) and humans (Alkire et al., 1997). For pentobarbital, they come from the barbiturate review of Slupe and Kirsch (2018). The two ketamine conditions sit apart. The lower dose (10 mg/kg) reaches the 110% upper end, because ketamine can raise cortical metabolism rather than suppress it (Langsjö et al., 2005); the higher dose (100 mg/kg) is assigned 75%. That non-monotonic metabolic response makes them unsuitable for a clean power-law fit, so they are omitted from the table below while remaining in the full nine-condition comparison that follows. The ninth condition, the pre-injection baseline of the paired pentobarbital series (Φ = 1.00, ξ = 22.0 ± 1.8), is also left out of the table; it returns in the paired analysis below.

Power-law fit: ν = −0.0002, R² < 0.001, p = 0.999. Bootstrap 95% confidence interval: [−0.09, +0.10]. ξ does not track metabolic rate.

Figure N2: The flatline: ξ does not track metabolic rate. Small markers are individual recordings; large edged markers are condition means with standard-error bars. Hue and shape indicate drug class; shade indicates dose. The dashed line marks the grand mean, ξ = 24.0 pixels; the fitted power law is flat (ν = −0.0002, p = 0.999). Isoflurane 2% (burst suppression) is the one condition that stands apart.

Susceptibility told the same story. χ_1D: ν = −0.04, p = 0.81. χ_2D: ν = −0.09, p = 0.71. Neither susceptibility variant scales with metabolic rate.

Figure N3: Susceptibility χ_2D versus metabolic rate. Like ξ, the integrated excess correlation shows no systematic scaling with energy throughput.

One condition stands out: isoflurane at 2% produces ξ = 34.7, about 46% longer than the baseline. It is the most informative point in the entire analysis, and we return to it below.

What Does Vary

A Kruskal-Wallis test (comparing groups without assuming a particular distribution shape) across all nine conditions reveals which observables respond to anesthesia and which do not:

The nine-condition correlation-length comparison did not reject its null (Kruskal-Wallis H = 13.2, p about 0.105). Non-rejection does not establish equal distributions or practical invariance. Susceptibility, correlation floor, amplitude, and spike rate show reported condition differences. A prespecified equivalence margin and corresponding interval or test would be needed to establish that spatial reach stays within a meaningful tolerance.

Figure N4: The reported condition comparison for spatial activity-correlation length does not reject equality at the 5-percent level. This is not an equivalence test or evidence of fixed anatomical reach.

Think of a road network. The roads have a characteristic length: how far you can drive before hitting the edge of town. In the cortex, the roads are the anatomical wiring, and the traffic is the activity that flows along it. Anesthesia changes the traffic (how much flows, how fast, whether it clumps or spreads) while leaving the roads untouched: it can change activity correlations without changing anatomical wiring. The observed spatial decay is therefore closer to a traffic-pattern scale than a direct measure of road length; the anatomical explanation remains unestablished.

Confirmation, One Mouse at a Time

Curic’s design included a paired component. Four mice were each recorded at baseline (pre-injection), thirty minutes after pentobarbital 80 mg/kg, and sixty minutes after. Each animal was measured under three metabolic states, eliminating between-subject variability. (The Friedman test below is the paired counterpart of the earlier Kruskal-Wallis comparison: the same question, asked within each animal rather than across groups.)

The paired analysis confirms the pooled result at the individual-animal level. In all four mice, pentobarbital drops the correlation floor (C_∞ decreases, Friedman p = 0.039) and raises correlation amplitude (Friedman p = 0.018). The drug peels away a background haze of global synchrony, leaving local correlation sharper, as fog clearing reveals the contours of a landscape. Spike rate rises in all four animals, the increase near the threshold of significance (Friedman p = 0.050).

The correlation length ξ: two mice up, two mice down. Friedman p = 0.78. Even within individual animals undergoing a threefold metabolic suppression, the spatial reach shows no consistent shift.

The Burst-Suppressed Cortex

Isoflurane at 2% is the exception that tests the rule. At this dose, the drug forces the cortex into burst suppression: synchronous oscillations where vast swaths of tissue fire in lockstep, followed by silent periods where nothing fires. The cortex has been pharmacologically seized. Think of a command economy where every factory produces the same product on the same schedule. The central authority has removed every alternative.

The numbers tell the story. Isoflurane 2% has the longest correlation length in the dataset (ξ = 34.7) and the highest correlation floor (C_∞ = 0.37, meaning distant pixel pairs retain a Pearson correlation of 0.37 regardless of separation).

Within the isoflurane dose series, the 2D susceptibility χ_2D increases as metabolic rate decreases: 1,821 → 2,205 → 2,970, a rise of 63% from the awake cortex to the burst-suppressed one. The power-law fit gives ν = −0.53 with R² = 0.96, though it rests on only three condition means (p = 0.12) and describes a trend rather than an established law. The negative sign means more suppression produces more apparent coordination. On the DCP’s reading, a negative ν is the signature of imposed coordination, the same sign the driven six-spin ring with a fixed Hamiltonian gives (−0.40), and here a drug supplies it on demand: the order in a burst-suppressed cortex is administered rather than grown.

Figure N5: The isoflurane dose series, from baseline (green) through 1% (light blue) to 2% (dark blue, burst suppression). All three curves decay with distance, and the 1% curve settles slightly below baseline at long range. At 2%, the burst-suppressed cortex holds the strongest long-range correlation, with the highest floor and the longest fitted decay length.

The resulting “coordination” is long-range, high-amplitude, and brittle. Burst suppression alternates between synchronous firing and silence: a state incompatible with sustained computation.

Compare the awake cortex, which sustains lower-amplitude, shorter-range coordination indefinitely. It adjusts its patterns moment by moment in response to sensory input, memory retrieval, and internal computation. Less striking in a snapshot; incomparably more capable over time.

Two Correlation Lengths in One Dataset

The dataset’s spatial activity correlations and avalanche distributions measure different aspects of dynamics. Their comparison leaves open how either maps to a dynamical correlation length.

The spatial correlation length is a fitted decay scale of activity correlation C(r). Anatomical connections constrain activity, but this observable is not tract length: coupling, driving, and noise can change correlations on fixed wiring. The reported dose-dependent ξ estimates therefore cannot be assigned directly to anatomical reach. A flat fitted relation with estimated metabolism leaves the causal explanation open.

Avalanche statistics describe cascade size and duration distributions. Their fitted slopes alone are not a dynamical correlation-length measurement; a cutoff or justified scaling relation would be needed.

One mechanism contributing to this traffic is ephaptic coupling (Chapter 17): each neuron’s firing generates an electric field that perturbs its neighbors without any synaptic connection. This may shape cascade propagation in real time. Activation patterns are either scale-free (near criticality, cascades of all sizes, long dynamical ξ) or truncated (away from criticality, only small cascades, short dynamical ξ). Cascade size is one rough indicator of reach: a cascade that dies after three pixels has coordinated three pixels, while one that crosses the cortex has coordinated the cortex.

Avalanche power-law exponents describe distribution slopes. A dynamical correlation length requires additional information, such as a cutoff, finite-size scaling, and a justified relation involving spatial or dynamic exponents. No such relation here derives the reported ν = +0.64 ± 0.07 from those slopes. That number is withdrawn from the dynamical-length comparison until the estimator and its assumptions are supplied.

The spatial power-law fit reports ν near zero against assigned metabolic rates. The avalanche analysis has not supplied a comparable length estimator, so a positive dynamical ν cannot currently be set beside that fit.

The Classification Sharpens

This result fits into the progression across scales. The quantum spin-chain exponents come from the author’s earlier six-spin driven-dissipative ring (exact Liouvillian), a separate run from the eight-spin chains in Chapter 17; the social values come from the World Values Survey trust regression (109 countries) and the firm-size analysis discussed there; the neural rows are the present analysis.

The reported exponents span distinct observables and estimators. Their signs do not by themselves identify whether a network’s anatomical backbone can reorganize. The proposed energy-to-adaptive-reach mechanism requires a validated dynamical-length measure and comparable flux calibration.

Behind the positive sign is a simple economy: every link in an adaptive network costs energy to keep, so more throughput can sustain more links, and longer ones. Friendships form when people find value in each other and fade when they do not; adaptive networks reorganize the same way, their links sustained only as long as the dynamics reward them.

The Curic dataset permits comparisons of spatial correlations and avalanche statistics in one physical system. It currently supplies no positive dynamical-length exponent or demonstrated within-brain boundary for the DCP.

What This Means for the Prediction

The ξ ~ Φ^ν prediction refers to the dynamical correlation length: the range over which activity patterns organize adaptively. The mechanisms involved can strengthen, weaken, form, or dissolve in response to energy flux, as trade routes open when commerce is profitable and close when it is not.

Adaptive topology could bring spatial and dynamical correlation lengths into closer correspondence: the roads reshape themselves to match the traffic. That correspondence depends on the interaction dynamics, adaptation timescale, and definitions of the two lengths. A fixed anatomical backbone may separate them, although fixed wiring alone does not require decoupling.

If the boundary holds, it sharpens the DCP: the principle would apply to systems whose coordination is maintained by invitation, where every link exists because the dynamics sustain it. Removing energy flux would then cause coordination to dissolve rather than merely fall silent.

Imposed Order Looks Stronger in a Snapshot

Imposed coordination looks more powerful than emergent coordination in a snapshot. The burst-suppressed cortex under isoflurane 2% produces the longest correlation length in the dataset. Measured by the range of synchronized activity at a single moment, coercion wins.

That impression dissolves when you ask what the coordinated system can do. The burst-suppressed cortex cannot process information, respond selectively to what it senses, form memories, generate predictions, or sustain computation. Every functional capacity has been sacrificed for the appearance of coordination.

The correlation length is one measurable axis of a difference that involves every aspect of neural processing.

The social parallel is an analogy, not a measurement. An authoritarian state can mobilize its entire population for a single purpose, producing impressive snapshots of coordinated action. A democracy, with its shorter “correlation length” of consensus, sustains adaptive governance over decades.

The Trust Attractor predicts that invitation-based governance recovers more readily from perturbation. Burst suppression offers a physiological analogy: synchronous firing alternates with silence, while the awake cortex sustains varied activity. Establishing greater thermodynamic stability in democratic governance requires a social dynamical model and its own evidence.

A snapshot measures reach. A trajectory measures resilience. The DCP is a statement about trajectories.

Data: Curic, D. et al. “Existence of multiple transitions of the critical state due to anesthetics.” Nat. Commun. 15, 7025 (2024). We thank Davor Curic for generously sharing the deconvolved calcium imaging data and for his expert guidance on spike detection from calcium indicators. Analysis scripts and results are available in the project repository.