Independent sets in bipartite graphs and approximating the partition function of the ferromagnetic Potts model
Abstract
This talk considers the problem of sampling an independent set uniformly at random from a bipartite graph (equivalently, the problem of approximately counting independent sets in a bipartite graph). I will start by discussing some natural Markov chain approaches to this problem, and show why these lead to slow convergence. It turns out that the problem is interesting in terms of computational complexity – in fact, it turns out to be equivalent to a large number of other problems, for example, approximating the partition function of the “ferromagnetic Ising model’’ (a 2-state particle model from statistical physics) in the presence of external fields (which are essentially vertex weights). These problems are all complete with respect to approximation-preserving reductions for a logically-defined complexity class, which means that if they can be approximated efficiently, so can the entire class. In recent work, we show some connections between this class of problems and the problem of approximating the partition function of the ``ferromagnetic Potts model’’ which is a generalisation of the Ising model—our result holds for q>2 spins. (This corresponds to the approximation problem for the Tutte polynomial in the upper quadrant
above the hyperbola q=2.) That result was presented in detail at a recent talk given by Mark Jerrum at Oxford’s one-day meeting in combinatorics. So I will just give a brief description (telling you what the Potts model is and what the result is) and then conclude with some more recently discovered connections to counting graph homomorphisms and approximating the cycle index polynomial.
Parallel stochastic simulation using graphics processing units for the Systems Biology Toolbox for MATLAB
Abstract
Graphics processing units (GPU) are well suited to decrease the
computational in-
tensity of stochastic simulation of chemical reaction systems. We
compare Gillespie’s
Direct Method and Gibson-Bruck’s Next Reaction Method on GPUs. The gain
of the
GPU implementation of these algorithms is approximately 120 times faster
than on a
CPU. Furthermore our implementation is integrated into the Systems
Biology Toolbox
for Matlab and acts as a direct replacement of its Matlab based
implementation.
Toposes in algebraic quantum theory
Abstract
Topology can be generalised in at least two directions: pointless
topology, leading ultimately to topos theory, or noncommutative
geometry. The former has the advantage that it also carries a logical
structure; the latter captures quantum settings, of which the logic is
not well understood generally. We discuss a construction making a
generalised space in the latter sense into a generalised space in the
former sense, i.e. making a noncommutative C*-algebra into a locale.
This construction is interesting from a logical point of view,
and leads to an adjunction for noncommutative C*-algebras that extends
Gelfand duality.
A model of crystal growth with corner regularization
Abstract
We investigate a dynamic model of two dimensional crystal growth
described by a forward-backward parabolic equation. The ill-posed
region of the equation describes the motion of corners on the surface.
We analyze a fourth order regularized version of this equation and
show that the dynamical behavior of the regularized corner can be
described by a traveling wave solution. The speed of the wave is found
by rigorous asymptotic analysis. The interaction between multiple
corners will also be presented together with numerical simulations.
This is joint work in progress with Fang Wan.
14:15
Numerical Approximation and BSDE representation for Switching Problems
Abstract
Hamadène and Jeanblanc provided a BSDE representation for the resolution of bi-dimensional continuous time optimal switching problems. For example, an energy producer faces the possibility to switch on or off a power plant depending on the current price of electricity and corresponding comodity. A BSDE representation via multidimensional reflected BSDEs for this type of problems in dimension larger than 2 has been derived by Hu and Tang as well as Hamadène and Zhang [2]. Keeping the same example in mind, one can imagine that the energy producer can use different electricity modes of production, and switch between them depending on the commodity prices. We propose here an alternative BSDE representation via the addition of constraints and artificial jumps. This allows in particular to reinterpret the solution of multidimensional reflected BSDEs in terms of one-dimensional constrained BSDEs with jumps. We provide and study numerical schemes for the approximation of these two type of BSDEs
14:00
Modelling Sterile Insect Techniques in variable Mosquito Populations
10:00
OCIAM meeting
Abstract
NO WORKSHOP - 09:45 coffee in DH Common Room for those attending the OCIAM Meeting
17:00
Elimination of wild ramification and local uniformization in arbitrary characteristic
Abstract
Elimination of wild ramification is used in the structure theory of valued function fields, with applications in areas such as local uniformization (i.e., local resolution of singularities) and the model theory of valued fields. I will give a survey on the role that Artin-Schreier extensions play in the elimination of wild ramification, and corresponding main theorems on the structure of valued function fields. I will show what these results tell us about local uniformization. I have shown that local uniformization is always possible after a separable extension of the function field of the algebraic variety (separable "alteration"). This was extended to the arithmetic case in joint work with Hagen Knaf. Recently, Michael Temkin has proved local uniformization by purely inseparable alteration.
Further, I will describe a classification of Artin-Schreier extensions with non-trivial defect. It can be used to improve one of the above mentioned main theorems ("Henselian Rationality"). This could be a key for a purely valuation theoretical proof of Temkin's result. On the other hand, the classification shows that separable alteration and purely inseparable alteration are just two ways to eliminate the critical defects. So the existence of these two seamingly "orthogonal" local uniformization results does not necessarily indicate that local uniformization without alteration is possible.
Old and new problems in the multiscale nonlinear dynamics of beams, strings, chains, rods, ropes, filaments and whips
A Primal-Dual Regularized Interior-Point Method for Convex Quadratic Programs
Abstract
Interior-point methods for linear and convex quadratic programming
require the solution of a sequence of symmetric indefinite linear
systems which are used to derive search directions. Safeguards are
typically required in order to handle free variables or rank-deficient
Jacobians. We propose a consistent framework and accompanying
theoretical justification for regularizing these linear systems. Our
approach is akin to the proximal method of multipliers and can be
interpreted as a simultaneous proximal-point regularization of the
primal and dual problems. The regularization is termed "exact" to
emphasize that, although the problems are regularized, the algorithm
recovers a solution of the original problem. Numerical results will be
presented. If time permits we will illustrate current research on a
matrix-free implementation.
This is joint work with Michael Friedlander, University of British Columbia, Canada
13:00
Girsanov's theorem, martingale representation and BSDE
Abstract
This talk I present a study of BSDEs with non-linear terms of quadratic growth by using Girsanov's theorem. In particular we are able to establish a non-linear version of the Cameron-Martin formula, which can be for example used to obtain gradient estimates for some non-linear parabolic equations.
Rotational Elasticity
Abstract
We consider a 3-dimensional elastic continuum whose material points
can experience no displacements, only rotations. This framework is a
special case of the Cosserat theory of elasticity. Rotations of
material points of the continuum are described mathematically by
attaching to each geometric point an orthonormal basis which gives a
field of orthonormal bases called the coframe. As the dynamical
variables (unknowns) of our theory we choose the coframe and a
density.
In the first part of the talk we write down the general dynamic
variational functional of our problem. In doing this we follow the
logic of classical linear elasticity with displacements replaced by
rotations and strain replaced by torsion. The corresponding
Euler-Lagrange equations turn out to be nonlinear, with the source
of this nonlinearity being purely geometric: unlike displacements,
rotations in 3D do not commute.
In the second part of the talk we present a class of explicit
solutions of our Euler-Lagrange equations. We call these solutions
plane waves. We identify two types of plane waves and calculate
their velocities.
In the third part of the talk we consider a particular case of our
theory when only one of the three rotational elastic moduli, that
corresponding to axial torsion, is nonzero. We examine this case in
detail and seek solutions which oscillate harmonically in time but
depend on the space coordinates in an arbitrary manner (this is a
far more general setting than with plane waves). We show [1] that
our second order nonlinear Euler-Lagrange equations are equivalent
to a pair of linear first order massless Dirac equations. The
crucial element of the proof is the observation that our Lagrangian
admits a factorisation.
[1] Olga Chervova and Dmitri Vassiliev, "The stationary Weyl
equation and Cosserat elasticity", preprint http://arxiv.org/abs/1001.4726
An introduction to the Ricci flow
Abstract
The aim of this talk is to get a feel for the Ricci flow. The Ricci flow was introduced by Hamilton in 1982 and was later used by Perelman to prove the Poincaré conjecture. We will introduce the notions of Ricci flow and Ricci soliton, giving simple examples in low dimension. We will also discuss briefly other types of geometric flows one can consider.
14:00
Normal subsystems of fusion systems
Abstract
There are two competing notions for a normal subsystem of a (saturated) fusion system. A recent theorem of mine shows how the two notions are related. In this talk I will discuss normal subsystems and their properties, and give some ideas on why this might be useful or interesting.
Hopf-Galois extensions and an exact sequence for H-Picard groups
Abstract
The topic of this talk is the representation theory of Hopf-Galois extensions. We consider the following questions.
Let H be a Hopf algebra, and A, B right H-comodule algebras. Assume that A and B are faithfully flat H-Galois extensions.
1. If A and B are Morita equivalent, does it follow that the subalgebras A^coH and B^coH of H-coinvariant elements are also Morita equivalent?
2. Conversely, if A^coH and B^coH are Morita equivalent, when does it follow that A and B are Morita equivalent?
As an application, we investigate H-Morita autoequivalences of the H-Galois extension A, introduce the concept of H-Picard group, and we establish an exact sequence linking the H-Picard group of A and
the Picard group of A^coH.(joint work with Stefaan Caenepeel)
Isotopy of Lagrangian submanifolds
Abstract
Lagrangian submanifolds are an important class of objects in symplectic geometry. They arise in diverse settings: as vanishing cycles in complex algebraic geometry, as invariant sets in integrable systems, as Heegaard tori in Heegaard-Floer theory and of course as "branes" in the A-model of mirror symmetry. We ask the difficult question: when are two Lagrangian submanifolds isotopic? Restricting to the simplest case of Lagrangian spheres in rational surfaces we will give examples where this question has a complete answer. We will also give some very pictorial examples (due to Seidel) illustrating how two Lagrangians can fail to be isotopic.
Infinitely many bifurcations for some nonlinear elliptic equations with rapidly growing nonlinearity.
15:45
15:45
Higher string topology
Abstract
The talk will begin with a brief account of the construction of string topology operations. I will point out some mysteries with the formulation of these operations, such as the role of (moduli) of surfaces, and pose some questions. The remainder of the talk will address these issues. In particular, I will sketch some ideas for a higher-dimensional version of string topology. For instance, (1) I will describe an E_{d+1} algebra structure on the (shifted) homology of the free mapping space H_*(Map(S^d,M^n)) and (2) I will outline how to obtain operations H_*(Map(P,M)) -> H_*(Map(Q,M)) indexes by a moduli space of zero-surgery data on a smooth d-manifold P with resulting surgered manifold Q.
Delzant and Potyagailo's hierarchical accessibility
Abstract
Take a group G and split it as the fundamental group of a graph of groups, then take the vertex groups and split them as fundamental groups of graphs of groups etc. If at some point you end up with a collection of unsplittable groups, then you have a hierarchy. Haken showed that for any 3-manifold M with an incompressible surface S, one can cut M along S and and then find other incompressible surfaces in M\S and cut again, and repeating this process one eventually obtains a collection of balls. Analogously, Delzant and Potyagailo showed that for any finitely presented group without 2-torsion and a certain sensible class E of subgroups of G, G admits a hierarchy where the edge groups of the splittings lie in E. I really like their proof and I will present it.
A priori estimates for the weak solutions to the Boltzmann equation with grazing collisions
Abstract
In this talk we consider the Boltzmann equation arising in gas dynamics with long-range interactions. Mathematically, it involves bilinear singular integral operators known as collision operators with non-cutoff collision kernels. As for the associated Cauchy problem, we develop a theory of weak solutions and present some of its a priori estimates related with physical quantities including the energy and moments.
On Fields over Fields
Abstract
Spectral analysis of the discrete Helmholtz operator preconditioned with a shifted Laplacian
Abstract
Shifted Laplace preconditioners have attracted considerable attention as
a technique to speed up convergence of iterative solution methods for the
Helmholtz equation. In the talk we present a comprehensive spectral
analysis of the discrete Helmholtz operator preconditioned with a shifted
Laplacian. Our analysis is valid under general conditions. The propagating
medium can be heterogeneous, and the analysis also holds for different types
of damping, including a radiation condition for the boundary of the computational
domain. By combining the results of the spectral analysis of the
preconditioned Helmholtz operator with an upper bound on the GMRES-residual
norm we are able to derive an optimal value for the shift, and to
explain the mesh-depency of the convergence of GMRES preconditioned
with a shifted Laplacian. We will illustrate our results with a seismic test
problem.
Joint work with: Yogi Erlangga (University of British Columbia) and Kees Vuik (TU Delft)
11:45
The Notion of a Stabilizer
Abstract
(IN: LADY MARGARET HALL)
As part of the Conference on Geometric Model Theory in honour of Professor Boris Zilber
09:00
One Two Three
Abstract
( IN: LADY MARGARET HALL)
As part of the Conference on Geometric Model Theory in honour of Professor Boris Zilber
14:45
14:30
Analytical and numerical solution of sedimentary delta formation
14:00
11:40
16:00
Valuation of P-adic L-functions and Multiplicative Euler Characterisitcs in Etale Cohomology
09:30