Dynamical Systems seminars

Uddipan Banik derives a dynamical account of dark-matter halo density profiles using the Boltzmann–Poisson equations for an expanding Universe populated by subhalos of many masses. A self-similar solution links density slope to halo mass growth; including gravitational scattering and subhalo capture gives gamma = 6/(s + 2). Initial collapse with s = 2 produces a prompt cusp of slope 3/2. Hierarchical growth connects the slope to the primordial power spectrum and recovers the Syer–White relation during matter domination. Stable clustering selects an inner Navarro–Frenk–White cusp of slope one, while stalled accretion approaches the outer slope of three. The framework connects these profiles to successive stages of halo assembly.

Taillte May investigates how unstable black-hole quasinormal-mode spectra affect observed gravitational-wave signals. Small changes in a gravitational potential can substantially shift the frequencies normally used to test astrophysical black holes. A simple time-domain model reveals a delay before those shifts appear, tied to the echo travel time between the perturbation and the main potential. In the perturbative regime, the frequency change after one echo agrees with the conventional frequency-domain prediction. Beyond that regime, the intermediate waveform cannot be represented by a single constant frequency shift.

Large recurrent networks are important models in several fields, including neuroscience, machine learning, physics, and applied mathematics. Yet their dynamics are difficult to study directly, because high-dimensional nonlinear systems can exhibit rich behavior that is hard to summarize in terms of individual trajectories. In this talk, I will discuss an approach that seeks to understand such dynamics through the structure of the network’s equilibria. I will focus on a random balanced network of threshold-linear units that undergoes a transition from a single stable equilibrium to extensive chaos as the disorder strength crosses a critical value. Using a combination of Kac–Rice theory, replica calculations, numerical root-finding, and dynamical mean-field theory, we show that the chaotic regime contains an exponentially large number of equilibria. These equilibria are all saddles, but with only a fractionally small number of unstable directions. Surprisingly, despite the completely random connectivity, the equilibria are not scattered randomly through phase space. Instead, they are strongly correlated and confined to a comparatively small region. The chaotic attractor lies within this same region, suggesting a direct geometric link between the organization of unstable equilibria and the collective structure of the dynamics. This picture helps explain why networks with extensive chaos can nevertheless display dynamics dominated by a relatively small number of collective modes. More broadly, the results suggest that the geometry of equilibria provides a useful complementary perspective to dynamical mean-field theory for understanding high-dimensional neural dynamics. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-05-27. Recording duration: 00:46:40.

Networks of interconnected neurons display diverse patterns of activity. Relating these patterns to the structure of the network is a central goal of theoretical neuroscience. Classic neural field and rate models have been powerful tools for this purpose due to their analytical tractability. Here, we show that the recently-developed combinatorial threshold-linear network (CTLN) model is a mean-field theory for excitatory-inhibitory Hawkes networks, with clustered connectivity, in an inhibition-stabilized regime. This mapping allows us to leverage powerful analytical results for CTLN networks to predict diverse macroscopic dynamics of clustered Hawkes networks, including metastability between various macroscopic fixed points, limit cycles, and chaotic attractors. We will then examine an extension of this approach to models with nonlinear dendritic dynamics, focusing on dendritic calcium spikes.We uncover a marked point process mean-field theory for these n etworks and use this to examine how somatic vs dendritic-targeting connectivity shapes the mean-field equilibrium phase diagram. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-05-13. Recording duration: 00:54:14.

How does a cat gracefully walk and suddenly freeze when spotting a mouse? In this talk, we look at how networks in the spinal cord generate movement. In particular, we address the fundamental yet poorly understood question of motor control: How can rhythmic movements, such as walking, be generated and stopped at any point in the cycle while posture is preserved? Since conventional models of spinal motor function rely on alternation between flexor and extensor modules, which are limited to just two phases, this question exposes the essential shortcoming of the conventional understanding: How can a system with only two phases generate and stop walking in any phase? To address this question and better understand the generation and stopping of motor activity, we use Neuropixels probes in the rat spinal cord during voluntary, freely moving locomotion. We utilize optogenetic activation of a brainstem nucleus to induce stopping. During locomotion, neuronal manifold activity exhibits robust rotational patterns that are topologically invariant with respect to speed (Linden 2022). Furthermore, this trajectory converges on a stable point-attractor precisely at the moment of arrest, and it persists until the movement is resumed. Through computational modeling, we propose that the walk-to-stop represents a bifurcation from a limit cycle to a fixed point attractor. We also propose a structural network mechanism for its physical implementation (Komi 2026). The structural mechanism entails a longitudinal projectome with a skewed Mexican hat topology, i.e., primarily local recurrent excitation and longer-range inhibition. Such a network can generate motor patterns via traveling waves, with frequency and amplitude controlled independently, and rhythm induced without requiring cellular pacemaker mechanisms. Together, our experimental observations support a new theory for the mechanism behind the generation of movement by networks in the spinal cord. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-03-25. Recording duration: 00:37:44.

