Neurocomputing

Wave recurrent neural networks (wRNNs) are biologically inspired recurrent architectures that use traveling-wave dynamics to support sequence learning and memory. However, their input-to-hidden pathway remains relatively simple compared with biological neurons, where dendrites perform nonlinear input integration. In this study, we introduce the Dendritic Wave Recurrent Neural Network (DWRNN), which augments the input pathway of the wRNN with nonlinear basal dendritic branches while preserving the original recurrent wave dynamics. We evaluate DW-RNN on a simple copy task, sequential MNIST (sMNIST), permuted sequential MNIST (psMNIST), and noisy sequential CIFAR-10 (nsCIFAR-10). On the copy task, DW-RNN shows learning behavior comparable to the standard wRNN, suggesting that dendritic input integration does not disrupt the recurrent wave-based memory mechanism. On the three sequential image-classification benchmarks, DW-RNN outperforms the standard wRNN, improving accuracy from 97.27 {+/-} 0.15% to 97.82 {+/-} 0.12% on sMNIST, from 96.74 {+/-} 0.17% to 96.92 {+/-} 0.10% on psMNIST, and from 54.30 {+/-} 0.79% to 55.65 {+/-} 0.55% on nsCIFAR-10. In addition to improving mean accuracy, DW-RNN exhibits lower across-seed variability on all three classification benchmarks, suggesting that dendritic input integration may improve the stability of wRNN training. Hidden-activity visualizations further show that DW-RNN preserves the characteristic traveling-wave patterns of the original wRNN. These results suggest that dendritic computation and traveling-wave recurrent dynamics provide complementary mechanisms for biologically inspired sequence learning.

Abnormal states such as erroneous behaviors are generally difficult to represent from neural data. However, such states are also known to have specific spatiotemporal features, indicating a feasibility of developing a method to focus on them. If a method can highlight these spatiotemporal features, it may effectively represent such abnormal states, helping evaluate abnormal brain functions. In the present study, we proposed the hierarchy of supported modules (HSM) to highlight spatiotemporal features that can represent abnormal states. HSM spatiotemporally transforms multidimensional neural time-series based on their spatiotemporal context. We evaluated HSM through decoding and similarity analyses using multiple publicly available datasets. In the HSM results, decoding accuracies were higher for erroneous behaviors than for normal behaviors, and similarities were lower between erroneous behaviors and normal behaviors than between normal behaviors, demonstrating the ability of HSM to capture the spatiotemporal features of erroneous behaviors. Surprisingly, many parts of these results were also present even before HSM learning, showing the virtue of HSM as a simple-to-use method. The proposed HSM method may help elucidate the mechanisms underlying erroneous behaviors.

Transcranial Electric stimulation (tES) is a safe and noninvasive technique increasingly used in treating brain disorders. Despite many studies on tES, there is still a lack of understanding of its mechanisms at the microscale network level. This is crucial for optimized parameter selection in therapy approaches such as the treatment of pharmacoresistant epilepsy. In this study, we made use of a recently published neuroinspired microscale model of the neocortex, known as NeoCoMM, and integrated a "Lambda E"-based model of tES. This updated version was used to investigate the acute effects of tES (tDCS and tACS), on the neural activity of various neuron types in both healthy and epileptic brain states. Results showed that in the case of healthy alpha and gamma rhythms, tACS induced electric field entrainment at the peak power frequency of the network oscillations as measured by the Local field Potentials (LFPs). This resonance-like entrainment was independent from the individual firing rate of cell types. For epileptic activity, tACS did not provide consistent results. Cathodal tDCS resulted in a promising decrease in hyperexcitable activity throughout simulations. These results advance our understanding of the impact of tES on network dynamics at both the extracellular and intracellular activity levels. Author summary

