L3

Tensors, also known as multiway arrays, have become ubiquitous as representations for operators or as convenient schemes for storing data. Yet, when it comes to compressing these objects or analyzing the data stored in them, the tendency is to ``flatten” or ``matricize” the data and employ traditional linear algebraic tools, ignoring higher dimensional correlations/structure that could have been exploited. Impediments to the development of equivalent tensor-based approaches stem from the fact that familiar concepts, such as rank and orthogonal decomposition, have no straightforward analogues and/or lead to intractable computational problems for tensors of order three and higher.

In this talk, we will review some of the common tensor decompositions and discuss their theoretical and practical limitations. We then discuss a family of tensor algebras based on a new definition of tensor-tensor products. Unlike other tensor approaches, the framework we derive based around this tensor-tensor product allows us to generalize in a very elegant way all classical algorithms from linear algebra. Furthermore, under our framework, tensors can be decomposed in a natural (e.g. ‘matrix-mimetic’) way with provable approximation properties and with provable benefits over traditional matrix approximation. In addition to several examples from recent literature illustrating the advantages of our tensor-tensor product framework in practice, we highlight interesting open questions and directions for future research.

PDE Workshop in Stability Analysis for Nonlinear PDEs will be running Monday 15th - Friday 19th August.

Location: L3, AWB

Our goal was to bring together leading experts in the stability analysis of nonlinear partial differential equations across multi-scale applications. Some of the topics to be addressed include: 

• Stability analysis of shock wave patterns of reflections/diffraction.
• Stability analysis of vortex sheets, contact discontinuities, and other characteristic discontinuities for multidimensional hyperbolic systems of conservation laws.
• Stability analysis of particle to continuum limits including the quantifying asymptotic/mean-field/large-time limits for pairwise interactions and particle limits for general interactions among multi-agent systems
• Stability analysis of asymptotic limits with emphasis on the vanishing viscosity limit of solutions from multidimensional compressible viscous to inviscid flows with large initial data.

The size-Ramsey number $\hat{r}(H)$ of a graph $H$ is the smallest number of edges a (host) graph $G$ can have, such that for any red/blue coloring of $G$, there is a monochromatic copy of $H$ in $G$. Recently, Conlon, Nenadov and Trujić showed that if $H$ is a graph on $n$ vertices and maximum degree three, then $\hat{r}(H) = O(n^{8/5})$, improving upon the bound of $n^{5/3 + o(1)}$ by Kohayakawa, Rödl, Schacht and Szemerédi. In our paper, we show that $\hat{r}(H)\leq n^{3/2+o(1)}$. While the previously used host graphs were vanilla binomial random graphs, we prove our result by using a novel host graph construction.
We also discuss why our bound is a natural barrier for the existing methods.
This is joint work with Kalina Petrova.

Detectors in collider experiments are modeled by light-ray operators in Quantum Field Theory. For example, energy detectors are certain null integrals of the stress-energy tensor, localized at an angle on the celestial sphere, where they collect quanta that escape in their direction. In this talk, I will discuss a series of work developing a nonperturbative, convergent operator product expansion (OPE) for light-ray operators in Conformal Field Theories (CFTs). Objects appearing in the expansion are more general light-ray operators, whose matrix elements can be computed by the generalized Lorentzian inversion formula. An important application is to event shapes in collider physics, which correspond to correlation functions of light-ray operators within the state created by the incoming particles. I will discuss some applications of the light-ray OPE in CFT, and mention some extensions to QCD which make contact with measurements at the LHC. Talk based primarily on [1905.01311] and [2010.04726].

We develop a method to investigate the geometric quantisation of a hypertoric variety from an equivariant viewpoint, in analogy with the equivariant Verlinde for Higgs bundles. We do this by first using the residual circle action on a hypertoric variety to construct its symplectic cut, resulting in a compact cut space which is needed for localisation. We introduce the notion of a moment polyptych associated to a hypertoric variety and prove that the necessary isotropy data can be read off from it. Finally, the equivariant Hirzebruch-Riemann-Roch formula is applied to the cut spaces and expresses the dimension of the equivariant quantisation space as a finite sum over the fixed-points. This is joint work with Johan Martens.

