Past

I will explain recent joint work with Matthew Emerton and David Savitt, in which we relate the geometry of various tamely potentially Barsotti--Tate deformation rings for two-dimensional Galois representations to the integral structure of the cohomology of Shimura curves. As a consequence, we establish some conjectures of Breuil regarding this integral structure.
The key technique is the Taylor--–Wiles--–Kisin patching argument, which,when combined with a new, geometric perspective on the Breuil–--Mezard conjecture, forges a tight link between the structure of cohomology (a global automorphic invariant) and local deformation rings (local Galois-theoretic invariants).

Composite membranes comprised of an ultra-thin coating film formed over a porous support membrane are the basis for state-of-the-art reverse osmosis (RO) and nanofiltration (NF) membranes, offering the possibility to independently optimize the support membrane and the coating film. However, limited information exists on transport through such composite membrane structures. Numerical calculations have been carried out in order to probe the impacts of the support membrane skin-layer pore size and porosity, support membrane bulk micro-porosity, and coating film thickness and morphology (i.e. surface roughness) on solvent and solute transport through composite membranes. Results suggest that the flux and rejection of a composite membrane may be fine-tuned, by adjusting support membrane skin layer porosity and pore size, independent of the properties of the coating film. Further, the water flux over the membrane surface is unevenly distributed, creating local ‘hot spots’ of high flux that may govern initial stages of membrane fouling and scaling. The analysis provides important insight on how the non-trivial interaction of support properties and film roughness may result in widely varying transport properties of the composite structure. In particular, the simulations reveal inherent trade-offs between flux, rejection and fouling propensity (the latter due to ‘hot spots’), which are purely consequences of geometrical factors, irrespective of materials chemistry.

I will be concerned with initial-boundary value problems for systems of conservation laws in one space variable. First, I will go over some of the most relevant features of these problems. In particular, I will stress that different viscous approximation lead, in general, to different limits.

Next, I will discuss possible ways of characterizing the limit of a given viscous approximation. Also, I will establish a uniqueness criterion that allows to conclude that the limit of a self-similar approximation introduced by Dafermos et al. actually coincide with the limit of the physical viscous approximation. Finally, if time allows I will mention consequences on the design of numerical schemes. The talk will be based on joint works with S. Bianchini, C. Christoforou and S. Mishra.

• Alfonso Bueno - Recent advances in mathematical modelling of cardiac tissue: A fractional step forward
• Matt Moore - Oblique water entry
• Matt Hennessy - Mathematical problems relating to organic solar cell production

Russian mathematician V.I.Arnold conjectured that convex, homogeneous bodies with less than four equilibria (also called mono-monostatic) may exist. Not only did his conjecture turn out to be true, the newly discovered objects show various interesting features. Our goal is to give an overview of these findings as well as to present some new results. We will point out that mono-monostatic bodies are neither flat, nor thin, they are not similar to typical objects with more equilibria and they are hard to approximate by polyhedra. Despite these "negative" traits, there seems to be strong indication that these forms appear in Nature due to their special mechanical properties.

Topological defects (TDs) are unavoidable consequence of continuous symmetry breaking phase transitions. They exhibit several universal features and often span apparently completely different systems. Particularly convenient testing ground to study basic physics of TDs are liquid crystals (LCs) due to their softness, liquid character and optical anisotropy. In the lecture I will present our recent theoretical studies of TDs in nematic LCs, which are of interest also to other branches of physics.

I will first focus on coarsening dynamics of TDs following the isotropic-nematic phase transition. Among others we have tested the validity of the Kibble-Zurek [1,2] prediction on the size of the so called protodomains, which was originally derived to estimate density of TDs as a function of inflation time in the early universe. Next I will consider nematic LC shells [3]. These systems are of interest because they could pave path to mm sized scaled crystals exhibiting different symmetries. Particular attention will be paid to curvature induced unbinding of pairs of topological defects. This process might play important role in membrane fission processes.  

[1] W.H. Zurek, Nature 317, 505 (1985).

[2] Z. Bradac et al.,  J.Chem.Phys 135, 024506 (2011)

[3] S. Kralj et al.,  Soft Matter 7, 670 (2011); 8, 2460  (2012).