In patients with epilepsy, seizures are associated with pathological neural synchronization. However, the preictal period preceding a seizure often exhibits reduced spatial synchronization compared to normal cognition. This observation aligns with the concept of the brain as a complex dynamical system, where a reduction in dimensionality and resilience can precede a phase transition. The Critical Brain Hypothesis suggests a connection between the loss of healthy scale-free behavior and various disorders, including epilepsy. Our study investigates preictal changes by utilizing network features, such as mean node degree and mean clustering coefficient, derived from thresholded correlation matrices of patient intracranial electrocorticographic electrode data. We observed a suppression of intermittent high-synchronization periods within the feature space during the minutes leading up to seizure onset. This constriction of the explored hypervolume in the preictal state indicates a breakdown in the brains ability to maintain normal coherence. We use these preictal changes to predict the probability of seizure onset using a Support Vector Machine algorithm. These discrete predictions can then be combined into real-time continuous seizure risk forecasts via Bayesian updating. This innovative and computationally lightweight approach has the potential to significantly improve upon static predictions, providing opportunities for more adaptable, quantitative, and interpretable tools for managing seizures.

Mapping hierarchical brain networks within traditional Euclidean space causes significant structural distortion, undermining neuroimaging diagnostic frameworks. While hyperbolic models like the Poincare ball preserve these nested topologies, they demand heavy computational overhead due to intricate Mobius operations and curved geodesics. This paper introduces a highly efficient non-Euclidean framework for analyzing neurocognitive decline utilizing the Beltrami-Klein ball model. By projecting hyperbolic geodesics as Euclidean straight lines, this approach converts complex distance calculations into simple dot products, radically reducing processing demands. We validated our methodology against state-of-the-art Poincare and Lorentz baselines using datasets for Schizophrenia, Parkinsons Disease, and Alzheimers Disease. The Klein-based framework demonstrates superior performance, delivering both higher diagnostic precision and accelerated processing velocities across all three neurocognitive disorders.

Background: Free-text EEG reports typically lack structure, hindering scalable analysis. We evaluate a large language model (LLM) pipeline to extract structured diagnostic labels and confidence levels from these reports. Methods: We developed a hierarchical annotation schema to classify EEG reports for four specific abnormality types using a four-point confidence scale. To establish ground truth, two certified EEG technicians annotated a diverse dataset of reports authored by neurologists with distinct writing styles. We then implemented a grammar-constrained Mistral-7B pipeline, iteratively prompt-tuned on a development set to mirror these expert annotations. The pipeline's effectiveness was evaluated against the human expert benchmark using core agreement (diagnostic accuracy) and certainty-adjusted agreement (confidence alignment), with classical NLP models serving as a secondary baseline. Results: Mistral-7B significantly outperformed baselines, achieving 96% accuracy for overall abnormality detection, approaching the human benchmark of 98%. Crucially, the model successfully identified rare epileptiform abnormalities where traditional models failed and generalized robustly across distinct reporting styles. While diagnostic accuracy was high, a performance gap persisted in certainty-adjusted agreement, indicating that accurately modeling nuanced clinical confidence remains a challenge. Conclusion: LLMs can effectively automate the extraction of structured diagnostic information from EEG reports with near-human accuracy and strong generalization. While confidence calibration requires further refinement, the combination of accurate classification and explainability makes this pipeline a promising tool for standardizing clinical data at scale. Keywords: Routine Clinical Electroencephalography; Large Language Models; Clinical NLP; Confidence Assessment; Explainable AI; Neurophysiological Evaluation

Transformer-based deep learning models have shown great potential for decoding visual EEG signals. However, their internal attention mechanisms are often evaluated primarily on optimization objectives, leaving their alignment with biological brain connectivity an open question. This study empirically evaluates how variations in EEG preprocessing strategies affect these attention representations using the Adaptive Thinking Mapper (ATM) model as a framework. We compared a baseline pipeline (MVNN only) against a comprehensive cleaning pipeline integrating ICA and notch filtering. The models were evaluated through cross-generalization, noise robustness, and spectral-temporal ablation analyses. Furthermore, we investigated the structural correspondence between the model's data-driven attention weights and neurophysiological reference networks (GPDC, PDC, and DTF) using Node Strength Correlation and Representational Similarity Analysis (RSA). The results show that the comprehensive preprocessing successfully suppresses non-neural artifacts, such as frontal noise and electrical interference, while maintaining comparable decoding accuracy and baseline robustness. Alignment analyses revealed that the broad spatial organization of the learned attention patterns remains highly stable across pipelines, capturing key directed connectivity dynamics with subtle, metric-dependent variations in global representational geometry. This work provides an empirical exploration into bridging data-driven attention weights with neurophysiological consistency, offering insights toward more transparent brain-computer interfaces.

