Convergence analysis of generalized modified splitting methods using multi-index series

We consider splitting methods for partial differential equations involving unbounded operators. For non-time-reversible dynamics, such as dissipative systems, negative splitting coefficients are generally not admissible because they require stepping backward in time, leading to an order barrier when all coefficients are required to be positive. We introduce generalized modified splitting methods to overcome this barrier. To analyze their convergence, we develop the corresponding multi-index series formalism, which provides a systematic framework for deriving order conditions. Using this formalism, we derive the order conditions and construct a generalized modified splitting method of order $6$. We also provide Python scripts that automate the generation and verification of order conditions, as well as the construction of new generalized modified splitting methods. Finally, we establish connections between the introduced multi-index series and related series formalisms from the literature, in particular, word series, Lie-Butcher series, and multi-index Butcher series.