University of Oxford

Gaussian fields are prevalent throughout mathematics and the sciences, for instance in physics (wave-functions of high energy electrons), astronomy (cosmic microwave background radiation) and probability theory (connections to SLE, random tilings etc). Despite this, the geometry of such fields, for instance the connectivity properties of level sets, is poorly understood. In this talk I will discuss methods of extracting geometric information about levels sets of a planar Gaussian field through discrete observations of the field. In particular, I will present recent work that studies three such discretisation schemes, each tailored to extract geometric information about the levels set to a different level of precision, along with some applications.

In this talk, I will talk about my current approach to model customer movements and in particular congestion inside supermarkets using queuing networks. As the research question for my project is ‘How should one design supermarkets to minimize congestion?’, I will then talk about my current progress in understanding how the network structure can affect this dynamics.

What makes cities successful? A complex systems approach to modelling urban economies

Urban centres draw a diverse range of people, attracted by opportunity, amenities, and the energy of crowds. Yet, while benefiting from density and proximity of people, cities also suffer from issues surrounding crime, congestion and density. Seeking to uncover the mechanisms behind the success of cities using novel tools from the mathematical and data sciences, this work uses network techniques to model the opportunity landscape of cities. Under the theory that cities move into new economic activities that share inputs with existing capabilities, path dependent industrial diversification can be described using a network of industries. Edges represent shared necessary capabilities, and are empirically estimated via flows of workers moving between industries. The position of a city in this network (i.e., the subnetwork of its current industries) will determine its future diversification potential. A city located in a central well-connected region has many options, but one with only few peripheral industries has limited opportunities.

We develop this framework to explain the large variation in labour formality rates across cities in the developing world, using data from Colombia. We show that, as cities become larger, they move into increasingly complex industries as firms combine complementary capabilities derived from a more diverse pool of workers. We further show that a level of agglomeration equivalent to between 45 and 75 minutes of commuting time maximizes the ability of cities to generate formal employment using the variety of skills available. Our results suggest that rather than discouraging the expansion of metropolitan areas, cities should invest in transportation to enable firms to take advantage of urban diversity.

This talk will be based on joint work with Eduardo Lora and Andres Gomez at Harvard University.

Hamilton-Jacobi-Bellman equations for dynamic pricing

I will discuss the Hamilton-Jacobi-Bellman (HJB) equation, which is a nonlinear, second-order, terminal value PDE problem. The equation arises in optimal control theory as an optimality condition.

Consider a dynamic pricing problem: over a given period, what is the best strategy to maximise revenues and minimise the cost of unsold items?

This is formulated as a stochastic control problem in continuous time, where we try to find a function that controls a stochastic differential equation based on the current state of the system.

The optimal control function can be found by solving the corresponding HJB equation.

I will present the solution of the HJB equation using a toy problem, for a risk-neutral and a risk-averse decision maker.

Cells adapt their metabolic state in response to changes in the environment.  I will present a systematic framework for the construction of flux graphs to represent organism-wide metabolic networks.  These graphs encode the directionality of metabolic fluxes via links that represent the flow of metabolites from source to target reactions.  The weights of the links have a precise interpretation in terms of probabilities or metabolite flow per unit time. The methodology can be applied both in the absence of a specific biological context, or tailored to different environmental conditions by incorporating flux distributions computed from constraint-based modelling (e.g., Flux-Balance Analysis). I will illustrate the approach on the central carbon metabolism of Escherichia coli, revealing drastic changes in the topological and community structure of the metabolic graphs, which capture the re-routing of metabolic fluxes under each growth condition.

By integrating Flux Balance Analysis and tools from network science, our framework allows for the interrogation of environment-specific metabolic responses beyond fixed, standard pathway descriptions.

The immune system finds very rare amounts of pathogens and responds against them in a timely and efficient manner. The time to find and respond against pathogens does not vary appreciably with the size of the host animal (scale invariant search and response). This is surprising since the search and response against pathogens is harder in larger animals.