The aim of this work is to validate patient-specific EEG baseline establishment using the e-norms method as a screening and retrospective-review tool for seizure detection in pediatric epilepsy. The method was applied to 247 seizure-free EEG recordings (263.92 hours) from 10 patients in the CHB-MIT Scalp EEG Database (ages 3-18). A composite stability metric combining first-derivative dynamics, spectral entropy, variance, and line length was computed per 2-second epoch across 23 channels. Patient-specific detection thresholds were derived from each patient's seizure-free baseline using a weighted statistical procedure. Performance was validated against 72 expert-annotated seizures (2,705 epochs) across 62 seizure files, with durations spanning 6 to 264 seconds (44-fold range). The results show that detection achieved 94.4% event-level sensitivity (68 of 72 seizures; 95% CI 86.6-97.8%) and 81.5% epoch-level sensitivity (2,204 of 2,705 epochs; 95% CI 80.0-82.9%). Eight of ten patients achieved 100% event-level sensitivity with epoch-level sensitivity ranging from 58.7% to 100.0%. Two patients showed partial event-level failures (CHB-15: 17 of 20; CHB-18: 5 of 6), with the four missed events attributable to two characterizable failure modes. Patient-specific thresholds ranged from 4.06 to 4.81 (mean 4.51 +/- 0.25); threshold variation did not correlate reliably with age or sex, confirming that no universal threshold could achieve comparable performance. Detection margins ranged from 0.88 to 1.24 times. Patient-specific e-norms achieves 94.4% event-level sensitivity for pediatric EEG seizure detection without requiring labeled seizure training data, exceeding published human expert inter-rater agreement (50-76%) and recent automated approaches in adult cohorts using behind-the-ear EEG and wearable ECG. Two characterizable failure modes account for the four missed events and inform appropriate clinical use. As a high-sensitivity screening tool complementary to real-time alarm systems, the method is ready for adult validation, prospective deployment, and head-to-head benchmarking.

Deep brain stimulation has demonstrated its therapeutic potential in modulating pathological oscillations associated with Parkinsons disease and epilepsy. However, its efficacy in treating disrupted theta-gamma phase-amplitude coupling seen in memory-related disorders, such as Alzheimers disease, remains poorly understood. While recent studies have targeted the entorhinal-hippocampal circuit, results remain inconsistent. This discrepancy stems from a lack of mechanistic understanding regarding how stimulation protocols affect this circuit. In this work, we present a reduced multicompartment model of the hippocampal CA1 area that reproduces theta-nested gamma oscillations characteristic of healthy neural activity during memory performance. The model comprises pyramidal, basket and OLM cells with simplified morphologies. We also incorporated CA3-to-CA1 axonal projections, providing a foundational framework for studying how stimulation-induced recruitment of afferent pathways modulates CA1 dynamics. By balancing computational efficiency with anatomical accuracy, our model enables systematic investigation of the effects of electrode placement and orientation, as well as stimulation amplitude and frequency on CA1 neural activity. We demonstrate that the excitatory response in CA1 is primarily driven by the recruitment of Schaffer collateral projections. Overall, this work provides a computationally efficient template for exploring diverse stimulation configurations and could be expanded for developing neuromodulatory strategies to restore physiological network dynamics. Author summaryDeep brain stimulation has shown success in treating Parkinsons disease by suppressing abnormal neural activity responsible for movement disorders. However, when applied to memory-related pathologies, such as Alzheimers disease, the therapeutic outcomes remain unpredictable, ranging from cognitive improvement to impairment. This discrepancy highlights a critical gap in our understanding of how stimulation protocols interact with neural dynamics of the targeted circuits. To address this, we developed a computationally efficient model of the hippocampus, which is involved in memory processes, in order to understand how deep brain stimulation might influence its activity. Our model maintains enough biological accuracy to capture essential memory-related neural activity while remaining lightweight enough for rapid execution and systematic exploration of different protocols. This computational efficiency allowed us to conduct systematic investigations of several stimulation configurations to study their effects on hippocampal dynamics. Overall, this model could provide a useful and computationally cost-efficient tool for exploring the mechanisms of deep brain stimulation and help optimize stimulation protocols aimed at alleviating memory disorders.

