Symplectic Reduction Research Articles

Conformal symplectic and constraint-preserving model order reduction of constrained conformal Hamiltonian systems

Dissipative and constrained dynamical systems model various physical phenomena comprising damping and constraining forces. These forces give rise to specific geometric properties of the dynamical system's governing equations. In this paper, we reduce and solve high-dimensional parameterized constrained and damped nonlinear differential equations with structure-preserving algorithms. These algorithms preserve the constrained equation's geometric properties, such as conformal symplecticness and phase-space, through model order reduction and numerical integration. We develop a reduced model based on an orthosymplectic reduced basis and establish the well-posedness and conformal symplecticness of the reduced model. Moreover, we derive a priori error bounds for state space and constraints errors. Subsequently, we utilize two RATTLE-based integrators to simulate the full and reduced models while preserving their geometric properties. Numerical experiments on constrained damped oscillators and a constrained lattice grid justify theoretical results and demonstrate the advantages of the structure-preserving computational methods for constrained dynamical systems.

Almost Commuting Scheme of Symplectic Matrices and Quantum Hamiltonian Reduction

Losev introduced the scheme X of almost commuting elements (i.e., elements commuting upto a rank one element) of g=sp(V)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}=\\mathfrak {sp}(V)$$\\end{document} for a symplectic vector space V and discussed its algebro-geometric properties. We construct a Lagrangian subscheme Xnil\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$X^{nil}$$\\end{document} of X and show that it is a complete intersection of dimension dim(g)+12dim(V)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ ext {dim}(\\mathfrak {g})+\\frac{1}{2}\ ext {dim}(V)$$\\end{document} and compute its irreducible onents. We also study the quantum Hamiltonian reduction of the algebra D(g)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {D}(\\mathfrak {g})$$\\end{document} of differential operators on the Lie algebra g\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathfrak {g}$$\\end{document} tensored with the Weyl algebra with respect to the action of the symplectic group, and show that it is isomorphic to the spherical subalgebra of a certain rational Cherednik algebra of Type C. We contruct a category Cc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {C}_c$$\\end{document} of D\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {D}$$\\end{document}-modules whose characteristic variety is contained in Xnil\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$X^{nil}$$\\end{document} and construct an exact functor from this category to the category O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}$$\\end{document} of the above rational Cherednik algebra. Simple objects of the category Cc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {C}_c$$\\end{document} are mirabolic analogs of Lusztig’s character sheaves.

Null Hamiltonian Yang–Mills theory: Soft Symmetries and Memory as Superselection

AbstractSoft symmetries for Yang–Mills theory are shown to correspond to the residual Hamiltonian action of the gauge group on the Ashtekar–Streubel phase space, which is the result of a partial symplectic reduction. The associated momentum map is the electromagnetic memory in the Abelian theory, or a nonlinear, gauge-equivariant, generalisation thereof in the non-Abelian case. This result follows from an application of Hamiltonian reduction by stages, enabled by the existence of a natural normal subgroup of the gauge group on a null codimension-1 submanifold with boundaries. The first stage is coisotropic reduction of the Gauss constraint, and it yields a symplectic extension of the Ashtekar–Streubel phase space (up to a covering). Hamiltonian reduction of the residual gauge action leads to the fully reduced phase space of the theory. This is a Poisson manifold, whose symplectic leaves, called superselection sectors, are labelled by the (gauge classes of the generalised) electric flux across the boundary. In this framework, the Ashtekar–Streubel phase space arises as an intermediate reduction stage that enforces the superselection of the electric flux at only one of the two boundary components. These results provide a natural, purely Hamiltonian, explanation of the existence of soft symmetries as a byproduct of partial symplectic reduction, as well as a motivation for the expected decomposition of the quantum Hilbert space of states into irreducible representations labelled by the Casimirs of the Poisson structure on the reduced phase space.

Monge–Ampère geometry and vortices

We introduce a new approach to Monge–Ampère geometry based on techniques from higher symplectic geometry. Our work is motivated by the application of Monge–Ampère geometry to the Poisson equation for the pressure that arises for incompressible Navier–Stokes flows. Whilst this equation constitutes an elliptic problem for the pressure, it can also be viewed as a non-linear partial differential equation connecting the pressure, the vorticity, and the rate-of-strain. As such, it is a key diagnostic relation in the quest to understand the formation of vortices in turbulent flows. We study this equation via an associated (higher) Lagrangian submanifold in the cotangent bundle to the configuration space of the fluid. Using our definition of a (higher) Monge–Ampère structure, we study an associated metric on the cotangent bundle together with its pull-back to the (higher) Lagrangian submanifold. The signatures of these metrics are dictated by the relationship between vorticity and rate-of-strain, and their scalar curvatures can be interpreted in a physical context in terms of the accumulation of vorticity, strain, and their gradients. We show explicity, in the case of two-dimensional flows, how topological information can be derived from the Monge–Ampère geometry of the Lagrangian submanifold. We also demonstrate how certain solutions to the three-dimensional incompressible Navier–Stokes equations, such as Hill’s spherical vortex and an integrable case of Arnol’d–Beltrami–Childress flow, have symmetries that facilitate a formulation of these solutions from the perspective of (higher) symplectic reduction.