Motivated by the problem of reconstructing the electron density of a molecule from pulsed X-ray diffraction images (about 10e+9 per reconstruction), we develop a framework for analyzing the convergence to invariant measures of random fixed point iterations built from mappings that, while expansive, nevertheless possess attractive fixed points.  Building on techniques that we have established for determining rates of convergence of numerical methods for inconsistent nonconvex
feasibility, we lift the relevant regularities to the setting of probability spaces to arrive at a convergence analysis for noncontractive Markov operators.  This approach has many other applications, for instance the analysis of distributed randomized algorithms.
We illustrate the approach on the problem of solving linear systems with finite precision arithmetic.

 

Noether’s theorem plays a fundamental role in modern physics by relating symmetries of a Lagrangian to conserved quantities of the Euler-Lagrange equations. In ideal fluid dynamics, the theorem relates the particle labeling symmetry to a Kelvin circulation law. Circulation is conserved for incompressible flows and, otherwise, is generated by advected variables through the momentum map due to a broken symmetry. We will introduce variational principles for fluid dynamics that constrain advection to be the sum of a smooth and geometric rough-in-time vector field. The corresponding rough Euler-Poincare equations satisfy a Kelvin circulation theorem and lead to a natural framework to develop parsimonious non-Markovian parameterizations of subgrid-scale dynamics.

Given a variety defined over a field of characteristic zero and an algebraically integrable foliation of corank less than or equal to two, we show the existence of a categorical quotient, defined on the non-empty open subset of algebraically smooth points, through which every invariant morphism factors uniquely. Some applications to quotients by connected groups will be discussed.
 

Combining the classical theory of optimal transport with modern operator splitting techniques, I will present a new numerical method for nonlinear, nonlocal partial differential equations, arising in models of porous media,materials science, and biological swarming. Using the JKO scheme, along with the Benamou-Brenier dynamical characterization of the Wasserstein distance, we reduce computing the solution of these evolutionary PDEs to solving a sequence of fully discrete minimization problems, with strictly convex objective function and linear constraint. We compute the minimizer of these fully discrete problems by applying a recent, provably convergent primal dual splitting scheme for three operators. By leveraging the PDE’s underlying variational structure, ourmethod overcomes traditional stability issues arising from the strong nonlinearity and degeneracy, and it is also naturally positivity preserving and entropy decreasing. Furthermore, by transforming the traditional linear equality constraint, as has appeared in previous work, into a linear inequality constraint, our method converges in fewer iterations without sacrificing any accuracy. We prove that minimizers of the fully discrete problem converge to minimizers of the continuum JKO problem as the discretization is refined, and in the process, we recover convergence results for existing numerical methods for computing Wasserstein geodesics. Simulations of nonlinear PDEs and Wasserstein geodesics in one and two dimensions that illustrate the key properties of our numerical method will be shown.

Directed graphs are a model for various phenomena in the
sciences. In topological data analysis particularly the advent of
applying topological tools to networks of brain neurons has spawned
interest in constructing topological spaces out of digraphs, developing
computational tools for obtaining topological information, and using
these to understand networks. At the end of the day, (homological)
computations of the spaces reveal something about the geometric
realisation, thereby losing the directionality information.

However, digraphs can also be associated with path algebras. We can now
consider applying Hochschild homology to extract information, hopefully
obtaining something more refined in terms of the combinatorics of the
directed edges and paths in the digraph. Unfortunately, Hochschild
homology tends to vanish beyond degree 1. We can overcome this by
considering different higher paths of simplices, and thus introduce
Hochschild homology of digraphs in higher degrees. Moreover, this
procedure gives an implementable persistence pipeline for network
analysis. This is a joint work with Luigi Caputi.