Understanding brain-behavior relationships requires models capturing the distributed, interactive, and multiscale nature of neural systems. Traditional univariate approaches and single-biomarker models are inherently limited in this context, as they fail to represent dependencies across regions and the hierarchical organization of brain networks. In this study, we propose a graph-based multivariable framework for brain imaging analysis that integrates key organizational principles of brain function-including segregation, integration, modularity, and temporal dynamics-within a unified graph neural network architecture. The framework represents brain data as hierarchical graphs, where node features encode regional activation and temporal variability, and graph structure captures interactions within and between functional modules. The proposed approach is evaluated using functional near-infrared spectroscopy (fNIRS) data as a case study, where subject-specific brain graphs are constructed from task-based recordings acquired shortly after cochlear implant activation to predict speech understanding outcomes one year later. Under leave-one-subject-out validation, the model demonstrates strong predictive performance (R = 0.73, p < 0.001), outperforming previously reported single-biomarker approaches. Perturbation-based analyses further show that predictions are driven by distributed patterns of activity and interaction across regions and modalities, rather than isolated features. These results illustrate the capability of the proposed framework to capture complex brain organization and highlight its potential as a generalizable platform for multivariable analysis and prediction in neuroimaging applications beyond the specific clinical use case considered here.

Learning new information in the presence of distracters and changing conditions requires the ability to adapt. In the brain, this adaptive capability has been linked to dynamic interactions between attention and working memory, which enable the selective filtering of irrelevant input while preserving behaviorally relevant information. Specific neural oscillations have been implicated in this process. Here, we introduce a phenomenological data-driven framework for oscillatory network modeling that learns condition-dependent coupling laws directly from neural recordings and enables inference of condition-dependent directed pathways. We apply our approach to magnetoen-cephalography (MEG) data collected while participants performed a working-memory task with and without distracters. Recall dynamics in the non-distracter condition are first modeled using a linear oscillatory network in which each region of interest is represented by two alpha-band harmonic oscillators. We use universal differential equations (UDE), an extension of neural differential equations, to capture distracter-induced changes in coupling laws. Symbolic regression is then used to interpret the modifications identified by UDE as nonlinear functions, and an additional method is proposed to identify the directed pathway from the newly emerging nonlinear terms in the dynamics of brain regions of interest. Despite inter-subject variability, working memory recall data from all four participants examined under distraction showed the emergence of a pathway from the dorsolateral prefrontal cortex (dlPFC) to the primary visual cortex (V1). This finding is consistent with the established role of the dlPFC in cognitive control and suggests that distracter processing recruits a directed interaction from prefrontal to visual regions. More broadly, our results illustrate that combining linear models whose parameters are learned from the data with universal differential equations augmented by interpretability methods enables the identification of condition-dependent coupling laws, their representation as interpretable mathematical functions, and the discovery of candidate directed pathways underlying adaptive changes in oscillatory networks without requiring strong prior assumptions about the underlying mechanisms.

Classification of motor imagery (MI) tasks through EEG is valuable in brain-computer interfacing and rehabilitation engineering. EEG channels selection for MI task classification is well discussed problem and is challenging due to its combinatorial nature. Most of the existing methods are subject and task-dependent. This paper introduces a subject-independent EEG channel selection. The proposed approach consists of two stages. First, we rank channels based on their divergence from a reference channel Cz. We hypothesize that channels less divergent from Cz are more relevant for MI task classification. In the second stage, we employ a three-stage feature selection and classification model to evaluate the selected channels. It consists of a bandpass filter, followed by common spatial pattern (CSP) filter and three classifiers viz. SVM, 1-NN and 5-NN. Two publicly available datasets viz. PhysioNet and BCI Competition III IVa datasets have been used to assess the method. It performs 15.21\% more than 3Cs and just 2.91\% less than all-channels accuracy with as few as 20/118 channels on BCI Competition data and 19.64\% more than 3Cs on the PhysioNet dataset with 16/64 channels. Empirical comparison implies that the method performs better than classical models such as CSP Rank, fishers rank, and normalized mutual information, significantly. Results support that our hypothesis that divergence between channels and a reference channel Cz can be used as a ranking measure for channel selection.

