hierarchical-critical-brain-dynamics - SKILL.md Agent Skill

name: hierarchical-critical-brain-dynamics description: "Hierarchical organization of critical brain dynamics. Analysis of how brain structure hierarchies interact with criticality hypothesis. Activation: hierarchical brain, critical dynamics, connectome hierarchy, brain criticality."

Hierarchical Organization of Critical Brain Dynamics

Investigating how the hierarchical organization of the brain interacts with the criticality hypothesis for collective neural dynamics.

Metadata

Source: arXiv:2604.21832v1
Authors: Gustavo G. Cambrainha, Daniel M. Castro, Leonardo L. Gollo, Pedro V. Carelli, Mauro Copelli
Published: 2026-04-23
Categories: q-bio.NC, physics.bio-ph, q-bio.QM
PDF: https://arxiv.org/pdf/2604.21832v1

Core Methodology

Research Question

How does the hierarchical organization of the brain (a fundamental structural principle) interact with brain criticality (a leading hypothesis for collective dynamics)? The study uses phenomenological renormalization group approaches applied to large-scale neuronal spiking activity from mouse visual cortex and hippocampus.

Key Findings

1. Hierarchy-Dependent Criticality

Signatures vary systematically along the known anatomical hierarchy
Measure-dependent organization: Static property exponents point in one direction, dynamic property exponents point in the opposite direction
Task modulation: Visual system signatures strongly modulated by engagement in visual tasks

2. Dynamic Reconstruction of Hierarchy

Correlations among criticality markers during active engagement can reconstruct anatomical hierarchy
Scaling exponents closely follow theoretically predicted relations
Exponents covary with hierarchical position

3. Multi-Region Analysis

Visual cortex: Clear hierarchical organization of critical signatures
Hippocampus: Similar hierarchical patterns despite different anatomical structure
Cross-regional consistency validates the framework

Analytical Framework

Phenomenological Renormalization Group (PRG)

The PRG approach coarse-grains neural activity across spatial scales to identify:

Scaling relations between different criticality markers
Hierarchy-dependent critical exponents
Task-modulated critical signatures

prg_analysis = {
    "static_exponents": ["tau", "alpha", "sigma"],  # Power-law exponents
    "dynamic_exponents": ["tau_t"],  # Temporal scaling
    "scaling_relations": {
        "tau": "avalanche size distribution",
        "alpha": "avalanche duration distribution", 
        "sigma": "size vs duration relation",
        "tau_t": "temporal correlation decay"
    },
    "hierarchy_gradient": "systematic variation along anatomical hierarchy"
}

Criticality Indicators by Type

criticality_markers = {
    "static_properties": {
        "avalanche_size": "Power-law fitting (tau)",
        "avalanche_duration": "Power-law fitting (alpha)",
        "size_duration_relation": "Exponent sigma"
    },
    "dynamic_properties": {
        "temporal_correlation": "Decay exponent tau_t",
        "branching_ratio": "sigma ≈ 1 (critical)",
        "susceptibility": "Divergence near critical point"
    },
    "task_modulation": "Engagement-dependent signature changes"
}

Implementation Guide

Connectome Analysis

import numpy as np
import networkx as nx
from scipy import stats

class HierarchicalCriticalityAnalyzer:
    """
    Analyze hierarchical organization and critical dynamics
    """
    
    def __init__(self, connectivity_matrix, node_hierarchy):
        self.adj = connectivity_matrix
        self.hierarchy = node_hierarchy
        self.graph = nx.from_numpy_array(connectivity_matrix)
        
    def compute_hierarchy_metrics(self):
        """
        Quantify hierarchical organization
        """
        # Trophic levels
        trophic = self._trophic_levels()
        
        # Hierarchical clustering
        clustering = self._hierarchical_clustering()
        
        # Fractal analysis
        fractal_dim = self._box_counting_dimension()
        
        return {
            "trophic_coherence": np.std(trophic),
            "hierarchy_height": max(trophic) - min(trophic),
            "fractal_dimension": fractal_dim,
            "modularity": nx.algorithms.community.modularity(
                self.graph, 
                nx.community.greedy_modularity_communities(self.graph)
            )
        }
    
    def analyze_avalanche_dynamics(self, spike_data, threshold):
        """
        Detect neural avalanches and test for criticality
        """
        # Detect avalanches
        avalanches = self._detect_avalanches(spike_data, threshold)
        
        # Size distribution
        sizes = [len(a) for a in avalanches]
        
        # Power-law fitting
        fit = self._powerlaw_fit(sizes)
        
        # Criticality tests
        branching = self._estimate_branching_ratio(avalanches)
        
        return {
            "tau": fit['exponent'],  # Power-law exponent
            "p_value": fit['p_value'],
            "branching_ratio": branching,
            "is_critical": 0.9 < branching < 1.1 and fit['p_value'] > 0.1
        }
    
    def cross_scale_analysis(self, spike_data, resolutions):
        """
        Analyze criticality across spatial scales
        """
        results = {}
        
        for res in resolutions:
            # Coarse-grain at this resolution
            coarse_data = self._coarse_grain(spike_data, res)
            
            # Analyze criticality
            crit = self.analyze_avalanche_dynamics(coarse_data, threshold=1)
            
            results[f"scale_{res}"] = crit
        
        return results

