Control theory and splitting methods
Pith reviewed 2026-05-23 23:11 UTC · model grok-4.3
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.
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
- 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.
Referee Report
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)
- [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.
- [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)
- 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.
- [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
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
-
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
-
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
Minor self-citation on auxiliary basis; central claims independent
specific steps
-
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
axioms (2)
- standard math Lie algebra rank condition is equivalent to small-time local controllability for analytic vector fields (standard result in control theory)
- domain assumption Flows of the vector fields f0 and f1 exist and are sufficiently smooth
Reference graph
Works this paper leans on
-
[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
work page 2022
-
[2]
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
work page 2016
-
[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
work page 2006
-
[4]
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
work page 2022
-
[5]
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
work page 2023
-
[6]
K. Beauchard and F. Marbach. Quadratic obstructions to small- time local controllability for scalar-input systems. J. Differential Equations , 264(5):3704–3774, 2018
work page 2018
-
[7]
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]
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
work page 2005
- [9]
- [10]
- [11]
-
[12]
E. Bronasco and A. Laurent. Hopf algebra structures for th e backward error analysis of ergodic stochastic differential equations. Submitted, 2024
work page 2024
- [13]
-
[14]
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
work page 2022
-
[15]
J. C. Butcher. An algebraic theory of integration methods. Math. Comp. , 26:79–106, 1972
work page 1972
-
[16]
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
work page 2009
- [17]
- [18]
-
[19]
K.-T. Chen. Iterated integrals and exponential homomorphism s. Proc. London Math. Soc. (3), 4:502–512, 1954
work page 1954
-
[20]
K.-T. Chen. Integration of paths, geometric invariants and a g eneralized Baker-Hausdorff formula. Ann. Math. (2) , 65:163–178, 1957
work page 1957
-
[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
work page 1958
-
[22]
S. A. Chin. Symplectic integrators from composite operator fa ctorizations. Physics Letters A, 226(6):344–348, 1997
work page 1997
-
[23]
W.-L. Chow. ¨Uber Systeme von linearen partiellen Differentialgleichungen erster O rdnung. Math. Ann. , 117:98–105, 1939
work page 1939
-
[24]
J.-M. Coron. Control and nonlinearity , volume 136. Providence, RI: American Mathematical Society (AMS), 2007
work page 2007
- [25]
-
[26]
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
work page 2015
-
[27]
E. Faou. Geometric numerical integration and Schr¨ odinger equatio ns. Zurich Lectures in Advanced Mathematics. European Mathematical Society (EMS), Z ¨ urich, 2012
work page 2012
-
[28]
M. Fliess. Fonctionnelles causales non lin´ eaires et ind´ etermin´ ees non commutatives. Bull. Soc. Math. France, 109(1):3–40, 1981
work page 1981
-
[29]
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
work page 2016
-
[30]
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
work page 1996
- [31]
-
[32]
M. Hall. A basis for free Lie rings and higher commutators in free g roups. Proc. Am. Math. Soc., 1:575–581, 1950
work page 1950
-
[33]
E. Hansen and A. Ostermann. Exponential splitting for unboun ded operators. Mathematics of computation , 78(267):1485–1496, 2009
work page 2009
-
[34]
E. Hansen and A. Ostermann. High order splitting methods for a nalytic semigroups exist. BIT Numerical Mathematics , 49(3):527–542, 2009
work page 2009
-
[35]
M. Hochbruck and A. Ostermann. Explicit exponential Runge-K utta methods for semilinear parabolic problems. SIAM J. Numer. Anal. , 43(3):1069–1090, 2005
work page 2005
-
[36]
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
work page 2000
-
[37]
P.-V. Koseleff. Calcul formel pour les m´ ethodes de Lie en m´ ecanique hamilt onienne. PhD thesis, Ecole Polytechnique, 1993
work page 1993
-
[38]
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
work page 2020
-
[39]
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
work page 2013
-
[40]
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
work page 1997
-
[41]
V. T. Luan and A. Ostermann. Exponential B-series: the stiff c ase. SIAM J. Numer. Anal. , 51(6):3431–3445, 2013
work page 2013
-
[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
work page 2008
-
[43]
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
work page 2011
-
[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
work page 1954
-
[45]
R. I. McLachlan and G. R. W. Quispel. Splitting methods. Acta Numer. , 11:341–434, 2002
work page 2002
-
[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
work page 2013
-
[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
work page 2008
-
[48]
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
work page 2017
-
[49]
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
work page 2018
-
[50]
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
work page 2006
- [51]
-
[52]
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
work page 1993
- [53]
- [54]
-
[55]
Q. Sheng. Solving linear partial differential equations by expone ntial splitting. IMA Journal of numerical analysis , 9(2):199–212, 1989
work page 1989
- [56]
- [57]
- [58]
-
[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
work page 1991
-
[60]
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
work page 1984
-
[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
work page 1978
- [62]
-
[63]
E. Witt. Die Unterringe der freien Lieschen Ringe. Math. Z. , 64:195–216, 1956. 35
work page 1956
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.