The first part of the talk will focus on using techniques from computer science to solve problems in immunology, specifically how the immune system achieves scale invariant search and response. I use machine learning techniques, ordinary differential equation models and spatially explicit agent based models to understand the dynamics of the immune system. I will talk about Hierarchical Bayesian non-linear mixed effects models to simulate immune response in different species.

The second part of the talk will focus on taking inspiration from the immune system to solve problems in computer science. I will talk about a model that describes the optimal architecture of the immune system and then show how architectures and strategies inspired by the immune system can be used to create distributed systems with faster search and response characteristics.

I argue that techniques from computer science can be applied to the immune system and that the immune system can provide valuable inspiration for robust computing in human engineered distributed systems.

The non-local Fisher KPP equation is used to model non-local interaction and competition in a population. I will discuss recent work on solutions of this equation with a compactly supported initial condition, which strengthens results on the spreading speed obtained by Hamel and Ryzhik in 2013. The proofs are probabilistic, using a Feynman-Kac formula and some ideas from Bramson's 1983 work on the (local) Fisher KPP equation.

Patlak-Keller-Segel equations 
\[
\begin{aligned}
u_t - L u &= - \mathop{\text{div}\,} (u \nabla v) \\
v_t - \Delta v &= u,
\end{aligned}
\]
where L is a dissipative operator, stem from mathematical chemistry and mathematical biology.
Their variants describe, among others, behaviour of chemotactic populations, including feeding strategies of zooplankton or of certain insects. Analytically, Patlak-Keller-Segel equations reveal quite rich dynamics and a delicate global smoothness vs. blowup dichotomy. 
We will discuss smoothness/blowup results for popular variants of the equations, focusing on the critical cases, where dissipative and aggregative forces seem to be in a balance. A part of this talk is based on joint results with Rafael Granero-Belinchon (Lyon).

Abstract: Triangulations of surfaces serve as important examples for cluster theory, with the natural operation of “diagonal flips” encoding mutation in cluster algebras and categories. In this talk we will focus on the combinatorics of mutation on marked surfaces with infinitely many marked points, which have gained importance recently with the rising interest in cluster algebras and categories of infinite rank. In this setting, it is no longer possible to reach any triangulation from any other triangulation in finitely many steps. We introduce the notion of mutation along infinite admissible sequences and show that this induces a preorder on the set of triangulations of a fixed infinitely marked surface. Finally, in the example of the completed infinity-gon we define transfinite mutations and show that any triangulation of the completed infinity-gon can be reached from any other of its triangulations via a transfinite mutation. The content of this talk is joint work with Karin Baur.

We consider the classic problem of establishing a statistical ranking of a set of n items given a set of inconsistent and incomplete pairwise comparisons between such items. Instantiations of this problem occur in numerous applications in data analysis (e.g., ranking teams in sports data), computer vision, and machine learning. We formulate the above problem of ranking with incomplete noisy information as an instance of the group synchronization problem over the group SO(2) of planar rotations, whose usefulness has been demonstrated in numerous applications in recent years. Its least squares solution can be approximated by either a spectral or a semidefinite programming (SDP) relaxation, followed by a rounding procedure. We perform extensive numerical simulations on both synthetic and real-world data sets (Premier League soccer games, a Halo 2 game tournament and NCAA College Basketball games) showing that our proposed method compares favorably to other algorithms from the recent literature.

We propose a similar synchronization-based algorithm for the rank-aggregation problem, which integrates in a globally consistent ranking pairwise comparisons given by different rating systems on the same set of items. We also discuss the problem of semi-supervised ranking when there is available information on the ground truth rank of a subset of players, and propose an algorithm based on SDP which recovers the ranks of the remaining players. Finally, synchronization-based ranking, combined with a spectral technique for the densest subgraph problem, allows one to extract locally-consistent partial rankings, in other words, to identify the rank of a small subset of players whose pairwise comparisons are less noisy than the rest of the data, which other methods are not able to identify.