A construction of commuting systems of integrable symplectic birational maps

We give a construction of completely integrable $(2m)$-dimensional Hamiltonian systems with cubic Hamilton functions. The construction depends on a constant skew-Hamiltonian matrix $A$, that is, a matrix satisfying $A^{\rm T}J=JA$, where $J$ is a non-degenerate skew-symmetric matrix defining the standard symplectic structure on the phase space $\mathbb R^{2m}$. Applying to any such system the so called Kahan-Hirota-Kimura discretization scheme, we arrive at a birational $(2m)$-dimensional map. We show that this map is symplectic with respect to a symplectic structure that is a perturbation of the standard symplectic structure on $\mathbb R^{2m}$, and possesses $m$ independent integrals of motion, which are perturbations of the original Hamilton functions and are in involution with respect to the invariant symplectic structure. Thus, this map is completely integrable in the Liouville-Arnold sense. Moreover, under a suitable normalization of the original $m$-tuples of commuting vector fields, their Kahan-Hirota-Kimura discretizations also commute and share the invariant symplectic structure and the $m$ integrals of motion.