Dynamical Systems

Marty Weissman presents a University of California, Santa Cruz course that develops mathematics for the life sciences through models of biological dynamics. Adapted from UCLA’s LS30 course, it integrates calculus, dynamical systems and linear algebra. The seminar describes the curriculum and its active learning discussion sections, where deliberate choices range from custom dynamical-system simulators to paper exercises and blackboard work. Drawing on two years of teaching, Weissman reflects on which approaches appear to support student learning and the development of the teaching team.

How do humans stabilize an inverted pendulum, and why does a balanced pole eventually fall? Drawing on 25 years of fingertip pole-balancing research, this talk examines neural correction delays: longer poles move more slowly relative to the nervous system’s response time. Delay-differential models can stabilize the upright position, yet skilled people still experience falls. The proposed explanation is microchaos arising from interactions among delay, sensory dead zones and frequency-dependent force encoding. A region of transient falling solutions lies next to stable microchaotic dynamics. Such microchaos is absent in virtual frontal-plane balancing tasks, while models of standing postural sway lack the corresponding transient regime. The comparison suggests that human falls, unlike pole falls, are more plausibly associated with medical events or slips and trips.

The brain adjusts its internal activity to direct visual and auditory attention and regulate responses to external stimuli. Appropriate endogenous activity may support normal cognition and behaviour, while its modulation can also produce dysfunction. Beginning with mathematical models of general anaesthesia, this talk examines how changes in noisy endogenous activity can fragment functional brain organization and account for loss of consciousness. Anaesthesia experiments in ferrets support the theoretical result. Further experimental examples are discussed to show how underlying internal noise may help explain observed transitions between brain states.

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.

Theoretical postdoctoral positions investigating hybrid quantum light-matter systems, including quantum-gas cavity QED and polaritonic chemistry; requires a PhD in physics with expertise in quantum optics and open quantum systems. Applications reviewed on a rolling basis until the hard deadline.

The brain adjusts its internal activity to direct visual and auditory attention and regulate responses to external stimuli. Appropriate endogenous activity may support normal cognition and behaviour, while its modulation can also produce dysfunction. Beginning with mathematical models of general anaesthesia, this talk examines how changes in noisy endogenous activity can fragment functional brain organization and account for loss of consciousness. Anaesthesia experiments in ferrets support the theoretical result. Further experimental examples are discussed to show how underlying internal noise may help explain observed transitions between brain states.

How do humans stabilize an inverted pendulum, and why does a balanced pole eventually fall? Drawing on 25 years of fingertip pole-balancing research, this talk examines neural correction delays: longer poles move more slowly relative to the nervous system’s response time. Delay-differential models can stabilize the upright position, yet skilled people still experience falls. The proposed explanation is microchaos arising from interactions among delay, sensory dead zones and frequency-dependent force encoding. A region of transient falling solutions lies next to stable microchaotic dynamics. Such microchaos is absent in virtual frontal-plane balancing tasks, while models of standing postural sway lack the corresponding transient regime. The comparison suggests that human falls, unlike pole falls, are more plausibly associated with medical events or slips and trips.

Navid Azizan presents methods that make neural models obey physical, safety and operational constraints at deployment. Hard-constrained neural networks, or HardNets, enforce input-dependent constraints by construction while preserving universal approximation within the feasible function class. Applications include models of chaotic dynamics with bounded trajectories, energy-constrained operator learning, safe reinforcement learning and control with formal guarantees. The talk then turns to enforcing constraints during sampling from pretrained diffusion and flow-matching models. Formulating generation as trajectory optimization allows receding-horizon control to guide outputs toward feasibility without retraining or imposing excessive restrictions on the sampling process. Examples from fluid dynamics, robot planning and control, PDE control and language-guided image editing illustrate how constraints can define admissible behavior while retaining expressive learning and generation.

Marty Weissman presents a University of California, Santa Cruz course that develops mathematics for the life sciences through models of biological dynamics. Adapted from UCLA’s LS30 course, it integrates calculus, dynamical systems and linear algebra. The seminar describes the curriculum and its active learning discussion sections, where deliberate choices range from custom dynamical-system simulators to paper exercises and blackboard work. Drawing on two years of teaching, Weissman reflects on which approaches appear to support student learning and the development of the teaching team.