In sensory systems, stimuli are represented through the diverse firing responses and receptive fields of neurons. These features emerge from the interaction between excitatory (E) and inhibitory (I) neuron populations within the network. Changes in sensory inputs alter this balance, leading to shifts in firing patterns and the input-output properties of individual neurons and the network. While these phenomena have been studied extensively with experiments and theory, the underlying principles for combining E and I inputs are still unclear. Here, the rules for probabilistically combining E and I inputs are derived that describe how neurons in a feedforward inhibitory circuit respond to stimuli. This simple model is broadly applicable, capturing a wide range of response features that would otherwise require multiple separate models and offers insights into the cellular and network mechanisms influencing the input-output properties of neurons, gain modulation, and the emergence of diverse temporal firing patterns. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-02-11. Recording duration: 00:48:34.

Many animals rely on persistent internal representations of continuous angular variables for working memory, motor control, and navigation. Theories have proposed that such representations are maintained by a class of recurrently connected networks called ring attractor networks. These networks rely on large numbers of neurons to maintain continuous and stable representations and to accurately integrate incoming signals. The head direction system of the fruit fly, however, seems to achieve these properties with a remarkably small network. These findings challenge our understanding of ring attractors and their putative implementation in neural circuits. In this talk, I will show analytically how small networks can overcome the constraints of their size to generate a ring attractor and are hence capable of stably maintaining an internal representation of a continuous, periodic variable. Further, I will show how ring attractors emerge in small threshold linear networks through the coordination of a discrete set of line attractors. More broadly, this work informs our understanding of the functional capabilities of small, discrete systems. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-02-04. Recording duration: 00:44:59.

Computations in neural circuits are often construed as being implemented through the coordinated dynamics of neurons. In this picture, the role of synaptic connectivity is to sculpt neuronal dynamics to implement computations of interest. Of course, synapses are not static but change on a variety of timescales, including fast timescales comparable to those of neurons. Thus, a more accurate view of computation in neural circuits may involve the coupled dynamics of neurons and synapses. This form of computation is closer to what is implemented by Transformers via an equivalence between ongoing synaptic plasticity and self-attention. I will first describe a nonlinear recurrent neural-network model with ongoing Hebbian dynamics of “fast” synapses atop unstructured “slow” synapses. I will then describe two computations implemented through neuronal-synaptic dynamics, which can be studied in this model using techniques including dynamical mean-field theory and random-matrix theory. First, there exists a novel phase termed “freezable chaos” in which a stable fixed point of neuronal dynamics is continuously destabilized by synaptic dynamics. This allows for the creation of a stable fixed point at any neuronal state visited by the network by halting synaptic plasticity. Second, I will describe an effect termed “persistent oscillations” in which, following stimulation by a periodic signal, a plastic network continues to autonomously reproduce a similar signal for a duration exceeding any intrinsic timescale in the system. Thus, ongoing Hebbian plasticity can provide a dynamic form of working memory, complementing the static form provided by freezable chaos. Ongoing experimental work suggests that this effect is realized in cortical organoids. Overall, this line of work suggests that synapses should be promoted to first-class dynamical degrees of freedom in our conceptual understanding of neural-circuit function. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2026-01-14. Recording duration: 00:43:14.

