Extend exponential-integrator error bounds to splitting discretizations on manifolds

Establish discretization error bounds for two-stage splitting schemes for Langevin-type dynamics on manifolds (including kinetic Langevin on Lie groups) that are analogous to the existing exponential-integrator or geometric Euler–Maruyama bounds, thereby enabling rigorous nonasymptotic analyses of such structure-preserving discretizations.

Background

The paper uses a two-stage operator splitting discretization to preserve Lie group structure, contrasting with more commonly analyzed exponential-integrator or geometric Euler–Maruyama schemes. While prior work has quantified discretization errors for manifold Langevin (without momentum), extending those bounds to splitting schemes remains technically challenging.

The authors note that existing tools tailored to exponential-integrator-based analysis have difficulties applying to their splitting method, highlighting a gap in discretization error theory for structure-preserving schemes on manifolds.