pith. sign in

arxiv: 2407.02127 · v2 · submitted 2024-07-02 · 🧮 math.NA · cs.NA· math.OC

Control theory and splitting methods

Pith reviewed 2026-05-23 23:11 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.OC
keywords splitting methodscontrol theoryLie algebra rank conditionsmall-time local controllabilitynumerical schemesforward flowsDirac controlsfree Lie algebra
0
0 comments X

The pith

Splitting methods achieve arbitrary order with only forward flows using complex coefficients for the reversible part.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper interprets numerical splitting methods for evolution equations with an irreversible vector field f0 as trajectories of a control-affine system driven by impulsive controls. This link to control theory yields the result that schemes of any order using solely forward flows of f0 exist when complex coefficients are allowed for f1. Equivalently, the Lie algebra rank condition becomes necessary and sufficient for small-time local controllability under complex-valued controls. For real coefficients the known order barriers are traced to bad Lie brackets that block controllability, and a free Lie algebra basis is used to identify when higher orders remain possible.

Core claim

There exist numerical schemes of arbitrary order involving only forward flows of f0, provided one allows complex coefficients for f1. Equivalently, for complex-valued controls, the Lie algebra rank condition is equivalent to small-time local controllability. For real-valued coefficients, order restrictions are linked to bad Lie brackets, and high-order methods are investigated via a basis of the free Lie algebra.

What carries the argument

Interpretation of splitting methods as control-affine trajectories with Dirac-mass controls, combined with the Lie algebra rank condition and the authors' free Lie algebra basis.

If this is right

  • Complex controls remove all order barriers for forward-only splitting schemes whenever the Lie algebra rank condition holds.
  • Real controls remain subject to order limits precisely when bad Lie brackets are present.
  • The free Lie algebra basis supplies explicit bracket combinations that allow construction of higher-order real schemes when they exist.
  • Classical splitting results can be recovered and extended by applying standard controllability tools to the impulsive-control formulation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The complex-coefficient construction could be tested numerically on low-dimensional non-reversible ODEs to confirm achievable orders.
  • The same control interpretation might suggest new splitting approaches for irreversible PDEs where only forward flows are stable.
  • Connections between the bad-bracket obstructions and existing barrier results in other numerical contexts remain open for exploration.

Load-bearing premise

The vector fields f0 and f1 generate a Lie algebra to which the rank condition and the free-Lie-algebra basis can be applied.

What would settle it

An explicit low-dimensional system satisfying the Lie algebra rank condition but admitting no arbitrary-order splitting scheme with complex coefficients would falsify the claimed equivalence.

read the original abstract

Our goal is to highlight some deep connections between numerical splitting methods and control theory. We consider evolution equations of the form $\dot{x} = f_0(x) + f_1(x)$, where $f_0$ encodes non-reversible dynamics, motivating schemes that involve only forward flows of $f_0$. In this context, a splitting method can be interpreted as a trajectory of the control-affine system $\dot{x}(t)=f_0(x(t))+u(t)f_1(x(t))$, associated with a control $u$ that is a finite sum of Dirac masses. The goal is then to find a control such that the flow generated by $f_0 + u(t)f_1$ is as close as possible to the flow of $f_0+f_1$. Using this interpretation and classical tools from control theory, we revisit well-known results on numerical splitting methods and prove several new ones. First, we show that there exist numerical schemes of arbitrary order involving only forward flows of $f_0$, provided one allows complex coefficients for $f_1$. Equivalently, for complex-valued controls, we prove that the Lie algebra rank condition is equivalent to small-time local controllability. Second, for real-valued coefficients, we show that the well-known order restrictions are linked to so-called "bad" Lie brackets from control theory, which are known to obstruct small-time local controllability. We investigate the conditions under which high-order methods exist, thanks to a basis of the free Lie algebra that we recently constructed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript interprets splitting methods for evolution equations of the form ẋ = f₀(x) + f₁(x) as trajectories of the control-affine system ẋ = f₀(x) + u(t)f₁(x) driven by Dirac-mass controls. It claims that, for complex-valued controls, arbitrary-order schemes exist that use only forward flows of f₀, and that the Lie algebra rank condition is equivalent to small-time local controllability. For real-valued controls, order barriers are linked to “bad” Lie brackets, and a recently constructed free-Lie-algebra basis is used to investigate when high-order real methods are possible.