Recurrent networks with balanced excitation-inhibition explain a wide range of neurophysiological observations, but can only implement a limited set of transformations on their input. On the other hand networks of firing-rate units with low-rank connectivity have universal computational capabilities, but do not work with spikes or generate noise self-consistently. Although empirical approaches to merge these two computational frameworks have been constructed, there is no established theory describing their unification. Here we develop such a theory. We study analytically and numerically networks with connectivity comprising random “strong”, and low-rank “weak” components. When the low-rank connectivity is slow, a well-defined notion of instantaneous firing rate emerges which implies universal computation as previously shown. However, the fact that such time-varying rates are the result of E-I balance has important implications. We show that internally or externally generated fluctuations along particular latent modes tend to break the E-I balance. Its maintenance is obtained through the emergence of a spontaneous coupling between the mean and the variance of the membrane potential and the norm of the latent state driving these modes. This leads to several predictions, the most counterintuitive of which is that coherent global fluctuations in subthreshold membrane potential (Vm) should coexist with desynchronized activity at constant firing rates when the dynamics of these modes is excited. To test our theory, we show that the coupling between the average Vm and the latent state adds new non-linear dimensions to the low-dimensional manifold of the network, which lead to a frequency doubling when the input to the network is periodic, a prediction that is borne out in population recordings from mouse V1. Our results unify two prevalent frameworks for cortical computation and clarify the relationship between computation, dynamics and geometry in circuits of spiking neurons. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-11-26. Recording duration: 00:39:26.

Minimising cortical prediction errors is thought to be a key computation underlying perception, action, and learning. Yet, how the cortex represents and uses uncertainty in this process remains unclear. In the first part of this talk, I will present a normative framework showing how uncertainty can modulate prediction error activity to yield uncertainty-modulated prediction errors (UPEs), hypothesised to be represented by layer 2/3 pyramidal neurons. We propose that these UPEs are computed through inhibitory mechanisms involving SST and PV interneurons. A circuit model demonstrates how cortical cell types can locally compute means, variances, and UPEs, leading to adaptive learning rates. In the second part, I will discuss how uncertainty modulation could be controlled by higher-level representations. We formally derived neural dynamics that minimise prediction errors under the assumption that cortical areas must not only predict the activity in other areas and sensory streams but also jointly project their inverse expected uncertainty about their predictions, which we call “confidence”. This yields a confidence-weighted integration of bottom-up and top-down signals, consistent with Bayesian principles, and predicts the existence of second-order errors that compare confidence with performance. We predict that these second-order errors propagate alongside classical prediction errors through the cortical hierarchy, and simulations demonstrate that this mechanism enables nonlinear classification within a single cortical area. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-11-12. Recording duration: 00:29:21.

We introduce a statistical physics framework for learning in neural architectures composed of single or interconnected asymmetric attractor networks. These systems can exhibit a manifold of global fixed points capable of implementing sophisticated input-output mappings, which we characterize analytically. Learning from extensive datasets is achieved through the stabilization of fixed points via a fully distributed and local learning process, implemented at the single-neuron level. This simple mechanism yields performance comparable to that of conventional feedforward deep neural networks trained using gradient-based methods. The effectiveness of the model stems from the dense and accessible manifolds of stable fixed points, which encode the internal representations of data. Unlike other approaches to deep learning without backpropagation, our method does not attempt to estimate gradients. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-06-18. Recording duration: 00:45:35.

For decades dynamical systems theory played a pivotal role in theoretical and computational neuroscience, as it links biophysical and biochemical processes to neural computation. In fact, dynamical systems are computationally universal. Rather than hand-crafting computational theories of neural function based on dynamical systems, recent developments in scientific machine learning (ML) and AI suggest that we may be able to infer such dynamical-computational models directly from neurophysiological and behavioral observations. This is called dynamical systems reconstruction (DSR), the learning of generative surrogate models of the underlying dynamics, including its long-term temporal and geometrical properties, from time series data. In my talk I will cover recent ML/AI architectures, training algorithms, and validation procedures for DSR. I will discuss specifically how recent AI architectures for DSR can integrate neuroscience data from multiple modalities (like multiple single-unit recordings and behavioral choices), across diverse time scales, and across many different animals and task designs, into a joint DSR model. This provides first steps toward dynamical systems based AI foundation models for neuroscience. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-04-23. Recording duration: 00:53:24.

