“…where W k (X, Y ) is a homogeneous Lie polynomial in X and Y of degree k [4]. The first few terms are There are several methods to compute W k [5,6,7,8]. In particular, a recursive algorithm has been proposed in [9] to express directly W k with the minimum number of independent commutators required at each degree k. Similar to the BCH formula, the Zassenhaus formula is useful in many different fields: q-analysis in quantum groups [10], quantum nonlinear optics [13], the Schrödinger equation in the semiclassical regime [12], and splitting methods in numerical analysis [11], etc.…”