Semester 1, 2024-25

High-order functional interactions in health and therapy: aging and transcranial ultrasound stimulation

Marylin Gatica, Northeastern University London

Monday 28 October 2024, 13:00-14:00
Watson B16

The brain interdependencies can be studied from either a structural or functional perspective. The former focuses typically on structural connectivity (SC), while the second considers statistical interactions (usually functional connectivity, FC). While SC is inherently pairwise because it describes white-matter fibers projecting from one region to another, FC is not limited to pairwise interdependencies. Despite this, FC analyses predominantly concentrate on pairwise statistics, usually neglecting the possibility of higher-order interactions. Moreover, the precise relationship between high-order and SC is largely unknown, partly due to the absence of mechanistic models that can efficiently map brain connectomics to functional connectivity.

To investigate these interlinked issues, we have built whole-brain computational models using anatomical and functional MRI data in two applications: healthy aging and transcranial ultrasound stimulation (TUS). We show that non-linear variations in the structural connectome can largely explain the differences in high-order functional interactions between age groups. Moreover, we showed the extent of perturbations in dynamical models to describe the high-order effects of TUS in two different brain targets.

Variational inference for stochastic differential equations driven by fractional Brownian motion

Manfred Opper, TU Berlin and University of Birmingham

Monday 4 November 2024, 13:00-14:00
Arts 201

Stochastic differential equations (SDE) driven by white noise are important models for stochastic dynamical systems in natural science and engineering. The statistical inference of the parameters of such models based on noisy observations has also attracted considerable interest in the machine learning community. Using Girsanov's change of measure approach one can apply powerful variational techniques to solve the inference problem. A limitation of standard SDE models is the fact that they show typically a fast decay of correlation functions. If one is interested in stochastic processes with a long-time memory, a well-known possibility is to replace the Brownian motion in the SDE by the so called fractional Brownian motion (fBM) which is no longer a Markov process. Unfortunately, variational inference for this case is much less straightforward. Our approach to this problem utilises a somewhat overlooked idea by Carmona and Coutin (1998) who showed that fBM can be exactly represented as an infinite-dimensional linear combination of Ornstein-Uhlenbeck processes with different time constants. Using an appropriate discretisation, we arrive at a finite dimensional approximation which is an 'ordinary' SDE model in an augmented space. For this new model we can apply (more or less) off-the shelve variational inference approaches.

A guide to graph-based learning

Jeremy Budd, University of Birmingham

Monday 4 November 2024, 14:00-15:00
Arts 201

In this talk, I will give an overview of the field of graph-based learning, a field that has matured over the last 15 years and is rich in both practical applications and theoretical underpinnings. The key idea of graph-based learning is to understand interrelated data as a graph, to solve variational problems and PDEs on that graph to analyse that data, and to study the limits of such models as the number of nodes goes to infinity. I will begin by motivating the approach and then will discuss the mathematical framework, three classic methods in the field, the nuances of implementing these methods, and finally the theoretical underpinnings of this field.

Static friction models, buckling and lift-off for a rod deforming on a cylinder

Rehan Shah, Queen Mary, University of London

Monday 18 November 2024, 13:00-14:00
Watson B16

In this talk, I will develop a comprehensive geometrically-exact theory for an end-loaded elastic rod constrained to deform on a cylindrical surface. By viewing the rod-cylinder system as a special case of an elastic braid, it will be shown that all forces and moments imparted by the deforming rod to the cylinder as well as all contact reactions can be obtained. This framework allows us to give a complete treatment of static friction consistent with force and moment balance. In addition to the commonly considered model of hard frictionless contact, I analyse two friction models in which the rod, possibly with intrinsic curvature, experiences either lateral or tangential friction. Applications of the theory include studying buckling of the constrained rod under compressive and torsional loads, finding critical loads to depend on Coulomb-like friction parameters, as well as the tendency of the rod to lift off the cylinder under further loading. The cylinder can also have arbitrary orientation relative to the direction of gravity. The cases of a horizontal and vertical cylinder, with gravity having only a lateral or axial component, are amenable to exact analysis, while numerical results map out the transition in buckling mechanism between the two extremes. Weight has a stabilising effect for near-horizontal cylinders, while for near-vertical cylinders it introduces the possibility of buckling purely due to self-weight. The results are relevant for many engineering and medical applications in which a slender structure winds inside or outside a cylindrical boundary.

Making a deposit: advances in our understanding of droplet evaporation and the coffee ring effect

Madeleine Moore, Loughborough University

Monday 25 November 2024, 13:00-14:00
Watson B16

The evaporation of liquid droplets has received significant research interest due to its fundamental significance in a variety of industrial and engineering applications such as inkjet printing, microscale and colloidal patterning, DNA microarray technologies and the manufacture of Q/OLEDs. One of the key reasons for this is the familiar ‘coffee-ring’ effect that refers to the ringlike stain left behind after a solute-laden droplet evaporates on a surface and its potential use in depositing specific patterns. While deceptively simple, there is a wealth of complexity in the problem, primarily embedded in the – potentially coupled – aspects of evaporation, the associated liquid flow and particle transport. These difficulties have limited the vast majority of existing models to only treating the simplest possible cases of asymptotically flat, circular droplets evaporating in isolation. This has dramatically limited their applicability in real-world contexts, in which these simplifications are generally broken. In this talk, we will discuss recent advances that attempt to broaden the existing theory with an eye on the ultimate goal of dynamically controlling the process to suit a specific application.

A novel use of pseudospectra in mathematical biology: Understanding HPA axis sensitivity

Catherine Drysdale, University of Birmingham

Monday 2 December 2024, 13:00-14:00
Watson Building B16

The Hypothalamic-Pituitary-Adrenal (HPA) axis is a major neuroendocrine system, and its dysregulation is implicated in various diseases. This system also presents interesting mathematical challenges for modelling. We consider a non-linear delay differential equation model and calculate pseudospectra of three different linearisations: a time-dependent Jacobian, linearisation around the limit cycle, and dynamic mode decomposition (DMD) analysis of Koopman operators (global linearisation). The time-dependent Jacobian provided insight into experimental phenomena, explaining why rats respond differently to perturbations during corticosterone secretion's upward versus downward slopes. We developed new mathematical techniques for the other two linearisations to calculate pseudospectra on Banach spaces and apply DMD to delay differential equations, respectively. These methods helped establish local and global limit cycle stability and study transients. Additionally, we discuss using pseudospectra to substantiate the model in experimental contexts and establish bio-variability via data-driven methods. This work is the first to utilise pseudospectra to explore the HPA axis.

Higher-order organisation of multivariate time series

Enrico Amico, University of Birmingham

Monday 2 December 2024, 14:00-15:00
Watson Building B16

Time series analysis has proven to be a powerful method to characterise several phenomena in biology, neuroscience and economics, and to understand some of their underlying dynamical features. Several methods have been proposed for the analysis of multivariate time series, yet most of them neglect the effect of non-pairwise interactions on the emerging dynamics. In this talk I will introduce a framework to characterise the temporal evolution of higher-order dependencies within multivariate time series. Using network analysis and topology, I will show that our framework robustly differentiates various spatiotemporal regimes of coupled chaotic maps. This includes chaotic dynamical phases and various types of synchronisation. Furthermore, using the higher-order co-fluctuation patterns in simulated dynamical processes as a guide, I will highlight and quantify signatures of higher-order patterns in data from brain functional activity, financial markets and epidemics. Overall, this approach sheds light on the higher-order organisation of multivariate time series, allowing a better characterization of dynamical group dependencies inherent to real-world data.