Past

Molecular Dynamics Simulations are a tool to study the behaviour

of atomic-scale systems. The simulations themselves solve the

equations of motion for hundreds to millions of particles over

thousands to billions of time steps. Due to the size of the

problems studied, such simulations are usually carried out on

large clusters or special-purpose hardware.

At a first glance, there is nothing much of interest for a

Numerical Analyst: the equations of motion are simple, the

integrators are of low order and the computational aspects seem

to focus on hardware or ever larger and faster computer

clusters.

The field, however, having been ploughed mainly by domain

scientists (e.g. Chemists, Biologists, Material Scientists) and

a few Computer Scientists, is a goldmine for interesting

computational problems which have been solved either badly or

not at all. These problems, although domain specific, require

sufficient mathematical and computational skill to make finding

a good solution potentially interesting for Numerical Analysts.

The proper solution of such problems can result in speed-ups

beyond what can be achieved by pushing the envelope on Moore's

Law.

In this talk I will present three examples where problems

interesting to Numerical Analysts arise. For the first two

problems, Constraint Resolution Algorithms and Interpolated

Potential Functions, I will present some of my own results. For

the third problem, using interpolations to efficiently compute

long-range potentials, I will only present some observations and

ideas, as this will be the main focus of my research in Oxford

and therefore no results are available yet.

According to the Ernst-Geroch reduction, in an axially symmetric stationary electrovac spacetime, the Einstein-Maxwell equations reduce to a harmonic map problem with singular boundary data. I will discuss the “regularity” of the reduced harmonic maps near the boundary and its implication on the regularity of the corresponding spacetimes.

We describe John Stalling's method of studying finitely generated free groups via graphs and moves on graphs called folds. We will then discuss how the theory can be extended to study the automorphism group of a finitely generated free group.

The representation theory of groups is surrounded by deep and difficult conjectures. In this talk we will take a tour of (some of) these problems, including Alperin's weight conjecture, Broué's conjecture, and Puig's finiteness conjecture.

A classic problem in invariant theory, often referred to as Hilbert's 14th problem, asks, when a group acts on a finitely generated commutative algebra by algebra automorphisms, whether the ring of invariants is still finitely generated. I shall present joint work with W. van der Kallen treating the problem for a Chevalley group over an arbitrary base. Progress on the corresponding problem of finite generation for rational cohomology will be discussed.

Algebraic classification of the Weyl tensor is an important tool for solving the Einstein equation. I shall review the classification for spacetimes of dimension greater than four, and recent progress in using it to construct new exact solutions. The higher-dimensional generalization of the Goldberg-Sachs theorem will be discussed.

The talk starts with the observation that many well-known systems of diffusive type

can be written as Wasserstein gradient flows. The aim of the talk is

to understand _why_ this is the case. We give an answer that uses a

connection between diffusive PDE systems and systems of Brownian

particles, and we show how the Wasserstein metric arises in this

context. This is joint work with Johannes Zimmer, Nicolas Dirr, and Stefan Adams.

We prove that when a sequence of Lévy processes $X(n)$ or a normed sequence of random walks $S(n)$ converges a.s. on the Skorokhod space toward a Lévy process $X$, the sequence $L(n)$ of local times at the supremum of $X(n)$ converges uniformly on compact sets in probability toward the local time at the supremum of $X$. A consequence of this result is that the sequence of (quadrivariate) ladder processes (both ascending and

descending) converges jointly in law towards the ladder processes of $X$. As an application, we show that in general, the sequence $S(n)$ conditioned to stay positive converges weakly, jointly with its local time at the future minimum, towards the corresponding functional for the limiting process $X$. From this we deduce an invariance principle for the meander which extends known results for the case of attraction to a stable law.

Defaults in a credit portfolio of many obligors or in an economy populated with firms tend to occur in waves. This may simply reflect their sharing of common risk factors and/or manifest their systemic linkages via credit chains. One popular approach to characterizing defaults in a large pool of obligors is the Poisson intensity model coupled with stochastic covariates, or the Cox process for short. A constraining feature of such models is that defaults of different obligors are independent events after conditioning on the covariates, which makes them ill-suited for modeling clustered defaults. Although individual default intensities under such models can be high and correlated via the stochastic covariates, joint default rates will always be zero, because the joint default probabilities are in the order of the length of time squared or higher. In this paper, we develop a hierarchical intensity model with three layers of shocks -- common, group-specific and individual. When a common (or group-specific) shock occurs, all obligors (or group members) face individual default probabilities, determining whether they actually default. The joint default rates under this hierarchical structure can be high, and thus the model better captures clustered defaults. This hierarchical intensity model can be estimated using the maximum likelihood principle. A default signature plot is invented to complement the typical power curve analysis in default prediction. We implement the new model on the US corporate bankruptcy data and find it far superior to the standard intensity model both in terms of the likelihood ratio test and default signature plot.

We consider the three dimensional Navier-Stokes equations with a large initial data and

we prove the existence of a global smooth solution. The main feature of the initial data

