Variational Quantum Algorithm for Solving the Liouvillian Gap

Abstract In open quantum systems, the Liouvillian gap characterizes the relaxation time toward the steady state. However, accurately computing this quantity is notoriously difficult due to the exponential growth of the Hilbert space and the non-Hermitian nature of the Liouvillian superoperator. In this work, we propose a variational quantum algorithm for efficiently estimating the Liouvillian gap. By utilizing the Choi–Jamiołkowski isomorphism, we reformulate the problem as finding the first excitation energy of an effective non-Hermitian Hamiltonian. Our method employs variance minimization with an orthogonality constraint to locate the first excited state and adopts a two-stage optimization scheme to enhance convergence. Moreover, to address scenarios with degenerate steady states, we introduce an iterative energy-offset scanning technique. Numerical simulations on the dissipative XXZ model confirm the accuracy and robustness of our algorithm across a range of system sizes and dissipation strengths. These results demonstrate the promise of variational quantum algorithms for simulating open quantum many-body systems on near-term quantum hardware.