Significance. If the central equivalences and constructions hold, the work supplies a control-theoretic explanation for known order restrictions in real splitting methods and opens a route to arbitrary-order schemes via complex coefficients. The explicit use of the authors’ prior free-Lie-algebra basis for the real case is a concrete technical contribution that could be reusable in other contexts.

major comments (2)
  1. [Abstract] Abstract (and § on complex controls): the asserted equivalence between LARC and STLC for complex-valued controls is presented as a new result obtained via classical tools, yet the manuscript does not indicate whether the proof requires any modification of the standard real-analytic controllability arguments or simply invokes them verbatim; an explicit comparison with the real case would strengthen the claim.
  2. [Real-coefficient case] Section on real coefficients and the free-Lie-algebra basis: the link between “bad” brackets and order restrictions is asserted, but the manuscript does not supply a concrete example (e.g., a specific vector-field pair and bracket that blocks order 3 or higher) showing how the basis construction directly yields the obstruction; without such an illustration the practical utility of the basis for splitting remains unclear.
minor comments (2)
  1. Notation: the distinction between the uncontrolled vector field f₀ and the controlled term u f₁ should be made uniform throughout; occasional reuse of f₀ + f₁ for both the target equation and the controlled system is confusing.
  2. [Abstract] The abstract states that the schemes involve “only forward flows of f₀”; a short remark clarifying whether this also excludes backward flows of f₁ (or treats them as complex) would help readers unfamiliar with the control interpretation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive suggestions. We address each major comment below and will make the indicated revisions.

read point-by-point responses
  1. Referee: [Abstract] Abstract (and § on complex controls): the asserted equivalence between LARC and STLC for complex-valued controls is presented as a new result obtained via classical tools, yet the manuscript does not indicate whether the proof requires any modification of the standard real-analytic controllability arguments or simply invokes them verbatim; an explicit comparison with the real case would strengthen the claim.

    Authors: The proof invokes the classical real-analytic controllability theorems verbatim; the complex case is reduced to an equivalent real system of doubled dimension to which the standard arguments apply without modification. We will insert a clarifying sentence in the abstract and the relevant section, together with a brief comparison noting that the real-valued case requires extra conditions precisely because of the obstructing 'bad' brackets. revision: yes

  2. Referee: [Real-coefficient case] Section on real coefficients and the free-Lie-algebra basis: the link between “bad” brackets and order restrictions is asserted, but the manuscript does not supply a concrete example (e.g., a specific vector-field pair and bracket that blocks order 3 or higher) showing how the basis construction directly yields the obstruction; without such an illustration the practical utility of the basis for splitting remains unclear.

    Authors: We agree that a concrete illustration would make the utility of the basis clearer. We will add a short example consisting of two explicit vector fields in which a bad bracket of weight 3 is identified by the free-Lie-algebra basis and directly prevents the existence of a real splitting method of order greater than 2. revision: yes

Circularity Check

1 steps flagged

Minor self-citation on auxiliary basis; central claims independent

specific steps
  1. self citation load bearing [Abstract, final sentence]
    "We investigate the conditions under which high-order methods exist, thanks to a basis of the free Lie algebra that we recently constructed."

    The sentence invokes the authors' own prior construction as the tool enabling the investigation of real-coefficient high-order schemes. While this is a minor auxiliary reference and does not affect the complex-coefficient equivalence result, it constitutes a self-citation whose load-bearing role for the real case is acknowledged in the text.

full rationale

The paper's core results establish arbitrary-order splitting schemes with complex coefficients via the equivalence of LARC and STLC for complex controls, relying on classical Lie-algebraic controllability tools. The sole self-citation appears in the secondary investigation of real-coefficient high-order methods, where the authors' prior free-Lie-algebra basis is invoked to link order restrictions to bad brackets. This citation is not load-bearing for the primary claims, introduces no self-referential definitions or fitted predictions, and does not reduce any stated equivalence or existence result to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper relies on standard results from nonlinear control theory and differential geometry; no free parameters or new postulated entities are introduced in the abstract.

axioms (2)
  • standard math Lie algebra rank condition is equivalent to small-time local controllability for analytic vector fields (standard result in control theory)
    Invoked when the abstract equates the rank condition with controllability for complex controls.
  • domain assumption Flows of the vector fields f0 and f1 exist and are sufficiently smooth
    Required for the splitting interpretation and for the control-affine system to be well-defined.