Introduction Understanding the mechanisms of human neuromaturation constitutes one of the fundamental questions of neuroscience. While it is well described that large-scale brain maturation is initiated within sensorimotor brain regions and progresses to associative cortex, the underlying developmental neurobiology remains to be fully characterized. Animal models have indicated that cortical inhibitory upregulation might be a driver of neurodevelopment. To investigate the hypothesis that cortical inhibitory upregulation plays a similar role in human neuromaturation, we developed a The Virtual Brain (TVB) based computational model (TVB-Child) to explore potential mechanisms of human neurodevelopment. Material and method We created neurodevelopmental dynamic brain network models capturing neurobiological maturation by using the large-scale brain simulator TVB and fitting brain network models to developmental functional MRI (fMRI) from the Human Connectome Project-Development (HCP-D) data set with 640 subjects with an age range of 6-21 years. Age-dependent trajectories in the fMRI data set were first analyzed by combined group-ICA/Dual Regression extracting subject-specific resting-state networks (RSN). Maturational topographical and topological redistribution of these networks were analyzed by linear and non-linear regression of RSN size and degree and strength centrality. Brain network models were fitted to the fMRI functional connectivity obtained from the HCP-D data set. Hypothesizing that cortical inhibition is a driver of neuromaturation, we analyzed spatiotemporal inhibition parameter gradients in the dynamic brain network model for the hypothesized significant correlations with fMRI RSN maturational trajectories. Results While during development frontoparietal (FP) and default mode network (DMN) grew and exhibited an increase in both degree and strength centrality, becoming dominant network hubs, the attention network underwent network pruning with a decrease in size and node degree. The primary sensory network changed little. For the fitted brain network models, we obtained a high degree of reproduction with correlation coefficients between empirical and simulated functional connectivities ranging between 0.80 and 0.95. Values of the feed forward inhibition model parameter wijFFI representing the strength of regional feedforward inhibitory input exhibited the most significant increase with age within the FP and DMN networks. A less pronounced, but significant, age-dependent increase of the inhibitory parameter values were seen in attention networks and no change within primary sensory networks. Conclusion Our study shows that high order (FP, DMN), attention and primary sensory networks exhibit distinct topographical and topological maturation trajectories. Moreover, brain network modeling revealed RSN-specific age-dependent inhibition trajectories, indicating that the model is able to reproduce and thus support candidate mechanisms of neurodevelopment.

Brain-computer interfaces face a fundamental trade-off between the signal fidelity and stimulation precision of noninvasive systems and the surgical burden and scalability of invasive systems. Non-invasive BCIs suffer from low signal quality and poor stimulation accuracy due to the skull barrier and the variability introduced by the scalp and skull. Existing invasive BCIs rely on traumatic surgical procedures or brain-penetrating electrodes, which limits their spatial extensibility, application, and patient acceptance. Here, we introduce a minimally invasive hybrid BCI architecture that uses the skull as a distributed interface layer rather than treating it solely as a barrier. The hybrid BCI comprises four integrated components: (1) the safe and smart micro-hole craniotomy; (2) distributed microelectrodes subcutaneously implanted in micro-holes in the skull with the distal end in contact with the dura; (3) an external bi-directional wearable headset for coupling, recording, stimulation, and channel selection; and (4) an AI-assisted planning and control agent. Animal studies have shown that micro-holes with a diameter of 300-800 m can be safely and conveniently prepared at any predefined locations across the skull without impairing the dura. In vivo experiments on rats demonstrate that the hybrid BCI with skull-implanted microelectrodes evidently increases resting-state spectral power and improves the signal-to-noise ratio of somatosensory and steady-state visual evoked responses compared to the scalp EEG; the computational modelling shows that distributed skull-dura microelectrodes can increase the intracranial electrical field strength and steer focused temporal-interference fields towards predefined deep brain targets. These findings will lay a solid foundation for future endeavors in wireless integration, safety evaluation and clinical benefits of the hybrid BCI. In summary, we propose the hybrid BCI as a distinct minimally invasive BCI paradigm with the great potential as a distributed, scalable, and upgradable neural interface that can expand the clinical application of minimally invasive BCI techniques.

