Abstract:Quantum computing holds great potential for simulating quantum systems, such as molecular systems, due to its inherent ability to represent and manipulate quantum states. However, simulating nonunitary dynamics of open quantum systems on a quantum computer is still challenging. Here, we present a scheme denoted as VQS-QJ-LTLME, which adopts the trajectory average of quantum jump Monte Carlo wave function variational evolution to evolve the system's density matrix and local thermalizing Lindblad operators to describe the system-environment interactions. This combination allows the quantum circuit to be initialized after a period of evolution, thereby eliminating the accumulation of errors. The VQS-QJ-LTLME algorithm requires only log_2 qubits for system size $n$ and holds a time complexity of O(Tn^3), making it particularly suitable for operation on noisy intermediate-scale quantum devices. To accelerate the Monte Carlo sampling process in VQS-QJ-LTLME, we introduce an efficient sampling method, named No-evolution sampling (NES). The VQS-QJ-LTLME aided by NES requires simulating only (n^3+1) trajectories on quantum computers, and then makes 10^6 ~ 10^7 samplings on classical computers in a few seconds. To demonstrate the performance of the VQS-QJ-LTLME algorithm, we simulate the dynamics of a two-level spin-boson model and a four-level reduced Fenna-Matthews-Olson system with real superconducting quantum computers and classical simulators. It is shown that the simulations produced by VQS-QJ-LTLME algorithm align closely with those obtained from the purely classical methods.

What problem does this paper attempt to address?

### What problem does this paper attempt to solve? This paper aims to solve the challenge of simulating non - unitary dynamics of open quantum systems on quantum computers. Specifically, the paper proposes a new algorithm named VQS - QJ - LTLME (Variational Quantum Simulation with Quantum Jump and Local Thermalizing Lindblad Master Equation) to achieve the following goals: 1. **Simulate non - unitary evolution**: - The evolution of open quantum systems (such as molecular systems interacting with the environment) is non - unitary, while quantum circuits are essentially unitary. This leads to difficulties in directly simulating open quantum systems on quantum computers. - VQS - QJ - LTLME solves this problem by combining variational quantum simulation (VQS), quantum jump Monte Carlo wave function evolution and local thermalizing Lindblad operators. 2. **Reduce error accumulation**: - In quantum computing, as the evolution time increases, errors will gradually accumulate, affecting the simulation accuracy. - By introducing self - initializing local thermalizing Lindblad operators, VQS - QJ - LTLME can re - initialize the quantum circuit after each quantum jump, thus reducing error accumulation. 3. **Improve resource utilization efficiency**: - Traditional methods require a large number of qubits and deeper quantum circuits to simulate density matrices. - VQS - QJ - LTLME only requires $\log_2(n)$ qubits (where $n$ is the system size), and has a time complexity of $O(Tn^3)$, making it particularly suitable for running on noisy intermediate - scale quantum (NISQ) devices. 4. **Accelerate Monte Carlo sampling**: - The paper introduces a new sampling method - non - evolution sampling (NES), which can perform a large number of samplings quickly on a classical computer, significantly reducing the dependence on quantum resources. - NES only needs to simulate $n^3 + 1$ trajectories on a quantum computer and then perform $10^6\sim10^7$ samplings on a classical computer, greatly improving efficiency. 5. **Verify algorithm performance**: - The paper demonstrates the effectiveness of the VQS - QJ - LTLME algorithm by simulating two benchmark systems - the two - level spin - Bose model (SBM) and the four - level simplified Fenna - Matthews - Olson (FMO) system. - The results show that the VQS - QJ - LTLME algorithm can accurately describe the evolution of the mixed - state density matrix of the system and its decoherence process, and its results are highly consistent with those obtained by classical methods. ### Summary The main contribution of this paper is to propose an innovative quantum algorithm VQS - QJ - LTLME, which can efficiently and accurately simulate the non - unitary dynamics of open quantum systems on quantum computers, solving the limitations of existing methods in terms of resource requirements and error control.

Integrating Self-initialized Local Thermalizing Lindblad Operators for Variational Quantum Algorithm with Quantum Jump: Implementation and Performance