Variational Approach to Entangled Non-Hermitian Open Systems

The pseudomode model effectively captures the nonperturbative dynamics of open quantum systems with significantly reduced degrees of freedom. However, it is limited by the exponential growth of the Hilbert space dimension. To overcome the computational challenges, we propose a novel method that combines the multiple Davydov Ansatz with the Choi-Jamiolkowski isomorphism. Within this framework, the Lindblad equation is transformed into the non-Hermitian Schrödinger equation in a double Hilbert space, with its dynamics determined using the time-dependent variational principle. Three cases are calculated to demonstrate the effectiveness of the proposed method. We first discuss how the Davydov Ansatz works for the model with a single pseudomode. Extending the method to multiple pseudomodes, we show that the Ansatz effectively circumvents the exponential growth of the Hilbert space. Additionally, the method is also capable of addressing potential intersections that emerge in multibath scenarios. This approach offers potential applicability to various types of pseudomode models and other dissipative systems, providing a promising tool for the studies of open quantum dynamics.