Telomeres, non-coding terminal structures of DNA strands, consist of repetitive long tandem repeats of a specific length. An absence of an enzyme, telomerase, in certain cellular structures requires an alternative telomerase-independent pathway for telomeric sequence length regulation. Besides linear telomeres other configurations such as telomeric circles and telomeric loops were experimentally observed. They are suspected to play an important role in a universal mechanism for stabilization of the ends of linear DNA that possibly dates back to pre-telomerase ages. We propose a mathematical model that captures biophysical interactions of various telomeric structures on a short time scale and that is able to reproduce experimental measurements in mtDNA of yeast. Moreover, the model opens up a couple of interesting mathematical problems such as validity of a quasi-steady state approximation and dynamic properties of discrete coagulation-fragmentation systems. We also identify and estimate key factors influencing the length distribution of telomeric circles, loops and strand invasions using numerical simulations.

For a given positive measure on a fixed domain, Gaussian quadrature routines can be defined via their property of integrating an arbitrary polynomial of degree $2n+1$ exactly using only $n+1$ quadrature nodes. In the special case of Gauss--Jacobi quadrature, this means that $$\int_{-1}^1 (1+x)^\alpha(1-x)^\beta f(x) dx = \sum_{j=0}^{n} w_j f(x_j), \quad \alpha, \beta > -1, $$ whenever $f(x)$ is a polynomial of degree at most $2n+1$. When $f$ is not a polynomial, but a function analytic in a neighbourhood of $[-1,1]$, the above is not an equality but an approximation that converges exponentially fast as $n$ is increased.

An undergraduate mathematician taking a numerical analysis course could tell you that the nodes $x_j$ are roots of the Jacobi polynomial $P^{\alpha,\beta}_{n+1}(x)$, the degree $n+1$ polynomial orthogonal with respect to the given weight, and that the quadrature weight at each node is related to the derivative $P'^{\alpha,\beta}_{n+1}(x_j)$. However, these values are not generally known in closed form, and we must compute them numerically... but how?

Traditional approaches involve applying the recurrence relation satisfied by the orthogonal polynomials, or solving the Jacobi matrix eigenvalue problem in the algorithm of Golub and Welsch, but these methods are inherently limited by a minimal complexity of $O(n^2)$. The current state-of-the-art is the $O(n)$ algorithm developed by Glasier, Liu, and Rokhlin, which hops from root to root using a Newton iteration evaluated with a Taylor series defined by the ODE satisfied by $P^{\alpha,\beta}_{n+1}$.

We propose an alternative approach, whereby the function and derivative evaluations required in the Newton iteration are computed independently in $O(1)$ operations per point using certain well-known asymptotic expansions. We shall review these expansions, outline the new algorithm, and demonstrate improvements in both accuracy and efficiency. 

Organisers: Hilary Priestley, Drew Moshier and Leo Cabrer.

This will be devoted to the applications of dualities to logic and algebra, focusing on general techniques. Thus it will seek to complement the specialised coverage in meetings devoted to, for example, modal logic, residuated structures and many-valued logics, or coalgebras. The featured topics for the Workshop will be drawn from completions of ordered structures, and applications; admissible rules, unification theory, interpolation and amalgamation; aspects of many-valued and substructural logics and ordered algebraic structures. Keynote speakers will be Leo Cabrer and Mai Gehrke.

The energy of a deformed crystal is calculated in the context of a central force lattice model in two dimensions. When the crystal shape is a lattice polygon, it is shown that the energy equals the bulk elastic energy, plus the boundary integral of a surface energy density, plus the sum over the vertices of a corner energy function. This is an exact result when the interatomic potential has finite range; for an infinite-range potential it is asymptotically valid as the lattice parameter tends to zero. The surface energy density is obtained explicitly as a function of the deformation gradient and boundary normal. The corner energy is found as an explicit function of the deformation gradient and the normals of the two facets meeting at the corner. A new bond counting approach is used, which reduces the problem to certain lattice point problems of number theory. The approach is then extended to more general convex regions with possibly curved boundary. The resulting surface energy density depends on the unit normal in a striking way. It is continuous at irrational directions, discontinuous at rational ones and nowhere differentiable. The method also yields an explicit interfacial energy for twin and phase boundaries.