Neural systems exhibit multiple firing states that reflect an organism's internal state and modulate the relationship between external environmental stimuli and behavior. Several studies have inferred these latent states by supplementing the traditional hidden Markov Model (HMM) with generalized linear models (GLMs) with non-Poisson behavioral observations. However, understanding the relationship between internal brain states and behavior also requires modeling the neural activity. Nonetheless, fitting multi-neuron GLM-HMMs is non-trivial due to high sparsity, collinearity, and low trial counts in neuronal datasets. Therefore, we built a robust multi-neuron GLM-HMM framework that uncovers latent states from population activity while incorporating the influence of time-stamped task variables and spike histories. To obtain reliable model parameters, we employ a modified expectation-maximization procedure. Specifically, we show that incorporating neuron-adaptive penalization in the maximization step overcomes the covariate co-linearity issues typical of time-stamped events and sparse spiking, yielding stable estimates of Poisson GLM coefficients. Furthermore, we incorporate a trust-region algorithm to ensure stable M-step convergence in the presence of ill-conditioned Hessians that can lead to unstable Newton-Raphson updates. We further demonstrate the utility of leave-one-out cross-validation analysis for evaluating model performance on datasets with low trial counts and without breaking their temporal structure. We evaluate our framework on three electrophysiological datasets from primates and rodents as they perform a decision-making task, demonstrate stable model convergence, and discuss the behavioral relevance of the inferred states.

Despite extensive study of cellular mechanisms underlying long-term potentiation, no single specific protein or gene has been identified which encodes an individual unit of information, or memory bit. Indeed, the brain engram remains a knowledge gap. The theory of exclusion led us to cancel one-by-one several unrealistic biological options, suggesting that the explanation resides somewhere else. Superposition of up to concentric 300 myelin layers, spiraled, and highly compacted wrapping a single axon and each wrap could host hundreds to thousands of niches, as memory cells, collectively consisting of a massive array of cells. The disjointed 3D spatial superposition allows storage of charges, nodes not facing from a layer to next. The thickness of a single myelin layer ranges from 7.0 to 20 nm. The dimension scale is approximately the exact dimensions of the charge trap, the tunnel and dielectric also equipping current AI microchips. Stored charges are positive ions, with similar effect whether charges are negative or positive charges creating an electromagnetic field. To write data, following an action potential, this voltage applies to the control gates of the myelin layers producing an ionic charge injection. This causes charges to gain energy and tunnel through the myelin layer across Ranvier nodes, via quantum tunneling, and deep into the concentric myelin multilayers. This is creating an insulated trapping of K+ ions isolated from the system. In a long white matter tract bundle, the near-perfect isolation of millions of axons within compressed myelin wrap-ion channel K+/Na+ systems provides quantum coherence and precision of asynchronous firing property. The injected ionic charges (K+) become physically stuck in traps within the myelin layers. The K+ ions may not move freely, completely trapped after AP ceases. Mirroring a single-bit, single-level-cell, a trapped ionic charge (ions K+) may represent a 1, while an empty cell (absence of K+) represents a 0. The trial-and-error process, with a Bayesian inference which may have also been the core evolution of the learning human brain. Based on selected mathematical equations, we analyzed the general scheme on how deep learning may be embedded in the brain

Low-cost and open-source neural recording systems are increasingly important for expanding access to electrophysiological research. However, many existing platforms still rely on specialized hardware or limited modularity, restricting flexibility for laboratories seeking customizable solutions. Here, we developed and evaluated a modular neural recording platform constructed entirely from commercially available components. Recordings were compared against the ground truth. The platform successfully recovered local field potential (LFP)-like waveforms in most conditions and detected spike-like activity during direct connection recordings. Principal component analysis and k-means clustering further demonstrated the ability to distinguish multiple simulated spike waveforms. Signal quality varied across configurations, with saline recordings and preamplifier integration introducing increased noise and reduced detectability. These findings demonstrate the feasibility of building affordable and modular electrophysiology systems using widely accessible hardware. Although the current implementation has limitations in sampling rate, noise performance, and in vivo validation, the presented framework provides a practical foundation for future customizable open-source neural recording.