Symplectic reduction of the sub-Riemannian geodesic flow for metabelian nilpotent groups

We consider nilpotent Lie groups for which the derived subgroup is abelian. We equip them with sub-Riemannian metrics and we study the normal Hamiltonian flow on the cotangent bundle. We show a correspondence between normal trajectories and polynomial Hamiltonians in some Euclidean space. We use the aforementioned correspondence to give a criterion for the integrability of the normal Hamiltonian flow. As an immediate consequence, we show that in Engel-type groups the flow of the normal Hamiltonian is integrable. For Carnot groups that are semidirect products of two abelian groups, we give a set of conditions that normal trajectories must fulfill to be globally length-minimizing. Our results are based on a symplectic reduction procedure.

Invariant Kähler potentials and symplectic reduction

For a proper Hamiltonian action of a Lie group $G$ on a Kahler manifold $(X,\omega)$ with momentum map $\mu$ we show that the symplectic reduction $\mu^{-1}(0)/G$ is a normal complex space. Every point in $\mu^{-1}(0)$ has a $G$-stable open neighborhood on which $\omega$ and $\mu$ are given by a $G$-invariant Kahler potential. This is used to show that $\mu^{-1}(0)/G$ is a Kahler space. Furthermore we examine the existence of potentials away from $\mu^{-1}(0)$ with both positive and negative results.

Algebraic and geometric reduction of multisymplectic manifolds

In this talk, we discuss an extension of the Marsden–Weinstein–Meyers symplectic reduction theorem to multisymplectic manifolds, and an adaptation of the Śniatycki–Weinstein, Dirac and Arms–Gotay–Cushman Poisson algebra reduction theorems to L∞ -algebras of multisymplectic observables. This is based on joint work with A Miti and L Ryvkin.

Optimization on the symplectic Stiefel manifold: SR decomposition-based retraction and applications

Numerous problems in optics, quantum physics, stability analysis, and control of dynamical systems can be brought to an optimization problem with matrix variable subjected to the symplecticity constraint. As this constraint nicely forms a so-called symplectic Stiefel manifold, Riemannian optimization is preferred, because one can borrow ideas from unconstrained optimization methods after preparing necessary geometric tools. Retraction is arguably the most important one which decides the way iterates are updated given a search direction. Two retractions have been constructed so far: one relies on the Cayley transform and the other is designed using quasi-geodesic curves. In this paper, we propose a new retraction which is based on an SR matrix decomposition. We prove that its domain contains the open unit ball which is essential in proving the global convergence of the associated gradient-based optimization algorithm. Moreover, we consider three applications—symplectic target matrix problem, symplectic eigenvalue computation, and symplectic model reduction of Hamiltonian systems—with various examples. The extensive numerical comparisons reveal the strengths of the proposed optimization algorithm.

The ABJM Amplituhedron

In this paper, we take a major step towards the construction and applications of an all-loop, all-multiplicity amplituhedron for three-dimensional planar N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 6 Chern-Simons matter theory, or the ABJM amplituhedron. We show that by simply changing the overall sign of the positive region of the original amplituhedron for four-dimensional planar N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 super-Yang-Mills (sYM) and performing a symplectic reduction, only three-dimensional kinematics in the middle sector of even-multiplicity survive. The resulting form of the geometry, combined with its parity images, gives the full loop integrand. This simple modification geometrically enforces the vanishing of odd-multiplicity cuts, and manifests the correct soft cuts as well as two-particle unitarity cuts. Furthermore, the so-called “bipartite structures” of four-point all-loop negative geometries also directly generalize to all multiplicities. We introduce a novel approach for triangulating loop amplituhedra based on the kinematics of the tree region, resulting in local integrands tailored to “prescriptive unitarity”. This construction sheds fascinating new light on the interplay between loop and tree amplituhedra for both ABJM and N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 sYM: the loop geometry demands that the tree region must be dissected into chambers, defined by the simultaneous positivity of maximal cuts. The loop geometry is then the “fibration” of the tree region. Using the new construction, we give explicit results of one-loop integrands up to ten points and two-loop integrands up to eight points by computing the canonical form of ABJM loop amplituhedron.

