OCCAM group meeting
Abstract
• Amy Smith presents: “Multiscale modelling of coronary blood flow derived from the microstructure”
• Laura Gallimore presents: “Modelling Cell Motility”
• Jean-Charles Seguis presents: “Coupling the membrane with the cytosol: a first encounter”
Climate: Methods for Model Comparison and Criticism in the context of data assimilation
Abstract
This will not be a normal workshop with a single scientist presenting an unsolved problem where mathematics may help. Instead it is more of a discussion meeting with a few speakers all interested in a single theme. So far we have:
Lenny Smith (LSE) on Using Empirically Inadequate Models to inform Your Subjective Probabilities: How might Solvency II inform climate change decisions?
Dan Rowlands (AOPP, Oxford) on "objective" climate forecasting;
Tim Palmer (ECMWF and AOPP, Oxford) on Constraining predictions of climate change using methods of data assimilation;
Chris Farmer (Oxford) about the problem of how to ascertain the error in the equations of a model when in the midst of probabilistic forecasting and prediction.
16:30
Models for dynamics of the human tear film
Abstract
We study dynamics from models for the human tear film in one and two dimensional domains.
The tear film is roughly a few microns thick over a domain on a centimeter scale; this separation of scales makes lubrication models desirable. Results on one-dimensional blinking domains are presented for multiple blink cycles. Results on two-dimensional stationary domains are presented for different boundary conditions. In all cases, the results are sensitive to the boundary conditions; this is intuitively satisfying since the tear film seems to be controlled primarily from the boundary and its motion. Quantitative comparison with in vivo measurement will be given in some cases. Some discussion of tear film properties will also be given, and results for non-Newtonian models will be given as available, as well as some future directions.
16:00
Entropy method for line-energies
Abstract
The aim of this talk is to analyze energy functionals concentrated on the jump set of 2D vector fields of unit length and of vanishing divergence.
The motivation of this study comes from thin-film micromagnetics where these functionals correspond to limiting wall-energies. The main issue consists in characterizing the wall-energy density (the cost function) so that the energy functional is lower semicontinuous (l.s.c.). The key point resides in the concept of entropies due to the scalar conservation law implied by our vector fields. Our main result identifies appropriate cost functions
associated to certain sets of entropies. In particular, certain power cost functions lead to l.s.c. energy functionals.
A second issue concerns the existence of minimizers of such energy functionals that we prove via a compactness result. A natural question is whether the viscosity solution is a minimizing configuration. We show that in general it is not the case for nonconvex domains.
However, the case of convex domains is still open. It is a joint work with Benoit Merlet, Ecole Polytechnique (Paris).
Nonlinear Eigenvalue Problems
Abstract
Nonlinear eigenvalue problem (NEP) is a class of eigenvalue problems where the matrix depends on the eigenvalue. We will first introduce some NEPs in real applications and some algorithms for general NEPs. Then we introduce our recent advances in NEPs, including second order Arnoldi algorithms for large scale quadratic eigenvalue problem (QEP), analysis and algorithms for symmetric eigenvalue problem with nonlinear rank-one updating, a new linearization for rational eigenvalue problem (REP).
Understanding the Second Mémoire of Évariste Galois: modern counterexamples to ancient writings
Abstract
The first part of Galois' Second Mémoire, less than three pages of manuscript written in 1830, is devoted to an amazing insight, far ahead of its time. Translated into modern mathematical language (and out of French), it is the theorem that a primitive soluble finite permutation group has prime-power degree. This, and Galois' ideas, and counterexamples to some of
them, will be my theme.
10:10
Normal forms for lattice polarized K3 surfaces and the Kuga-Satake Hodge Conjecture II (HoRSe seminar)
Establishing Complexity of Problems Parameterized Above Average
Abstract
In the Max Acyclic Subdigraph problem we are given a digraph $D$ and ask whether $D$ contains an acyclic subdigraph with at least $k$ arcs. The problem is NP-complete and it is easy to see that the problem is fixed-parameter tractable, i.e., there is an algorithm of running time $f(k)n$ for solving the problem, where $f$ is a computable function of $k$ only and $n=|V(D)|$. The last result follows from the fact that the average number of arcs in an acyclic subdigraph of $D$ is $m/2$, where $m$ is the number of arcs in $D$. Thus, it is natural to ask another question: does $D$ have an acyclic subdigraph with at least $m/2 +k$ arcs?
Mahajan, Raman and Sikdar (2006, 2009), and by Benny Chor (prior to 2006) asked whether this and other problems parameterized above the average are fixed-parameter tractable (the problems include Max $r$-SAT, Betweenness, and Max Lin). Most of there problems have been recently shown to be fixed-parameter tractable.
Methods involved in proving these results include probabilistic inequalities, harmonic analysis of real-valued
functions with boolean domain, linear algebra, and algorithmic-combinatorial arguments. Some new results obtained in this research are of potential interest for several areas of discrete mathematics and computer science. The examples include a new variant of the hypercontractive inequality and an association of Fourier expansions of real-valued functions with boolean domain with weighted systems of linear equations over $F^n_2$.
I’ll mention results obtained together with N. Alon, R. Crowston, M. Jones, E.J. Kim, M. Mnich, I.Z. Ruzsa, S. Szeider, and A. Yeo.
Normal forms for lattice polarized K3 surfaces and the Kuga-Satake Hodge Conjecture (HoRSe seminar)
Modelling cell migration in the mouse embryo
Abstract
The visceral endoderm (VE) is an epithelium of approximately 200 cells
encompassing the early post-implantation mouse embryo. At embryonic day
5.5, a subset of around 20 cells differentiate into morphologically
distinct tissue, known as the anterior visceral endoderm (AVE), and
migrate away from the distal tip, stopping abruptly at the future
anterior. This process is essential for ensuring the correct orientation
of the anterior-posterior axis, and patterning of the adjacent embryonic
tissue. However, the mechanisms driving this migration are not clearly
understood. Indeed it is unknown whether the position of the future
anterior is pre-determined, or defined by the movement of the migrating
cells. Recent experiments on the mouse embryo, carried out by Dr.
Shankar Srinivas (Department of Physiology, Anatomy and Genetics) have
revealed the presence of multicellular ‘rosettes’ during AVE migration.
We are developing a comprehensive vertex-based model of AVE migration.
In this formulation cells are treated as polygons, with forces applied
to their vertices. Starting with a simple 2D model, we are able to mimic
rosette formation by allowing close vertices to join together. We then
transfer to a more realistic geometry, and incorporate more features,
including cell growth, proliferation, and T1 transitions. The model is
currently being used to test various hypotheses in relation to AVE
migration, such as how the direction of migration is determined, what
causes migration to stop, and what role rosettes play in the process.
12:00
Characterization and Rigidity of the Kerr-Newman Solution
Abstract
A celebrated result in mathematical general relativity is the uniqueness of the Kerr(-Newman) black-holes as regular solutions to the stationary and axially-symmetric Einstein(-Maxwell) equations. The axial symmetry can be removed if one invokes Hawking's rigidity theorem. Hawking's theorem requires, however, real analyticity of the solution. A recent program of A. Ionescu and S. Klainerman seeks to remove the analyticity requirement in the vacuum case. They were able to show that any smooth extension of "Kerr data" prescribed on the horizon, satisfying the Einstein vacuum equations, must be Kerr, using a characterization of Kerr metric due to M. Mars. In this talk I will give a characterization for the Kerr-Newman metric, and extend the rigidity result to cover the electrovacuum case.
Global regular solutions to the Navier-Stokes equations in a cylinder with slip boundary conditions
Abstract
We consider the motion of a viscous incompressible fluid described by
the Navier-Stokes equations in a bounded cylinder with slip boundary
conditions. Assuming that $L_2$ norms of the derivative of the initial
velocity and the external force with respect to the variable along the
axis of the cylinder are sufficiently small we are able to prove long
time existence of regular solutions. By the regular solutions we mean
that velocity belongs to $W^{2,1}_2 (Dx(0,T))$ and gradient of pressure
to $L_2(Dx(0,T))$. To show global existence we prolong the local solution
with sufficiently large T step by step in time up to infinity. For this purpose
we need that $L_2(D)$ norms of the external force and derivative
of the external force in the direction along the axis of the cylinder
vanish with time exponentially.
Next we consider the inflow-outflow problem. We assume that the normal
component of velocity is nonvanishing on the parts of the boundary which
are perpendicular to the axis of the cylinder. We obtain the energy type
estimate by using the Hopf function. Next the existence of weak solutions is
proved.
15:45
On spaces of homomorphisms and spaces of representations
Abstract
The subject of this talk is the structure of the space of homomorphisms from a free abelian group to a Lie group G as well as quotients spaces given by the associated space of representations.
These spaces of representations admit the structure of a simplicial space at the heart of the work here.
Features of geometric realizations will be developed.
What is the fundamental group or the first homology group of the associated space in case G is a finite, discrete group ?
This deceptively elementary question as well as more global information given in this talk is based on joint work with A. Adem, E. Torres, and J. Gomez.
15:45
Hermite polynomial aliasing in Gaussian quadrature
Abstract
A representation of Hermite polynomials of degree 2n + 1, as sum of an element in the polynomial ideal generated by the roots of the Hermite polynomial of degree n and of a reminder, suggests a folding of multivariate polynomials over a finite set of points. From this, the expectation of some polynomial combinations of random variables normally distributed is computed. This is related to quadrature formulas and has strong links with designs of experiments.
This is joint work with G. Pistone
14:15
Ergodicity of Hypoelliptic SDEs Driven by Fractional Brownian Motion
Abstract
We demonstrate that stochastic differential equations (SDEs) driven by fractional Brownian motion with Hurst parameter H > 1/2 have similar ergodic properties as SDEs driven by standard Brownian motion. The focus in this article is on hypoelliptic systems satisfying H\"ormander's condition. We show that such systems satisfy a suitable version of the strong Feller property and we conclude that they admit a unique stationary solution that is physical in the sense that it does not "look into the future".
The main technical result required for the analysis is a bound on the moments of the inverse of the Malliavin covariance matrix, conditional on the past of the driving noise.