Introduction to tensor numerical methods in higher dimensions
Abstract
Tensor numerical methods provide the efficient separable representation of multivariate functions and operators discretized on large $n^{\otimes d}$-grids, providing a base for the solution of $d$-dimensional PDEs with linear complexity scaling in the dimension, $O(d n)$. Modern methods of separable approximation combine the canonical, Tucker, matrix product states (MPS) and tensor train (TT) low-parametric data formats.
\\
\\
The recent quantized-TT (QTT) approximation method is proven to provide the logarithmic data-compression on a wide class of functions and operators. Furthermore, QTT-approximation makes it possible to represent multi-dimensional steady-state and dynamical equations in quantized tensor spaces with the log-volume complexity scaling in the full-grid size, $O(d \log n)$, instead of $O(n^d)$.
\\
\\
We show how the grid-based tensor approximation in quantized tensor spaces applies to super-compressed representation of functions and operators (super-fast convolution and FFT, spectrally close preconditioners) as well to hard problems arising in electronic structure calculations, such as multi-dimensional convolution, and two-electron integrals factorization in the framework of Hartree-Fock calculations. The QTT method also applies to the time-dependent molecular Schr{\"o}dinger, Fokker-Planck and chemical master equations.
\\
\\
Numerical tests are presented indicating the efficiency of tensor methods in approximation of functions, operators and PDEs in many dimensions.
\\
\\
Stochastic Control Representations for Penalized Backward Stochastic Differential Equations
Abstract
In this talk, We show that both reflected BSDE and its associated penalized BSDE admit both optimal stopping representation and optimal control
representation. We also show that both multidimensional reflected BSDE and its associated multidimensional penalized BSDE admit optimal switching representation. The corresponding optimal stopping problems for penalized BSDE have the feature that one is only allowed to stop at Poisson arrival times.
12:00
Quadratic interaction functional and structure of solutions to hyperbolic conservation laws
Abstract
The proof of several properties of solutions of hyperbolic systems of conservation laws in one space dimension (existence, stability, regularity) depends on the existence of a decreasing functional, controlling the nonlinear interactions of waves. In a special case (genuinely nonlinear systems) the interaction functional is quadratic, while in the general case it is cubic. Several attempts to prove the existence of a a quadratic functional also in the most general case have been done. I will present the approach we follow in order to prove this result, an some of its implication we hope to exploit.
\\ \\
Work in collaboration with Stefano Modena.
16:00
Symbolic dynamics: taking another look at complex quadratic maps
Abstract
Complex dynamical systems have been very well studied in recent years, in particular since computer graphics now enable us to peer deep into structures such as the Mandlebrot set and Julia sets, which beautifully illustrate the intricate dynamical behaviour of these systems. Using new techniques from Symbolic Dynamics, we demonstrate previously unknown properties of a class of quadratic maps on their Julia sets.
10:30
Existence and construction of expanders (expanders are good for you)
A model for a protein oscillator in Myxococcus xanthus
Abstract
Cell polarity in the rod-shaped bacterium Myxococcus xanthus is crucial for the direction of movement of individual cells. Polarity is governed by a regulatory system characterized by a dynamic spatiotemporal oscillation of proteins between the opposite cell poles. A mathematical framework for a minimal macroscopic model is presented which produces self-sustained regular oscillations of the protein concentrations. The mathematical model is based on a reaction-diffusion PDE system and is independent of external triggers. Necessary conditions on the reaction terms leading to oscillating solutions are derived theoretically. Possible scenarios for protein interaction are numerically tested for robustness against parameter variation. Finally, possible extensions of the model will be addressed.
17:00
Relatively hyperbolic groups, mapping class groups and random walks
Abstract
I will discuss similarities and differences between the geometry of
relatively hyperbolic groups and that of mapping class groups.
I will then discuss results about random walks on such groups that can
be proven using their common geometric features, namely the facts that
generic elements of (non-trivial) relatively hyperbolic groups are
hyperbolic, generic elements in mapping class groups are pseudo-Anosovs
and random paths of length $n$ stay $O(\log(n))$-close to geodesics in
(non-trivial) relatively hyperbolic groups and
$O(\sqrt{n}\log(n))$-close to geodesics in mapping class groups.
Limit method in extremal combinatorics
Abstract
Razborov's flag algebras provide a formal system
for operating with asymptotic inequalities between subgraph densities,
allowing to do extensive "book-keeping" by a computer. This novel use
of computers led to progress on many old problems of extremal
combinatorics. In some cases, finer structural information can be
derived from a flag algebra proof by by using the Removal Lemma or
graph limits. This talk will overview this approach.
14:15
Ensemble-based occurrence probability of extreme events from TIGGE [TBC] POSTPONED TO NEXT TERM
Symmetry operators and conserved quantities for fields on Kerr
Self-gravitating elastic bodies
Abstract
Self-gravitating elastic bodies provide models for extended
objects in general relativity. I will discuss constructions of static
and rotating self-gravitating bodies, as well as recent results on the
initial value problem for self-gravitating elastic bodies.
The complexity of group presentations, manifolds, and the Andrews-Curtis conjecture
Abstract
Nonnegative local martingales, Novikow's and Kazamaki's criteria, and the distribution of explosion times
Abstract
I will give a new proof for the famous criteria by Novikov and Kazamaki, which provide sufficient conditions for the martingale property of a nonnegative local martingale. The proof is based on an extension theorem for probability measures that can be considered as a generalization of a Girsanov-type change of measure.
In the second part of my talk I will illustrate how a generalized Girsanov formula can be used to compute the distribution of the explosion time of a weak solution to a stochastic differential equation
Poisson random forests and coalescents in expanding populations.
Abstract
Let (V, ≥) be a finite, partially ordered set. Say a directed forest on V is a set of directed edges [x,y> with x ≤ y such that no vertex has indegree greater than one.
Thus for a finite measure μ on some partially ordered measurable space D we may define a Poisson random forest by choosing a set of vertices V according to a Poisson point process weighted by the number of directed forests on V, and selecting a directed forest uniformly.
We give a necessary and sufficient condition for such a process to exist and show that the process may be expressed as a multi-type branching process with type space D.
We build on this observation, together with a construction of the simple birth death process due to Kurtz and Rodrigues [2011] to develop a coalescent theory for rapidly expanding populations.
14:15
Fenchel-Nielsen coordinates from spectral networks
Abstract
16:30
Recent applications of and trends in model theory.
Abstract
There are many recent points of contact of model theory and other
parts of mathematics: o-minimality and Diophantine geometry, geometric group
theory, additive combinatorics, rigid geometry,... I will probably
emphasize long-standing themes around stability, Diophantine geometry, and
analogies between ODE's and bimeromorphic geometry.
Modelling of irreversible deformations near the tip of a crack in a porous domain containing oil and gas
Abstract
Thermomechanical processes observed in deformable solids under intensive dynamic or quasi-static loadings consist of coupled mechanical, thermal and fracturing stages. The fracturing processes involve formation, motion and interaction of defects in crystals, phase transitions, breaking of bonds between atoms, accumulation of micro-structural damages (pores, cracks), etc. Irreversible deformations, zones of adiabatic shear micro-fractures are caused by these processes. Dynamic fracturing is a complicated multistage process, which includes appearance, evolution and confluence of micro-defects and formation of embryonic micro-cracks and pores that can grow and lead to the breaking-up of bodies with formation of free surfaces. This results in a need to use more advanced mathematical and numerical techniques.
This talk presents modelling of irreversible deformation near the tip of a crack in a porous domain containing oil and gas during the hydraulic fracturing process. The governing equations for a porous domain containing oil and gas are based on constructing a mathematical model of thermo-visco-elasto-plastic media with micro-defects (micro-pores) filled with another phase (e.g., oil or/and gas). The micro-pores can change their size during the process of dynamical irreversible deformation. The existing pores can expand or collapse. The model was created by using fundamental thermodynamic principles and, therefore, it is a thermodynamically consistent model. All the processes (i.e., irreversible deformation, fracturing, micro-damaging, heat transfer) within a porous domain are strongly coupled. An explicit normalized-corrected meshless method is used to solve the resulting governing PDEs. The flexibility of the proposed technique allows efficient running using a great number of micro- and macro- fractures. The results are presented, discussed and future studies are outlined.
Modelling chronic diseases and their consequences into the future reliably and usefully
Abstract
We wish to discuss the role of Modelling in Health Care. While risk factor prevalences vary and change with time it is difficult to anticipate the change in disease incidence that will result without accurately modelling the epidemiology. When detailed study of the prevalence of obesity, tobacco and salt intake, for example, are studied clear patterns emerge that can be extrapolated into the future. These can give rise to estimated probability distributions of these risk factors across age, sex, ethnicity, social class groups etc into the future. Micro simulation of individuals from defined populations (eg England 2012) can then estimate disease incidence, prevalence, death, costs and quality of life. Thus future health and other needs can be estimated, and interventions on these risk factors can be simulated for their population effect. Health policy can be better determined by a realistic characterisation of public health. The Foresight microsimulation modelling of the National Heart Forum (UK Health Forum) will be described. We will emphasise some of the mathematical and statistical issues associated with so doing.
Multiplicity in difference geometry
Abstract
The study of difference algebraic geometry stems from the efforts of Macintyre and Hrushovski to
count the number of solutions to difference polynomial equations over fields with powers of Frobenius.
We propose a notion of multiplicity in the context of difference algebraic schemes and prove a first principle
of preservation of multiplicity. We shall also discuss how to formulate a suitable intersection theory of difference schemes.
Feedback and embryonic stem cell fate commitment
Abstract
Self-renewal and pluripotency of mouse embryonic stem (ES) cells are controlled by a complex transcriptional regulatory network (TRN) which is rich in positive feedback loops. A number of key components of this TRN, including Nanog, show strong temporal expression fluctuations at the single cell level, although the precise molecular basis for this variability remains unknown. In this talk I will discuss recent work which uses a genetic complementation strategy to investigate genome-wide mRNA expression changes during transient periods of Nanog down-regulation. Nanog removal triggers widespread changes in gene expression in ES cells. However, we found that significant early changes in gene expression were reversible upon re-induction of Nanog, indicating that ES cells initially adopt a flexible “primed” state. Nevertheless, these changes rapidly become consolidated irreversible fate decisions in the continued absence of Nanog. Using high-throughput single cell transcriptional profiling we observed that the early molecular changes are both stochastic and reversible at the single cell level. Since positive feedback commonly gives rise to phenotypic variability, we also sought to determine the role of feedback in regulating ES cell heterogeneity and commitment. Analysis of the structure of the ES cell TRN revealed that Nanog acts as a feedback “linchpin”: in its presence positive feedback loops are active and the extended TRN is self-sustaining; while in its absence feedback loops are weakened, the extended TRN is no longer self-sustaining and pluripotency is gradually lost until a critical “point-of-no-return” is reached. Consequently, fluctuations in Nanog expression levels transiently activate different sub-networks in the ES cell TRN, driving transitions between a (Nanog expressing) feedback-rich, robust, self-perpetuating pluripotent state and a (Nanog-diminished), feedback-depleted, differentiation-sensitive state. Taken together, our results indicate that Nanog- dependent feedback loops play a central role in controlling both early fate decisions at the single cell level and cell-cell variability in ES cell populations.
How frequently does the Hasse principle fail?
Abstract
Counter-examples to the Hasse principle are known for many families of geometrically rational varieties. We discuss how often such failures arise for Chatelet surfaces and certain higher-dimensional hypersurfaces. This is joint work with Regis de la Breteche.
Centers and traces of categorified affine Hecke algebras (or, some tricks with coherent complexes on the Steinberg variety)
Abstract
The bounded coherent dg-category on (suitable versions of) the Steinberg stack of a reductive group G is a categorification of the affine Hecke algebra in representation theory. We discuss how to describe the center and universal trace of this monoidal dg-category. Many of the techniques involved are very general, and the description makes use of the notion of "odd micro-support" of coherent complexes. This is joint work with Ben-Zvi and Nadler.
Supersymmetry and Morse Theory
Abstract
Morse theory gives an estimate of the dimensions of the cohomology groups of a manifold in terms of the critical points of a function.
One can do better and compute the cohomology in terms of this function using the so-called Witten complex.
Already implicit in work of Smale in the fifties, it was rediscovered by Witten in the eighties using techniques from (supersymmetric) quantum field theories.
I will explain Witten's (heuristic) arguments and describe the Witten complex.
Deflating characters of symmetric groups and Foulkes’ Conjecture
Abstract
The symmetric group S_{mn} acts naturally on the collection of set partitions of a set of size mn into n sets each of size m, and the resulting permutation character is the Foulkes character. These characters are the subject of the longstanding Foulkes Conjecture. In this talk, we define a deflation map which sends a character of the symmetric group S_{mn} to a character of S_n. The values of the images of the irreducible characters under this map are described combinatorially in a rule which generalises two well-known combinatorial rules in the representation theory of symmetric groups, the Murnaghan-Nakayama formula and Young's rule. We use this in a new algorithm for computing irreducible constituents of Foulkes characters and verify Foulkes’ Conjecture in some new cases. This is joint work with Anton Evseev (Birmingham) and Mark Wildon (Royal Holloway).
Optimization meets Statistics: Fast global convergence for high-dimensional statistical recovery
Abstract
Many methods for solving high-dimensional statistical inverse problems are based on convex optimization problems formed by the weighted sum of a loss function with a norm-based regularizer.
\\
Particular examples include $\ell_1$-based methods for sparse vectors and matrices, nuclear norm for low-rank matrices, and various combinations thereof for matrix decomposition and robust PCA. In this talk, we describe an interesting connection between computational and statistical efficiency, in particular showing that the same conditions that guarantee that an estimator has good statistical error can also be used to certify fast convergence of first-order optimization methods up to statistical precision.
\\
\\
Joint work with Alekh Agarwahl and Sahand Negahban Pre-print (to appear in Annals of Statistics)
\\
http://www.eecs.berkeley.edu/~wainwrig/Papers/AgaNegWai12b_SparseOptFul…
Robust Portfolio Optimization under Heavy Tailed Returns
Abstract
We consider the problem of optimizing a portfolio of medium to low frequency
quant strategies under heavy tailed distributions. Approaching this problem by modelling
returns through mixture distributions, we derive robust and relative robust methodologies
and discuss conic optimization approaches to solving these models.
12:00
1D Burgers Turbulence as a model case for the Kolmogorov Theory
Abstract
The Kolmogorov 1941 theory (K41) is, in a way, the starting point for all
models of turbulence. In particular, K41 and corrections to it provide
estimates of small-scale quantities such as increments and energy spectrum
for a 3D turbulent flow. However, because of the well-known difficulties
involved in studying 3D turbulent flow, there are no rigorous results
confirming or infirming those predictions. Here, we consider a well-known
simplified model for 3D turbulence: Burgulence, or turbulence for the 1D
Burgers equation. In the space-periodic case with a stochastic white in
time and smooth in space forcing term, we give sharp estimates for
small-scale quantities such as increments and energy spectrum.
11:00
"Small rigid subsets of the reals"
Abstract
A topological space is called rigid if its only autohomeomorphism is the identity map. Using the Axiom of Choice it is easy to construct rigid subsets of the real line R, but sets constructed in this way always have size continuum. I will explore the question of whether it is possible to have rigid subsets of R that are small, meaning that their cardinality is smaller than that of the continuum. On the one hand, we will see that forcing can be used to produce models of ZFC in which such small rigid sets abound. On the other hand, I will introduce a combinatorial axiom that can be used to show the consistency with ZFC of the statement "CH fails but every rigid subset of R has size continuum". Only a working knowledge of basic set theory (roughly what one might remember from C1.2b) and topology will be assumed.
Self-similar groups
Abstract
Self-similarity is a fundamental idea in many areas of mathematics. In this talk I will explain how it has entered group theory and the links between self-similar groups and other areas of research. There will also be pretty pictures.
10:30
The McKay Correspondence
Abstract
I will give an introduction to The McKay Correspondence, relating the irreducible representations of a finite subgroup Γ ≤ SL2 (C), minimal resolutions of the orbit space C2 /Γ, and affine Dynkin diagrams.
Multiscale challenges and the hybrid method for stochastic simulation of biochemical systems
Abstract
Complex systems emerging from many biochemical applications often exhibit multiscale and multiphysics (MSMP) features: The systems incorporate a variety of physical processes or subsystems across a broad range of scales. A typical MSMP system may come across scales with macroscopic, mesoscopic and microscopic kinetics,
deterministic and stochastic dynamics, continuous and discrete state space, fastscale and slow-scale reactions, and species of both large and small populations. These complex features present great challenges in the modeling and simulation practice. The goal of our research is to develop innovative computational methods and rigorous fundamental theories to answer these challenges. In this talk we will start with introduction of basic stochastic simulation algorithms for biochemical systems and multiscale
features in the stochastic cell cycle model of budding yeast. With detailed analysis of these multiscale features, we will introduce recent progress on simulation algorithms and computational theories for multiscale stochastic systems, including tau-leaping methods, slow-scale SSA, and the hybrid method.
Hochschild invariants of matrix factorizations and classical singularity theory/derived critical loci
Bootstrap percolation on infinite trees
Abstract
While usual percolation concerns the study of the connected components of
random subgraphs of an infinite graph, bootstrap percolation is a type of
cellular automaton, acting on the vertices of a graph which are in one of
two states: `healthy' or `infected'. For any positive integer $r$, the
$r$-neighbour bootstrap process is the following update rule for the
states of vertices: infected vertices remain infected forever and each
healthy vertex with at least $r$ infected neighbours becomes itself
infected. These updates occur simultaneously and are repeated at discrete
time intervals. Percolation is said to occur if all vertices are
eventually infected.
As it is often difficult to determine precisely which configurations of
initially infected vertices percolate, one often considers a random case,
with each vertex infected independently with a fixed probability $p$. For
an infinite graph, of interest are the values of $p$ for which the
probability of percolation is positive. I will give some of the history
of this problem for regular trees and present some new results for
bootstrap percolation on certain classes of randomly generated trees:
Galton--Watson trees.
Functors between categories of IndCoherent complexes, deRham (Borel-Moore/co-) homology, and matrix factorizations
12:00
Mathematical modelling with fully anisotropic diffusion
Abstract
***** PLEASE NOTE THIS SEMINAR WILL TAKE PLACE ON TUESDAY 19TH FEBRUARY *****
With "fully anisotropic" I describe diffusion models of the form u_t=\nabla \nabla (D(x) u), where the diffusion tensor appears inside both derivatives. This model arises naturally in the modeling of brain tumor spread and wolf movement and other applications. Since this model does not satisfy a maximum principle, it can lead to interesting spatial pattern formation, even in the linear case. I will present a detailed derivation of this model and discuss its application to brain tumors and wolf movement. Furthermore, I will present an example where, in the linear case, the solution blows-up in infinite time; a quite surprising result for a linear parabolic equation (joint work with K.J. Painter and M. Winkler).
A Lagrangian approach for nonhomogeneous incompressible fluids
Abstract
In this talk we focus on the incompressible Navier–Stokes equations with variable
density. The aim is to prove existence and uniqueness results in the case of a discontinuous
initial density (typically we are interested in discontinuity along an interface).
In the first part of the talk, by making use of Fourier analysis techniques, we establish the existence of global-in-time unique solutions in a critical
functional framework, under some smallness condition over the initial data,
In the second part, we use another approach to avoid the smallness condition over the nonhomogeneity : as a matter of fact, one may consider any density bounded
and bounded away from zero and still get a unique solution. The velocity is required to have subcritical regularity, though.
In all the talk, the Lagrangian formulation for describing the flow plays a key role in the analysis.
A continuum of exponents for the rate of escape of random walks on groups
Abstract
Abstract: A central question in the theory of random walks on groups is how symmetries of the underlying space gives rise to structure and rigidity of the random walks. For example, for nilpotent groups, it is known that random walks have diffusive behavior, namely that the rate of escape, defined as the expected distance of the walk from the identity satisfies E|Xn|~=n^{1/2}. On nonamenable groups, on the other hand we have E|Xn| ~= n. (~= meaning upto (multiplicative) constants )
In this work, for every 3/4 <= \beta< 1 we construct a finitely generated group so that the expected distance of the simple random walk from its starting point after n steps is n^\beta (up to constants). This answers a question of Vershik, Naor and Peres. In other examples, the speed exponent can fluctuate between any two values in this interval.
Previous examples were only of exponents of the form 1-1/2^k or 1 , and were based on lamplighter (wreath product) constructions.(Other than the standard beta=1/2 and beta=1 known for a wide variety of groups) In this lecture we will describe how a variation of the lamplighter construction, namely the permutational wreath product, can be used to get precise bounds on the rate of escape in terms of return probabilities of the random walk on some Schreier graphs. We will then show how groups of automorphisms of rooted trees, related to automata groups , can be constructed and analyzed to get the desired rate of escape. This is joint work with Balint Virag of the University of Toronto. (Paper available at http://arxiv.org/abs/1203.6226)
No previous knowledge of random walks,automaton groups or wreath products is assumed.
Rough paths, controlled distributions, and nonlinear SPDEs
Abstract
Abstract: Hairer recently had the remarkable insight that Lyons' theory of rough paths can be used to make sense of nonlinear SPDEs that were previously ill-defined due to spatial irregularities. Since rough path theory deals with the integration of functions defined on the real line, the SPDEs studied by Hairer have a one-dimensional spatial index variable. I will show how to combine paraproducts, a notion from functional analysis, with ideas from the theory of controlled rough paths, in order to develop a formulation of rough path theory that works in any index dimension. As an application, I will present existence and uniqueness results for an SPDE with multidimensional spatial index set, for which previously it was not known how to describe solutions. No prior knowledge of rough paths or paraproducts is required for understanding the talk. This is joint work with Massimiliano Gubinelli and Peter Imkeller.
14:15