Hierarchical Criticality Model

class HierarchicalIsingModel:
    """
    Hierarchical Ising model for brain dynamics
    """
    
    def __init__(self, hierarchy_depth, branching_factor):
        self.depth = hierarchy_depth
        self.branching = branching_factor
        self.build_hierarchy()
        
    def build_hierarchy(self):
        """
        Construct hierarchical connectivity
        """
        self.nodes = []
        for level in range(self.depth):
            n_nodes = self.branching ** level
            self.nodes.append({
                'level': level,
                'count': n_nodes,
                'coupling': 1.0 / (level + 1)  # Decreasing coupling with level
            })
    
    def simulate(self, temperature, timesteps):
        """
        Monte Carlo simulation
        """
        # Initialize spins
        spins = np.random.choice([-1, 1], size=self.total_nodes())
        
        # Metropolis dynamics
        for t in range(timesteps):
            for i in range(len(spins)):
                # Compute local field
                h = self.local_field(i, spins)
                
                # Metropolis update
                delta_E = 2 * spins[i] * h
                if delta_E < 0 or np.random.random() < np.exp(-delta_E / temperature):
                    spins[i] *= -1
        
        return spins
    
    def compute_susceptibility(self, temperatures):
        """
        Find critical temperature from susceptibility peak
        """
        susceptibilities = []
        
        for T in temperatures:
            spins = self.simulate(T, 10000)
            magnetization = np.mean(spins)
            variance = np.var(spins)
            susceptibilities.append(variance / T)
        
        # Critical temperature at peak
        T_c = temperatures[np.argmax(susceptibilities)]
        
        return T_c, susceptibilities

Key Insights

1. Hierarchical Structure Enables Criticality

Modularity creates locally stable dynamics
Inter-module connections enable global coordination
Hierarchy depth controls critical scaling

2. Criticality Optimizes Hierarchical Processing

Maximized dynamic range at each level
Efficient information routing between modules
Adaptability through critical fluctuations

3. Multi-Scale Information Processing

Lower levels: Fast, local computations
Higher levels: Slow, integrative processing
Critical dynamics facilitate cross-scale communication

Theoretical Implications

Brain Organization Principles

Structural Hierarchies: Evolutionarily conserved
Critical Dynamics: Generic property of complex networks
Reciprocal Optimization: Structure and dynamics co-evolve

Information Processing Benefits

Efficiency: Minimal wiring cost
Robustness: Graceful degradation
Adaptability: Context-dependent routing
Capacity: Maximized information storage

Applications

1. Brain Disease Modeling

Disruption of Hierarchy: Disconnection syndromes
Loss of Criticality: Epilepsy, anesthesia
Altered Scaling: Neurodegenerative diseases

2. Artificial Neural Networks

Architectural Design: Hierarchical organization
Training Dynamics: Critical initialization
Information Routing: Attention mechanisms

3. Brain-Computer Interfaces

Optimal Stimulation: Exploiting critical dynamics
Signal Decoding: Multi-scale analysis
Closed-Loop Control: Feedback at critical point

Related Skills

brain-network-controllability
neural-critical-dynamics-theory
griffiths-phase-brain-criticality
brain-criticality-hypothesis-assessment

References

Cambrainha, G.G. et al. (2026). Hierarchical organization of critical brain dynamics. arXiv:2604.21832.
Beggs, J.M. & Plenz, D. (2003). Neuronal avalanches in neocortical circuits.
Kaiser, M. (2007). Brain architecture: A design for natural computation.

Implementation Status

Theoretical framework
Connectome data analysis
Computational modeling
Clinical applications
BCI integration