The use of electronic circuits to model neural systems goes back to C. Mead and is present in models, from leaky-integrate-and-fire to Hodking-Huxley. Simulating neural network with analog hardware is attractive: it allows to implement neurocomputations in real time without discretization approximations, it has perfect simulation-time scaling with system size, and it provides ready-to-deploy neuromorphic circuit for applications. There are implementations in CMOS technology, however, they are complex, require sophisticated fabrication facilities and, most important, suffer from significant device mismatch. In a radically different approach, based on the concept of memristors, we introduce a neuro-synaptic circuit of unprecedented simplicity, with readily available cheap off-the-shelf electronic components, that can quantitatively reproduce textbook theoretical neuron and synaptic current models. Our neuron circuits can avoid the mismatch problem and are easily tuneable at bio-compatible time-scales. We first introduce a voltage-gated conductance bursting neuron model that produces spike traces that bare striking similarity to experimental recordings. We then introduce synaptic current circuits and show the modularity of our method implementing neurocomputing primitives of basic network motifs, including CPGs. With this "theoretical hardware" approach we show: (i) that neuron adaptation and self-excitation can be viewed as a self-consistent dynamical problem; (ii) that a dynamical memory can be minimally implemented with a single recursive spiking neuron; (iii) that an adaptive membrane current reveals a connection between bursting and the driven harmonic oscillator, perhaps pointing to a neural correlate of the pendular limb motion. Finally we discuss the limitation of the approach to networks of mid-size and its potential application for brain-machine-interfaces, robotics and AI. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-03-19. Recording duration: 00:43:46.

In the brain, while early sensory areas encode and process external inputs rapidly, higher-association areas are endowed with slow dynamics to benefit information accumulation over time. This property raises the question of why diverse timescales are well localized rather than being mixed up across the cortex, despite high connection density and an abundance of feedback loops that support reliable signal propagation. In this talk, we will address this question by analyzing a large-scale network model of the primate cortex, and we identify a novel dynamical regime termed "interference-free propagation". In this regime, the mean components of the synaptic currents to each downstream area are imbalanced to ensure signals to propagate reliably, while the temporally fluctuating components of the synaptic inputs governed by upstream areas' timescales are largely canceled out, leading to the localization of its own timescale in each downstream area. Our result provides new insights into the operational regime of the cortex, leading to the coexistence of hierarchical timescale localization and reliable signal propagation. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-02-26. Recording duration: 00:48:02.

Recent advances in techniques for monitoring and perturbing neural populations have greatly enhanced our ability to study circuits in the brain. In particular, two-photon holographic optogenetics now enables precise photostimulation of experimenter-specified groups of individual neurons, while simultaneous two-photon calcium imaging enables the measurement of ongoing and induced activity across the neural population. Despite the enormous space of potential photostimulation patterns and the time-consuming nature of photostimulation experiments, very little algorithmic work has been done to determine the most effective photostimulation patterns for identifying the neural population dynamics. Here, I will discuss ongoing development of active learning techniques to efficiently select which neurons to stimulate such that the resulting neural responses will best inform a dynamical model of the neural population activity. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-02-05. Recording duration: 00:42:29.

Functional connectivity has been the focus of many research groups aiming to study the interaction between cells and brain regions. A standard method for analyzing connectivity is to statistically compare pairwise interactions between cells or brain regions across behavioral states or conditions. This methodology ignores the intrinsic properties of functional connectivity as a multivariate and dynamic signal, expressing the correlational configuration of the network. In this talk, I will present a geometric approach, combining Graph Theory and Riemannian Geometry to build "a graph of graphs" and extract the latent dynamics of the overall correlational structure. Using this approach, we formulate the statistical relations between network dynamics and spontaneous behavior as a second-order Taylor’s expansion. Our analysis shows that fast fluctuations in functional connectivity of large-scale cortical networks are closely linked to variations in behavioral metrics related to the arousal state. We further expand this methodology to longer time scales to study the effect of dopamine on network dynamics in the primary motor cortex (M1) during learning. We developed a series of analysis methods indicating that as animals learn to perform a motor task, the network of pyramidal neurons in layer 2-3 gradually and monotonically reorganizes toward an "expert" configuration. Our results highlight the critical role of dopamine in driving synaptic plasticity: Blocking dopaminergic neurotransmission locally in M1 prevented motor learning at the behavioral level and concomitantly halted plasticity changes in network activity and in functional connectivity. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-01-29. Recording duration: 00:27:08.