Quantum generalized Calogero–Moser systems from free Hamiltonian reduction

The one-dimensional system of particles with a 1/x2 repulsive potential is known as the Calogero–Moser system. Its classical version can be generalized by substituting the coupling constants with additional degrees of freedom, which span the so(N) or su(N) algebra with respect to Poisson brackets. We present the quantum version of this generalized model. As the classical generalization is obtained by a symplectic reduction of a free system, we present a method of obtaining a quantum system along similar lines. The reduction of a free quantum system results in a Hamiltonian, which preserves the differences in dynamics of the classical system depending on the underlying, orthogonal or unitary, symmetry group. The orthogonal system is known to be less repulsive than the unitary one, and the reduced free quantum Hamiltonian manifests this trait through an additional attractive term ∑i<j−ħ2(xi−xj)2, which is absent when one performs the straightforward canonical Dirac quantization of the considered system. We present a detailed and rigorous derivation of the generalized quantum Calogero–Moser Hamiltonian, we find the spectra and wavefunctions for the number of particles N=2,3, and we diagonalize the Hamiltonian partially for a general value of N.

Symplectic model reduction of Hamiltonian systems using data-driven quadratic manifolds

This work presents two novel approaches for the symplectic model reduction of high-dimensional Hamiltonian systems using data-driven quadratic manifolds. Classical symplectic model reduction approaches employ linear symplectic subspaces for representing the high-dimensional system states in a reduced-dimensional coordinate system. While these approximations respect the symplectic nature of Hamiltonian systems, linear basis approximations can suffer from slowly decaying Kolmogorov N-width, especially in wave-type problems, which then requires a large basis size. We propose two different model reduction methods based on recently developed quadratic manifolds, each presenting its own advantages and limitations. The addition of quadratic terms to the state approximation, which sits at the heart of the proposed methodologies, enables us to better represent intrinsic low-dimensionality in the problem at hand. Both approaches are effective for issuing predictions in settings well outside the range of their training data while providing more accurate solutions than the linear symplectic reduced-order models.

Adaptive symplectic model order reduction of parametric particle-based Vlasov–Poisson equation

High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on a set of parameters. In this work, we derive reduced order models for the semi-discrete Hamiltonian system resulting from a geometric particle-in-cell approximation of the parametric Vlasov–Poisson equations. Since the problem’s nondissipative and highly nonlinear nature makes it reducible only locally in time, we adopt a nonlinear reduced basis approach where the reduced phase space evolves in time. This strategy allows a significant reduction in the number of simulated particles, but the evaluation of the nonlinear operators associated with the Vlasov–Poisson coupling remains computationally expensive. We propose a novel reduction of the nonlinear terms that combines adaptive parameter sampling and hyper-reduction techniques to address this. The proposed approach allows decoupling the operations having a cost dependent on the number of particles from those that depend on the instances of the required parameters. In particular, in each time step, the electric potential is approximated via dynamic mode decomposition (DMD) and the particle-to-grid map via a discrete empirical interpolation method (DEIM). These approximations are constructed from data obtained from a past temporal window at a few selected values of the parameters to guarantee a computationally efficient adaptation. The resulting DMD-DEIM reduced dynamical system retains the Hamiltonian structure of the full model, provides good approximations of the solution, and can be solved at a reduced computational cost.

Reductions: precontact versus presymplectic

We show that contact reductions can be described in terms of symplectic reductions in the traditional Marsden–Weinstein–Meyer as well as the constant rank picture. The point is that we view contact structures as particular (homogeneous) symplectic structures. A group action by contactomorphisms is lifted to a Hamiltonian action on the corresponding symplectic manifold, called the symplectic cover of the contact manifold. In contrast to the majority of the literature in the subject, our approach includes general contact structures (not only co-oriented) and changes the traditional view point: contact Hamiltonians and contact moment maps for contactomorphism groups are no longer defined on the contact manifold itself, but on its symplectic cover. Actually, the developed framework for reductions is slightly more general than purely contact, and includes a precontact and presymplectic setting which is based on the observation that there is a one-to-one correspondence between isomorphism classes of precontact manifolds and certain homogeneous presymplectic manifolds.

A QP perspective on topology change in Poisson–Lie T-duality

We describe topological T-duality and Poisson–Lie T-duality in terms of QP (differential graded symplectic) manifolds and their canonical transformations. Duality is mediated by a QP-manifold on doubled non-abelian ‘correspondence’ space, from which we can perform mutually dual symplectic reductions, where certain canonical transformations play a vital role. In the presence of spectator coordinates, we show how the introduction of bibundle structure on correspondence space realises changes in the global fibration structure under Poisson–Lie duality. Our approach can be directly translated to the worldsheet to derive dual string current algebras. Finally, the canonical transformations appearing in our reduction procedure naturally suggest a Fourier–Mukai integral transformation for Poisson–Lie T-duality.