\[

%\large

We study nonnegative radial solutions to the problem

\begin{equation*}

\left\{

\begin{split}

-\Delta u = \lambda K(\left|x \right|) f(u), \quad x \in \Omega

\\u = 0 \quad \qquad \quad \qquad \mbox{if } \left|x \right| = r_0

\\u \rightarrow0 \quad \qquad \quad \qquad \mbox{as } \left|x \right|\rightarrow\infty,

\end{split} \right.

\end{equation*}

where $\lambda$ is a positive parameter, $\Delta u=\mbox{div} \big(\nabla u\big)$ is the Laplacian of $u$,

$\Omega=\{x\in\ \mathbb{R}^{n}; n \textgreater 2, \left|x \right| \textgreater r_0\}$ and $K$ belongs to a class of functions such that $\lim_{r\rightarrow \infty}K(r)=0$. For classes of nonlinearities $f$ that are negative at the origin and sublinear at $\infty$ we discuss existence and uniqueness results when $\lambda$ is large.

\]

• James Herterich         -           Mathematical modelling of water purification
• Paul Roberts               -           Mathematical models of retinal oxygen distribution
• Stephen O'Keeffe       -           Mathematical modelling of growth and stability in biological structures
• Andrey Melnik             -           Dynamics of anisotropic remodelling in elastic tissues

We consider two independent random sequences of length n.
We consider optimal alignments according to a scoring function S.
We show that when the scoring function S is chosen at random
then with probability 1, the frequency of the aligned letter pairs
converges to a unique distribution as n goes to infinity. We also show
some concentration of measure phenomena.

• Yujia Chen - Solving Surface Poisson's Equation via the Closest Point Method
• Alex Lewis - Modelling liquid crystal devices
• Georgina Lang - Modelling of Brain Tissue Swelling

Chief Organiser: Shari Levine.  Main speakers: Alexander Arhangel'skii, Alan Dow, Aisling McCluskey, Jan van Mill, Frank Tall, Vladimir Tkachuk

Contact for further information: @email

Respiration is a redox reaction in which oxidation of a substrate (often organic) is coupled to the reduction of a terminal electron acceptor (TEA) such as oxygen. Iron oxides in various mineral forms are abundant in sediments and sedimentary rocks, and many subsurface microbes have the ability to respire using Fe(III) as the TEA in anoxic conditions. This process is environmentally important in the degradation of organic substrates and in the redox-cycling of iron. But low mineral solubility limits the bioavailability of Fe(III), which microbes access primarily through reductive dissolution. For aqueous nutrients, expressions for microbial growth and nutrient uptake rates are standardly based on Monod kinetics. We address the question of what equivalent description is appropriate when solid phase Fe(III) is the electron acceptor.

Given a diffusion in R^n, we prove a small-noise expansion for its density. Our proof relies on the Laplace method on Wiener space and stochastic Taylor expansions in the spirit of Benarous-Bismut. Our result applies (i) to small-time asymptotics and (ii) to the tails of the distribution and (iii) to small volatility of volatility.

We shall study applications of this result to stochastic volatility models, recovering the Berestycki- Busca-Florent formula (using (i)), the Gulisashvili-Stein expansion (from (ii)) and Lewis' expansions (using (iii)).

This is a joint work with J.D. Deuschel (TU Berlin), P. Friz (TU Berlin) and S. Violante (Imperial College London).

A pseudofinite field is a perfect pseudo-algebraically closed (PAC) field which

has $\hat{\mathbb{Z}}$ as absolute Galois group. Pseudofinite fields exists and they can

be realised as ultraproducts of finite fields. A group $G$ is geometrically

represented in a theory $T$ if there are modles $M_0\prec M$ of $T$,

substructures $A,B$ of $M$, $B\subset acl(A)$, such that $M_0\le A\le B\le M$

and $Aut(B/A)$ is isomorphic to $G$. Let $T$ be a complete theory of

pseudofinite fields. We show that, geometric representation of a group whose order

is divisibly by $p$ in $T$ heavily depends on the presence of $p^n$'th roots of unity

in models of $T$. As a consequence of this, we show that, for almost all

completions of the theory of pseudofinite fields, over a substructure $A$, algebraic

closure agrees with definable closure, if $A$ contains the relative algebraic closure

of the prime field. This is joint work with Ehud Hrushovski.

The computational analysis of a mathematical model describing a complex system is often based on the following roadmap: first, an experiment is conceived, in which the measured data are (either directly or indirectly) related to the input data of the model equations; second, such equations are computationally solved to provide iconographic reconstructions of the unknown physical or physiological parameters of the system; third, the reconstructed images are utilized to validate the model or to inspire appropriate improvements. This talk will adopt such framework to investigate three applied problems, respectively in solar physics, neuroscience and physiology. The solar physics problem is concerned with the exploitation of hard X-ray data for the comprehension of energy transport mechanisms in solar flares. The neuroscientific problem is the one to model visual recognition in humans with the help of a magnetocencephalography experiment. Finally, the physiological problem investigates the kinetics of the kidney-bladder system by means of nuclear data.

While a general framework of approximating the solution to Hamilton-Jacobi-Bellman (HJB) equations by difference methods is well established, and efficient numerical algorithms are available for one-dimensional problems, much less is known in the multi-dimensional case. One difficulty is the monotone approximation of cross-derivatives, which guarantees convergence to the viscosity solution. We propose a scheme combining piecewise freezing of the policies in time with a suitable spatial discretisation to establish convergence for a wide class of equations, and give numerical illustrations for a diffusion equation with uncertain parameters. These equations arise, for instance, in the valuation of financial derivatives under model uncertainty.

This is joint work with Peter Forsyth.

A method of defining non-equilibrium entropy for a chaotic dynamical system is proposed which, unlike the usual method based on Boltzmann's principle $S = k\log W$, does not involve the concept of a macroscopic state. The idea is illustrated using an example based on Arnold's `cat' map. The example also demonstrates that it is possible to have irreversible behaviour, involving a large increase of entropy, in a chaotic system with only two degrees of freedom.

This is intended as an introductory talk about one of the most

important (and most geometric) aspects of Geometric Group Theory. No

prior knowledge of any maths will be assumed.

We present some recent results concerning domain wall motion in one-dimensional nanowires, including the existence, velocity and stability of travelling-wave and precessing solutions.  We consider the case of unixial anisotropy, characteristic of wires with symmetrical (e.g., circular) cross-section, as opposed to strongly anisotropic geometries (films and strips) that have received greater attention.  This is joint work with Arseni Goussev and Valeriy Slastikov.

Organisers: Hilary Priestley, Drew Moshier and Leo Cabrer.

This will be dedicated principally to extensions of duality theory beyond zero-dimensional structures and to its application in novel settings. Topics that are likely to feature include duality for bilattice-based structures and associated semantics; extensions to compact Hausdorff spaces, bitopological duality, and duality for continuous data; applications to coalgebraic logic. We shall be seeking two-way interaction between those focused on a particular application and those who are seeking to extend the theory. Keynote speakers will be Mike Mislove and Drew Moshier. Samson Abramsky will be away from Oxford fromJune 12, but we are grateful for his offer to give a talk on June 11. We are also pleased to announce that, through the good offices of Georg Gottlob (Oxford Department of Computer Science), we are able to include within W1 a tutorial lecture on the applications of bilattice semantics to computer science; this will be given by Ofer Arieli.

 Nuclear fusion offers the prospect of abundant clean energy production, but the physical and engineering challenges are very great. In nuclear fusion reactors, the fuel is in the form of a plasma (charged gas) which is confined at high temperature and density using a toroidal magnetic field. This configuration is susceptible to turbulence, which transports heat out of the plasma and prevents fusion. It is believed that rotating the plasma suppresses turbulence, but experiments are expensive and even modest numerical simulation requires hundreds of thousands of CPU hours. We present a numerical technique for one of the five phase-space dimensions that both improves the accuracy of the calculation and greatly reduces the resolution required.

This is joint work with Dani Wise and builds on his earlier

work. Let M be a compact oriented irreducible 3-manifold which is neither a

graph manifold nor a hyperbolic manifold. We prove that the fundamental

group of M is virtually special. This means that it virtually embeds in a

right angled Artin group, and is in particular linear over Z.

The Schramm-Loewner evolution (SLE(\kappa)) is a family of random fractal curves that arise in a natural way as scaling limits of interfaces in critical models in statistical physics. The SLE curves are constructed by solving the Loewner differential equation driven by a scaled Brownian motion. We will give an overview of some of the almost sure properties of SLE curves, concentrating in particular on properties that can be derived by studying the the geometry of growing curve locally at the tip. We will discuss a multifractual spectrum of harmonic measure at the tip, regularity in the capacity parameterization, and continuity of the curves as the \kappa-parameter is varied while the driving Brownian motion sample is kept fixed.

This is based on joint work with Greg Lawler, and with Steffen Rohde and Carto Wong.

An intriguing, and largely open, question in mathematical fluid dynamics is whether solutions of the Navier-Stokes equations converge in some sense to a solution of the Euler equations in the zero viscosity limit. In fact this convergence could conceivably fail due to oscillations and concentrations occuring in the sequence.

In the late 1980s, R. DiPerna and A. Majda extended the classical concept of Young measure to obtain a notion of measure-valued solution of the Euler equations, which records precisely these oscillation and concentration effects. In this talk I will present a result recently obtained in joint work with L. Székelyhidi, which states that any such measure-valued solution is generated by a sequence of distributional solutions of the Euler equations.

The result is interesting from two different viewpoints: On the one hand, it emphasizes the huge flexibility of the concept of weak solution for Euler; on the other hand, it provides an example of a characterization theorem for Young measures in the tradition of D. Kinderlehrer and P. Pedregal where the differential constraint on the generating sequence does not satisfy the constant rank condition.

There are two main statistical mechanical models of ferromagnetism: the simpler and better-understood Ising model, and the more realistic and more challenging classical Heisenberg model, where the spins are in the 2-sphere instead of in {-1,+1}. In dimensions one and two, the classical Heisenberg model with nearest-neighbor interactions has no phase transition, but in three dimensions it has been intractable.

To shed some light on the qualitative behavior of the 3D Heisenberg model, we use the versatile tools of mean-field theory and Stein's method in recent work with Elizabeth Meckes, studying the Heisenberg model on a complete graph with the number of vertices going to infinity. Our results include detailed descriptions of the magnetization, the empirical spin distribution, the free energy, and a second-order phase transition.

A map betweem metric spaces is a bilipschitz homeomorphism if it

is Lipschitz and has a Lipschitz inverse; a map is a bilipschitz embedding

if it is a bilipschitz homeomorphism onto its image. Given metric spaces

X and Y, one may ask if there is a bilipschitz embedding X--->Y, and if

so, one may try to find an embedding with minimal distortion, or at least

estimate the best bilipschitz constant. Such bilipschitz embedding

problems arise in various areas of mathematics, including geometric group

theory, Banach space geometry, and geometric analysis; in the last 10

years they have also attracted a lot of attention in theoretical computer

science.

The lecture will be a survey bilipschitz embedding in Banach spaces from

the viewpoint of geometric analysis.

Abstract: In this talk, I will discuss the proportionality between tree amplitudes and the ultraviolet divergences in their one-loop corrections in Yang-Mills and (N < 4) Super Yang-Mills theories in four-dimensions. From the point of view of local perturbative quantum field theory, i.e. Feynman diagrams, this proportionality is straightforward: ultraviolet divergences at loop-level are absorbed into coefficients of local operators/interaction vertices in the original tree-amplitude. Ultraviolet divergences in loop amplitudes are also calculable through on-shell methods. These methods ensure manifest gauge-invariance, even at loop-level (no ghosts), at the expense of manifest locality. From an on-shell perspective, the proportionality between the ultraviolet divergences the tree amplitudes is thus not guaranteed. I describe systematic structures which ensure proportionality, and their possible connections to other recent developments in the field.

• Savina Joseph - Current generation in solar cells
• Shengxin Zhu - Spectral distribution, smoothing effects and smoothness matching for radial basis functions
• Ingrid von Glehn - Solving surface PDEs with the closest point method

The section conjecture of Grothendieck's anabelian geometry speculates about a description of the set of rational

points of a hyperbolic curve over a number field entirely in terms of profinite groups and Galois theory.

In the talk we will discuss local to global aspects of the conjecture, and what can be achieved when sections with

additional group theoretic properties are considered.

Complex networks have been used to model almost any

real-world complex systems. An especially important

issue regards how to related their structure and dynamics,

which contributes not only for the better understanding of

such systems, but also to the prediction of important

dynamical properties from specific topological features.

In this talk I revise related research developed recently

in my group. Particularly attention is given to the concept

of accessibility, a new measurement integrating topology

and dynamics, and the relationship between frequency of

visits and node degree in directed modular complex

networks. Analytical results are provided that allow accurate

prediction of correlations between structure and dynamics

in systems underlain by directed diffusion. The methodology

is illustrated with respect to the macaque cortical network.