The Magnus Expansion for Non-Hermitian Hamiltonians

The Magnus expansion offers a powerful geometric method to express a time-ordered exponential operator solution, typical of time-dependent linear differential equations, as an ordinary operatorial exponential. This representation possesses advantageous theoretical and qualitative features, such as preserving the underlying Lie algebraic structure of the system while remaining an exact solution to the original differential equation. For any finite-dimensional Hermitian Hamiltonian governing a quantum system, the standard Magnus expansion guarantees a manifestly unitary representation at every finite order of truncation. However, this foundational property is no longer preserved if the underlying Hamiltonian is non-Hermitian. In this work, we derive a generalized, unitarized expansion framework that systematically restores and maintains the property of geometric unitarity for all bounded finite-dimensional non-Hermitian Hamiltonians.