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.