Closed dynamical recursion equations for correlation functions and their application for the construction of the Liouvillian spectrum in Lindbladian systems

For an open quantum system described by the Lindblad equation, full characterization of its dynamics typically needs the knowledge of the Liouvillian spectrum and correlation functions. Solving the Liouvillian spectrum and correlation functions are usually formidable tasks, and most previous studies are constrained to simple models and lower-order correlations. In this work, we derive a closed form of dynamical recursion equations for all even-order correlation functions associated with the quadratic Liouvillian superoperator, thus extending Wick's theorem and enabling the reconstruction of the corresponding Liouvillian spectrum. Furthermore, we study the quartic Liouvillian superoperator, establishing the necessary and sufficient conditions for the closure of second-order correlation functions. Building on this, the dynamical expressions for even-order correlation functions under quadratic dissipation are derived, culminating in a method for constructing the Liouvillian spectrum in this extended framework.