Capturing naturalistic behavioral dynamics is essential for understanding social interaction in ecologically valid settings. Existing investigations of naturalistic social interaction rely on time-aggregated analysis methods better suited for task-based experiments, which lose the complex, moment-to-moment dynamics exhibited in naturalistic settings. The emerging field of topological data analysis (TDA) provides new tools to characterize fine-grained dynamics in time-series data that cannot be captured by time-averaged methods. The present work utilizes Temporal Mapper, a recently developed TDA specifically tailored to analyzing dynamical systems. Temporal Mapper characterizes complex temporal dynamics as transition networks, where nodes are stable states and edges are transitions between states. Originally designed for human neural time series analysis, here we demonstrate the utility of Temporal Mapper to capture rich animal postural dynamics during naturalistic social interaction. We utilized an existing dataset with 12 video recording sessions of two domestic ferrets (Mustela putorius furo) during naturalistic interaction and tracked the postures of animals during social interaction. Ferrets were chosen due to their strong social-cognitive skills and rich postural dynamics for investigating social behavior via posture estimation. Temporal Mapper was then used to represent the postural dynamics as transition networks for each recording session. Here, we found that posture states are significantly smaller and more widespread during active social interaction compared to non-social activities. Additionally, the number of sequential postural states before transitioning to new behaviors is more consistent during active social interaction than non-social activities. Together, our findings suggest that social activity has a broad range of unstable postural states arranged in consistent sequences. Our method, Temporal Mapper, allows for network structure analysis of complex naturalistic data, applicable for characterizing rich dynamics in different species, scales, and paradigms.

Personalized whole-brain modeling aims to transform treatment planning for neurological disorders by enabling patient-specific simulations of brain network dynamics. Neural mass models (NMMs) offer a tractable compromise between biophysical detail and computational cost and can be directly linked to macroscopic observables such as EEG. However, scaling NMMs to whole-brain networks with realistic connectivity, conduction delays, and cortical surface resolution--and fitting them to individual patient data--imposes computational demands that existing frameworks cannot meet at clinically relevant timescales. Here we introduce BraiNN, a JAX-based Python framework for large-scale neural mass modeling that achieves speedups of up to two to three orders of magnitude over existing tools by leveraging GPU/TPU-accelerated, XLA-compiled array computation. BraiNN combines a region-level Jansen-Rit network with a subject-specific cortical surface mesh of coupled neural mass models and biophysically grounded EEG forward modeling via reciprocity-based lead fields. Its fully differentiable computational graph enables a hybrid personalization pipeline that pairs Bayesian optimization for global parameter exploration with gradient-based refinement, completing EEG-driven spectral fitting of an eight-dimensional parameter space in approximately 2-3 hours on a single consumer GPU--compared to multiple days with conventional neural mass modeling software. Numerical verification against established benchmarks confirms that BraiNN faithfully reproduces canonical synchronization and bifurcation dynamics of Jansen-Rit networks. By reducing the time requirements for personalizing a high-detail whole-brain surface model from days to a few hours on consumer-grade hardware, BraiNN brings personalized brain network modeling closer to practical use in clinical contexts. We anticipate that BraiNN will serve as a foundation for patient-specific digital twins and EEG-guided neuromodulation planning.

Neural decoding can be viewed as a representation learning problem in which neural activity is mapped into an intermediate representation before downstream reconstruction. The choice of intermediate representation influences both performance and learning difficulty. Here we developed a novel framework for studying how intermediate representation choice influences downstream learning and reconstruction. As a proof-of-concept, we instantiated our framework using diffusion latent representations extracted from different diffusion timesteps for neural speech decoding. Component-wise evaluation showed that reconstruction performance differed substantially across diffusion timesteps, with teacher-forced Word Error Rates of 44.7%, 7.5%, and 3.5% for different latent models. These results demonstrate that diffusion latent representations can serve as effective intermediate representations for learning from neural activity, but that their effectiveness depends strongly on the selected diffusion timestep. More broadly, our framework provides a basis for systematically studying how intermediate representation choice influences downstream learning and reconstruction.