Networks of excitatory and inhibitory (EI) neurons form a canonical circuit in the brain. Classical theoretical analyses of dynamics in EI networks have revealed key principles such as EI balance or paradoxical responses to external inputs. These seminal results assume that synaptic strengths depend on the type of neurons they connect but are otherwise statistically independent. However, recent synaptic physiology datasets have uncovered connectivity patterns that deviate significantly from independent connection models. Simultaneously, studies of task-trained recurrent networks have emphasized the role of connectivity structure in implementing neural computations. Despite these findings, integrating detailed connectivity structures into mean-field theories of EI networks remains a substantial challenge. In this talk, I will outline a theoretical approach to understanding dynamics in structured EI networks by employing a low-rank approximation based on an analytical computation of the dominant eigenvalues of the full connectivity matrix. I will illustrate this approach by investigating the effects of pair-wise connectivity motifs on linear dynamics in EI networks. Specifically, I will present recent results demonstrating that an over-representation of chain motifs induces a strong positive eigenvalue in inhibition-dominated networks, generating a potential instability that challenges classical EI balance criteria. Furthermore, by examining the effects of external input, we found that chain motifs can, on their own, induce paradoxical responses, wherein an increased input to inhibitory neurons leads to a counterintuitive decrease in their activity through recurrent feedback mechanisms. Altogether, our theoretical approach opens new avenues for relating recorded connectivity structures with dynamics and computations in biological networks. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2025-01-22. Recording duration: 00:47:27.

Classical views suggest that memories are stored in assemblies of excitatory neurons that become strongly interconnected during learning. However, recent experimental and theoretical results have challenged this view, leading to the hypothesis that memories are encoded in assemblies containing both excitatory (E) and inhibitory (I) neurons. Understanding the effects of these E-I assemblies on memory function is therefore essential. Using a biologically constrained model of an olfactory memory network, I will first describe how introducing E-I assemblies reorganizes odor-evoked activity patterns in neural state space. Indeed, the “geometry” of neural activity provides valuable insights about the computational properties of neural networks. I will then describe the behavior of networks with E-I assemblies upon partial manipulation of inhibitory neurons. Finally, I will discuss recent experimental data supporting predictions of the model. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2024-12-18. Recording duration: 00:35:13.

John Hopfield and Geoffrey Hinton explain how ideas from physics helped establish computational models that store patterns and learn from data. Hopfield develops a physical perspective on collective computation, including associative memories in which network dynamics recover stored patterns from incomplete or distorted inputs. Hinton examines Boltzmann machines, connecting probabilistic neural networks and learning to concepts from statistical physics. Together, the lectures link energy landscapes, interacting units and stochastic behaviour with mechanisms for representation and computation. They provide the scientific background to the neural-network contributions recognised by the 2024 physics prize and the conceptual route from models of collective systems to machine learning.

Continuous attractors offer a unique class of solutions for storing continuous-valued variables in recurrent system states for indefinitely long time intervals. Unfortunately, continuous attractors suffer from severe structural instability in general---they are destroyed by most infinitesimal changes of the dynamical law that defines them. This fragility limits their utility especially in biological systems as their recurrent dynamics are subject to constant perturbations. We observe that the bifurcations from continuous attractors in theoretical neuroscience models display various structurally stable forms. Although their asymptotic behaviors to maintain memory are categorically distinct, their finite-time behaviors are similar. We build on the persistent manifold theory to explain the commonalities between bifurcations from and approximations of continuous attractors. Fast-slow decomposition analysis uncovers the existence of a persistent slow manifold that survives the seemingly destructive bifurcation, relating the flow within the manifold to the size of the perturbation. Moreover, this allows the bounding of the memory error of these approximations of continuous attractors. Finally, we train recurrent neural networks on analog memory tasks to support the appearance of these systems as solutions and their generalization capabilities. Therefore, we conclude that continuous attractors are functionally robust and remain useful as a universal analogy for understanding analog memory. Presented in the van Vreeswijk Theoretical Neuroscience Seminar series (formerly WWTNS) on 2024-11-27. Recording duration: 00:59:01.