Symdyn: An automated algebraic solution for high-order quantum systems

Abstract

Many important quantum systems are characterized by Hamiltonians expressible as a linear combination of time-independent generators of a finite Lie algebra, H ̂ ( t ) = ∑ l = 1 L η l ( t ) g ̂ l . The Wei-Norman framework provides a powerful analytical method for computing the exact time-evolution operator in a factorized form, U ̂ ( t ) = ∏ l = 1 L e Λ l ( t ) g ̂ l . By exploiting the underlying algebraic structure, the method yields a system of differential equations—the Wei-Norman equations—which determine the time-evolution operator coefficients Λ l ( t ) from the Hamiltonian ones η l ( t ) . Although robust and formally extendable, applying the method by hand becomes increasingly difficult or even intractable for high-order systems ( L ≥ 6 ). In this work, we introduce , a Python library that automates the application of this method. Given the commutators between generators, the library symbolically derives the Wei-Norman equations and efficiently computes nested commutators, similarity transformations, and the differential equations that yield Baker-Campbell-Hausdorff-like relations. We validate the library by recovering known results for quantum systems associated with the Lie algebras su ( 1 , 1 ) , su ( 2 ) , sl ( 2 ) , and so ( 2 , 1 ) , and further demonstrate its capabilities by deriving the Wei-Norman equations for a system of two coupled time-dependent harmonic oscillators involving 11 generators. Additionally, we specialize the library to handle the Lie group SU ( N ) using a generic Cartan-Weyl basis for its algebra, and demonstrate its versatility by deriving the Wei-Norman equations for su ( 2 ) , su ( 3 ) , and su ( 4 ) , including the explicit construction of a set of quantum gates for universal quantum computation.