is that it varies slowly in the vertical direction and has a norm which blows up as the

small parameter goes to zero. In the language of geometrical optics, this type of

initial data can be seen as the ``ill prepared" case. Using analytical-type estimates

and the special structure of the nonlinear term of the equation we obtain the existence

of a global smooth solution generated by this large initial data. This talk is based on a

work in collaboration with J.-Y. Chemin and I. Gallagher and on a joint work with Z.

Zhang.

The aim is to explore whether we can extend the work of PM Woodward first published many years ago, to see if we can extract more information than we do to date from our radar returns. A particular interest is in the information available for target recognition, which requires going beyond Woodward's assumption that the target has no internal structure.

CFD is an indispensible part of the design process for all major gas turbine components. The growth in the use of CFD from single-block structured mesh steady state solvers to highly resolved unstructured mesh unsteady solvers will be described, with examples of the design improvements that have been achieved. The European Commission has set stringent targets for the reduction of noise, emissions and fuel consumption to be achieved by 2020. The application of CFD to produce innovative designs to meet these targets will be described. The future direction of CFD towards whole engine simulations will also be discussed.

I will survey the theory of quasiHamiltonian spaces, a.k.a. group valued moment maps. In rough correspondence with historical development, I will first show how they emerge from the study of loop group representations, and then how they arise as a special case of "presymplectic realizations" in Dirac geometry.

In this talk I will introduce and study opers over a smooth projective curve X defined over a field of positive characteristic. I will describe a bijective correspondence between the set of stable vector bundles E over X such that the pull-back F^*(E) under the Frobenius

map F of X has maximal Harder-Narasimhan polygon and the set of opers having zero p-curvature. These sets turn out to be finite, which allows us to derive dimensions of certain Quot-schemes and certain loci of stable Frobenius-destabilized vector bundles over X.

Fix a positive integer $k$, and consider the class of all graphs which do not have $k+1$  vertex-disjoint cycles.  A classical result of Erdos and P\'{o}sa says that each such graph $G$ contains a blocker of size at most $f(k)$.  Here a {\em blocker} is a set $B$ of vertices such that $G-B$ has no cycles.

We give a minor extension of this result, and deduce that almost all such labelled graphs on vertex set $1,\ldots,n$ have a blocker of size $k$.  This yields an asymptotic counting formula for such graphs; and allows us to deduce further properties of a graph $R_n$ taken uniformly at random from the class: we see for example that the probability that $R_n$ is connected tends to a specified limit as $n \to \infty$.

There are corresponding results when we consider unlabelled graphs with few disjoint cycles. We consider also variants of the problem involving for example disjoint long cycles.

This is joint work with Valentas Kurauskas and Mihyun Kang.

Gupta et al. introduced a very general multi-commodity flow problem in which the cost of a given flow solution on a graph $G=(V,E)$ is calculated by first computing the link loads via a load-function l, that describes the load of a link as a function of the flow traversing the link, and then aggregating the individual link loads into a single number via an aggregation function.

We show the existence of an oblivious routing scheme with competitive ratio $O(\log n)$ and a lower bound of $\Omega(\log n/\logl\og n)$ for this model when the aggregation function agg is an $L_p$-norm.

Our results can also be viewed as a generalization of the work on approximating metrics by a distribution over dominating tree metrics and the work on minimum congestion oblivious. We provide a convex combination of trees such that routing according to the tree distribution approximately minimizes the $L_p$-norm of the link loads. The embedding techniques of Bartal and Fakcharoenphol et al. [FRT03] can be viewed as solving this problem in the $L_1$-norm while the result on congestion minmizing oblivious routing solves it for $L_\infty$. We give a single proof that shows the existence of a good tree-based oblivious routing for any $L_p$-norm.

Certain elastic structures and growing tissues (leaves, flowers or marine invertebrates) exhibit residual strain at free equilibria. We intend to study this phenomena through an elastic growth variational model. We will first discuss this model from a differential geometric point of view: the growth seems to change the intrinsic metric of the tissue to a new target non-flat metric. The non-vanishing curvature is the cause of the non-zero stress at equilibria.

We further discuss the scaling laws and $\Gamma$-limits of the introduced 3d functional on thin plates in the limit of vanishing thickness. Among others, given special forms of growth tensors, we rigorously derive the non-Euclidean versions of Kirchhoff and von Karman models for elastic non-Euclidean plates. Sobolev spaces of isometries and infinitesimal isometries of 2d Riemannian manifolds appear as the natural space of admissible mappings in this context. In particular, as a side result, we obtain an equivalent condition for existence of a $W^{2,2}$ isometric immersion of a given $2$d metric on a bounded domain into $\mathbb R3$.

Lattices in semisimple Lie groups have been studied from the point of view of number theory, algebraic groups, topology and geometry, and geometric group theory. The Fragestellung of one line of investigation is to what extent the properties of the lattice determine, and are determined by, the properties of the group. This talk reviews a number of results about lattices, and in particular looks at Mostow--Margulis rigidity.