pith-pipeline@v0.9.0 · 5820 in / 1283 out tokens · 64710 ms · 2026-05-23T23:11:16.222332+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages

  1. [1]

    M. J. H. Al-Kaabi, K. Ebrahimi-Fard, D. Manchon, et al. Post-Lie M agnus expansion and BCH-recursion. SIGMA. Symmetry, Integrability and Geometry: Methods and A pplications, 18:023, 2022. 31

  2. [2]

    Alamo and J

    A. Alamo and J. M. Sanz-Serna. A technique for studying strong and weak local errors of splitting stochastic integrators. SIAM J. Numer. Anal. , 54(6):3239–3257, 2016

  3. [3]

    A. D. Bandrauk, E. Dehghanian, and H. Lu. Complex integration s teps in decomposition of quantum exponential evolution operators. Chemical physics letters , 419(4-6):346–350, 2006

  4. [4]

    Beauchard, J

    K. Beauchard, J. Le Borgne, and F. Marbach. Growth of struc ture constants of free Lie algebras relative to Hall bases. Journal of Algebra , 612:281–378, 2022

  5. [5]

    Beauchard, J

    K. Beauchard, J. Le Borgne, and F. Marbach. On expansions fo r nonlinear systems error estimates and convergence issues. Comptes Rendus. Math´ ematique, 361(G1):97–189, 2023

  6. [6]

    Beauchard and F

    K. Beauchard and F. Marbach. Quadratic obstructions to small- time local controllability for scalar-input systems. J. Differential Equations , 264(5):3704–3774, 2018

  7. [7]

    Beauchard and F

    K. Beauchard and F. Marbach. A unified approach of obstructio ns to small-time local con- trollability for scalar-input systems. arXiv preprint arXiv:2205.14114 , 2022

  8. [8]

    Blanes and F

    S. Blanes and F. Casas. On the necessity of negative coefficients for operator splitting schemes of order higher than two. Appl. Numer. Math. , 54(1):23–37, 2005

  9. [9]

    Blanes, F

    S. Blanes, F. Casas, P. Chartier, and A. Murua. Optimized high-o rder splitting methods for some classes of parabolic equations. Mathematics of Computation , 82(283):1559–1576, 2013

  10. [10]

    Blanes, F

    S. Blanes, F. Casas, and A. Escorihuela-Tom` as. Applying splitt ing methods with complex coefficients to the numerical integration of unitary problems. J. Comput. Dyn. , 9(2):85–101, 2022

  11. [11]

    Blanes, F

    S. Blanes, F. Casas, and A. Murua. Splitting methods for differe ntial equations. arXiv preprint arXiv:2401.01722, 2024

  12. [12]

    Bronasco and A

    E. Bronasco and A. Laurent. Hopf algebra structures for th e backward error analysis of ergodic stochastic differential equations. Submitted, 2024

  13. [13]

    Y. A. Bronsard, Y. Bruned, G. Maierhofer, and K. Schratz. S ymmetric resonance based integrators and forest formulae. arXiv preprint arXiv:2305.16737 , 2023

  14. [14]

    Bruned and K

    Y. Bruned and K. Schratz. Resonance-based schemes for dis persive equations via decorated trees. In Forum of Mathematics, Pi , volume 10, page e2. Cambridge University Press, 2022

  15. [15]

    J. C. Butcher. An algebraic theory of integration methods. Math. Comp. , 26:79–106, 1972

  16. [16]

    Castella, P

    F. Castella, P. Chartier, S. Descombes, and G. Vilmart. Splitting methods with complex times for parabolic equations. BIT Numer. Math. , 49(3):487–508, 2009

  17. [17]

    N. K. Chada, B. Leimkuhler, D. Paulin, and P. A. Whalley. Unbiased kinetic Langevin Monte Carlo with inexact gradients. arXiv:2311.05025, 2024

  18. [18]

    Chambers

    J. Chambers. Symplectic integrators with complex time steps. The Astronomical Journal , 126(2):1119, 2003

  19. [19]

    K.-T. Chen. Iterated integrals and exponential homomorphism s. Proc. London Math. Soc. (3), 4:502–512, 1954

  20. [20]

    K.-T. Chen. Integration of paths, geometric invariants and a g eneralized Baker-Hausdorff formula. Ann. Math. (2) , 65:163–178, 1957

  21. [21]

    K.-T. Chen, R. Fox, and R. Lyndon. Free differential calculus. I V: The quotient groups of the lower central series. Ann. Math. (2) , 68:81–95, 1958. 32

  22. [22]

    S. A. Chin. Symplectic integrators from composite operator fa ctorizations. Physics Letters A, 226(6):344–348, 1997

  23. [23]

    W.-L. Chow. ¨Uber Systeme von linearen partiellen Differentialgleichungen erster O rdnung. Math. Ann. , 117:98–105, 1939

  24. [24]

    J.-M. Coron. Control and nonlinearity , volume 136. Providence, RI: American Mathematical Society (AMS), 2007

  25. [25]

    Curry, K

    C. Curry, K. Ebrahimi-Fard, and B. Owren. The Magnus expans ion and post-Lie algebras. Mathematics of Computation , 89(326):2785–2799, 2020

  26. [26]

    Ebrahimi-Fard, A

    K. Ebrahimi-Fard, A. Lundervold, and H. Z. Munthe-Kaas. On t he Lie enveloping algebra of a post-Lie algebra. J. Lie Theory , 25(4):1139–1165, 2015

  27. [27]

    E. Faou. Geometric numerical integration and Schr¨ odinger equatio ns. Zurich Lectures in Advanced Mathematics. European Mathematical Society (EMS), Z ¨ urich, 2012

  28. [28]

    M. Fliess. Fonctionnelles causales non lin´ eaires et ind´ etermin´ ees non commutatives. Bull. Soc. Math. France, 109(1):3–40, 1981

  29. [29]

    Fløystad and H

    G. Fløystad and H. Munthe-Kaas. Pre- and post-Lie algebras: the algebro-geometric view. In Computation and Combinatorics in Dynamics, Stochastics an d Control: The Abel Symposium, Rosendal, Norway, August 2016 , pages 321–367. Springer, 2018

  30. [30]

    Goldman and T

    D. Goldman and T. J. Kaper. N th-order operator splitting sche mes and nonreversible systems. SIAM journal on numerical analysis , 33(1):349–367, 1996

  31. [31]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner. Geometric numerical integration , volume 31 of Springer Series in Computational Mathematics . Springer-Verlag, Berlin, second edition, 2006. Structure-preserving algorithms for ordinary differential equat ions

  32. [32]

    M. Hall. A basis for free Lie rings and higher commutators in free g roups. Proc. Am. Math. Soc., 1:575–581, 1950

  33. [33]

    Hansen and A

    E. Hansen and A. Ostermann. Exponential splitting for unboun ded operators. Mathematics of computation , 78(267):1485–1496, 2009

  34. [34]

    Hansen and A

    E. Hansen and A. Ostermann. High order splitting methods for a nalytic semigroups exist. BIT Numerical Mathematics , 49(3):527–542, 2009

  35. [35]

    Hochbruck and A

    M. Hochbruck and A. Ostermann. Explicit exponential Runge-K utta methods for semilinear parabolic problems. SIAM J. Numer. Anal. , 43(3):1069–1090, 2005

  36. [36]

    Iserles, H

    A. Iserles, H. Z. Munthe-Kaas, S. P. Nørsett, and A. Zanna. Lie-group methods. In Acta numerica, 2000, volume 9 of Acta Numer., pages 215–365. Cambridge Univ. Press, Cambridge, 2000

  37. [37]

    P.-V. Koseleff. Calcul formel pour les m´ ethodes de Lie en m´ ecanique hamilt onienne. PhD thesis, Ecole Polytechnique, 1993

  38. [38]

    Laurent and G

    A. Laurent and G. Vilmart. Exotic aromatic B-series for the stu dy of long time integrators for a class of ergodic SDEs. Math. Comp. , 89(321):169–202, 2020

  39. [39]

    Leimkuhler and C

    B. Leimkuhler and C. Matthews. Rational construction of stoc hastic numerical methods for molecular sampling. Appl. Math. Res. Express. AMRX , 2013(1):34–56, 2013

  40. [40]

    L´ opez-Marcos, J

    M. L´ opez-Marcos, J. M. Sanz-Serna, and R. D. Skeel. Explicit symplectic integrators using Hessian–vector products. SIAM Journal on Scientific Computing , 18(1):223–238, 1997. 33

  41. [41]

    V. T. Luan and A. Ostermann. Exponential B-series: the stiff c ase. SIAM J. Numer. Anal. , 51(6):3431–3445, 2013

  42. [42]

    C. Lubich. From quantum to classical molecular dynamics: reduced mode ls and numerical analysis. Zurich Lectures in Advanced Mathematics. European Mathematic al Society (EMS), Z¨ urich, 2008

  43. [43]

    Lundervold and H

    A. Lundervold and H. Munthe-Kaas. Hopf algebras of formal d iffeomorphisms and numerical integration on manifolds. In Combinatorics and physics , volume 539 of Contemp. Math., pages 295–324. Amer. Math. Soc., Providence, RI, 2011

  44. [44]

    W. Magnus. On the exponential solution of differential equation s for a linear operator. Com- munications on pure and applied mathematics , 7(4):649–673, 1954

  45. [45]

    R. I. McLachlan and G. R. W. Quispel. Splitting methods. Acta Numer. , 11:341–434, 2002

  46. [46]

    H. Z. Munthe-Kaas and A. Lundervold. On post-Lie algebras, L ie–Butcher series and moving frames. Foundations of Computational Mathematics , 13:583–613, 2013

  47. [47]

    H. Z. Munthe-Kaas and W. M. Wright. On the Hopf algebraic stru cture of Lie group integra- tors. Found. Comput. Math. , 8(2):227–257, 2008

  48. [48]

    Murua and J

    A. Murua and J. M. Sanz-Serna. Word series for dynamical sys tems and their numerical integrators. Found. Comput. Math. , 17(3):675–712, 2017

  49. [49]

    Murua and J

    A. Murua and J. M. Sanz-Serna. Hopf algebra techniques to ha ndle dynamical systems and numerical integrators. In Computation and combinatorics in dynamics, stochastics an d control, volume 13 of Abel Symp., pages 629–658. Springer, Cham, 2018

  50. [50]

    Prosen and I

    T. Prosen and I. Piˇ zorn. High order non-unitary split-step d ecomposition of unitary operators. Journal of Physics A: Mathematical and General , 39(20):5957, 2006

  51. [51]

    Rashevski

    P. Rashevski. About connecting two points of a completely nonh olonomic space by admissible curve. Uch. Zapiski Ped. Inst. Libknechta , 2:83–94, 1938

  52. [52]

    Reutenauer

    C. Reutenauer. Free Lie algebras , volume 7 of London Mathematical Society Monographs. New Series . The Clarendon Press, Oxford University Press, New York, 1993. Oxford Science Publications

  53. [53]

    R¨ oßler

    A. R¨ oßler. Rooted tree analysis for order conditions of stoch astic Runge-Kutta methods for the weak approximation of stochastic differential equations. Stoch. Anal. Appl. , 24(1):97–134, 2006

  54. [54]

    Rowlands

    G. Rowlands. A numerical algorithm for Hamiltonian systems. Journal of Computational Physics, 97(1):235–239, 1991

  55. [55]

    Q. Sheng. Solving linear partial differential equations by expone ntial splitting. IMA Journal of numerical analysis , 9(2):199–212, 1989

  56. [56]

    Shirshov

    A. Shirshov. On the bases of a free Lie algebra. Algebra Logika, 1(1):14–19, 1962

  57. [57]

    Sussmann

    H. Sussmann. Lie brackets and local controllability: a sufficient c ondition for scalar-input systems. SIAM J. Control Optim. , 21(5):686–713, 1983

  58. [58]

    Sussmann

    H. Sussmann. A general theorem on local controllability. SIAM J. Control Optim. , 25(1):158– 194, 1987

  59. [59]

    M. Suzuki. General theory of fractal path integrals with applic ations to many-body theories and statistical physics. Journal of Mathematical Physics , 32(2):400–407, 1991. 34

  60. [60]

    Takahashi and M

    M. Takahashi and M. Imada. Monte Carlo calculation of quantum systems. II. Higher order correction. Journal of the Physical Society of Japan , 53(11):3765–3769, 1984

  61. [61]

    G. Viennot. Alg` ebres de Lie libres et mono ¨ ıdes libres, volume 691 of Lecture Notes in Mathe- matics. Springer, Berlin, 1978. Bases des alg` ebres de Lie libres et factor isations des mono ¨ ıdes libres

  62. [62]

    Wisdom, M

    J. Wisdom, M. Holman, and J. Touma. Symplectic correctors. Fields Institute Communica- tions, 10:217, 1996

  63. [63]

    E. Witt. Die Unterringe der freien Lieschen Ringe. Math. Z. , 64:195–216, 1956. 35