Abstract
We develop a data-driven framework for identifying non-Markovian dynamical equations of motion for open quantum systems. Starting from the Nakajima--Zwanzig formalism, we vectorize the reduced density matrix into a four-dimensional state vector and cast the dynamics as a Volterra integro-differential equation with an operator-valued memory kernel. The learning task is then formulated as a constrained optimization problem over the admissible operator space, where correlation functions are approximated by rational functions using Pad\'e approximants. We establish well-posedness of the learnin