Time quasi-periodic vortex patches of Euler equation in the plane

We prove the existence of time quasi-periodic vortex patch solutions of the 2 $d$ -Euler equations in $\mathbb{R}^{2} $ , close to uniformly rotating Kirchhoff elliptical vortices, with aspect ratios belonging to a set of asymptotically full Lebesgue measure. The problem is reformulated into a quasi-linear Hamiltonian equation for a radial displacement from the ellipse. A major difficulty of the KAM proof is the presence of a zero normal mode frequency, which is due to the conservation of the angular momentum. The key novelty to overcome this degeneracy is to perform a perturbative symplectic reduction of the angular momentum, introducing it as a symplectic variable in the spirit of the Darboux-Carathéodory theorem of symplectic rectification, valid in finite dimension. This approach is particularly delicate in a infinite dimensional phase space: our symplectic change of variables is a nonlinear modification of the transport flow generated by the angular momentum itself. This is the first time such an idea is implemented in KAM for PDEs. Other difficulties are the lack of rotational symmetry of the equation and the presence of hyperbolic/elliptic normal modes. The latter difficulties -as well as the degeneracy of a normal frequency- are absent in other vortex patches problems which have been recently studied using the formulation introduced in this paper.

Symplectic Model Reduction of Hamiltonian Systems on Nonlinear Manifolds and Approximation with Weakly Symplectic Autoencoder

Classical model reduction techniques project the governing equations onto linear subspaces of the high-dimensional state-space. For problems with slowly decaying Kolmogorov- -widths such as certain transport-dominated problems, however, classical linear-subspace reduced-order models (ROMs) of low dimension might yield inaccurate results. Thus, the concept of classical linear-subspace ROMs has to be extended to more general concepts, like model order teduction (MOR) on manifolds. Moreover, as we are dealing with Hamiltonian systems, it is crucial that the underlying symplectic structure is preserved in the reduced model, as otherwise it could become unphysical in the sense that the energy is not conserved or stability properties are lost. To the best of our knowledge, existing literature addresses either MOR on manifolds or symplectic model reduction for Hamiltonian systems, but not their combination. In this work, we bridge the two aforementioned approaches by providing a novel projection technique called symplectic manifold Galerkin (SMG), which projects the Hamiltonian system onto a nonlinear symplectic trial manifold such that the reduced model is again a Hamiltonian system. We derive analytical results such as stability, energy-preservation, and a rigorous a posteriori error bound. Moreover, we construct a weakly symplectic deep convolutional autoencoder as a computationally practical approach to approximate a nonlinear symplectic trial manifold. Finally, we numerically demonstrate the ability of the method to achieve higher accuracy than (non-)structure-preserving linear-subspace ROMs and non-structure-preserving MOR on manifold techniques.

Algebraic Symplectic Reduction and Quantization of Singular Spaces

The algebraic method of singular reduction is applied for non regular group action on manifolds which provides singular symplectic spaces. The problem of deformation quantization of the singular surfaces is the focus. For some examples of singular Poisson spaces the deformation quantization is explicitly constructed. In is shown that for the flat phase space with the classical moment map and the orthogonal group action the deformation quantization converges for the entire arguments of exponential type.

Algebraic approach to the frozen formalism problem of time

The long-standing problem of time in canonical quantum gravity is the source of several conceptual and technical issues. Here, recent mathematical results are used to provide a consistent algebraic formulation of dynamical symplectic reduction that avoids difficult requirements such as the computation of a complete set of Dirac observables or the construction of a physical Hilbert space. In addition, the new algebraic treatment makes it possible to implement a consistent realization of the gauge structure off the constraint surface. As a consequence, previously unrecognized consistency conditions are imposed on deparametrization---the method traditionally used to unfreeze evolution in completely constrained systems. A detailed discussion of how the new formulation extends previous semiclassical results shows that an internal time degree of freedom need not be semiclassical in order to define a consistent quantum evolution.

Symplectic reduction and a Darboux–Moser–Weinstein theorem for Lie algebroids

Symplectic reduction and a Darboux–Moser–Weinstein theorem for Lie algebroids

Partial Dirac cohomology and tempered representations

The tempered representations of a real reductive Lie group G are naturally partitioned into series associated with conjugacy classes of Cartan subgroups H of G. We define partial Dirac cohomology, apply it for geometric construction of various models of these H–series representations, and show how this construction fits into the framework of geometric quantization and symplectic reduction.