Simulating Open Quantum System Dynamics on NISQ Computers with Generalized Quantum Master Equations

Yuchen Wang,Ellen Mulvihill,Zixuan Hu,Ningyi Lyu,Saurabh Shivpuje,Yudan Liu,Micheline B. Soley,Eitan Geva,Victor S. Batista,Sabre Kais

DOI: https://doi.org/10.1021/acs.jctc.3c00316

2023-08-08

Abstract:We present a quantum algorithm based on the Generalized Quantum Master Equation (GQME) approach to simulate open quantum system dynamics on noisy intermediate-scale quantum (NISQ) computers. This approach overcomes the limitations of the Lindblad equation, which assumes weak system-bath coupling and Markovity, by providing a rigorous derivation of the equations of motion for any subset of elements of the reduced density matrix. The memory kernel resulting from the effect of the remaining degrees of freedom is used as input to calculate the corresponding non-unitary propagator. We demonstrate how the Sz.-Nagy dilation theorem can be employed to transform the non-unitary propagator into a unitary one in a higher-dimensional Hilbert space, which can then be implemented on quantum circuits of NISQ computers. We validate our quantum algorithm as applied to the spin-boson benchmark model by analyzing the impact of the quantum circuit depth on the accuracy of the results when the subset is limited to the diagonal elements of the reduced density matrix. Our findings demonstrate that our approach yields reliable results on NISQ IBM computers.

Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to simulate the dynamics of open quantum systems on quantum computers. Specifically, the paper proposes a quantum algorithm based on the Generalized Quantum Master Equation (GQME) to overcome the limitations of the Lindblad equation. The Lindblad equation assumes that the coupling between the system and the environment is weak and satisfies the Markov approximation, while the GQME method can provide a strict derivation of the equations of motion for any sub - intensively reduced density matrix elements. Through this method, the dynamic behavior of open quantum systems on noisy intermediate - scale quantum (NISQ) computers can be simulated more accurately.

### Main Contributions
1. **Treatment of Non - Unitary Propagators**:
- The paper addresses a major challenge in quantum simulation, namely that the time - evolution operator is non - unitary while quantum gates are unitary. To this end, the authors use the Sz. - Nagy unitary dilation theorem to convert the non - unitary propagator into a unitary propagator in the extended Hilbert space, which can be implemented in a quantum circuit.

2. **Application of the Generalized Quantum Master Equation**:
- Through the GQME method, the paper provides a more general framework to describe the dynamics of open quantum systems without relying on the Markov and Born approximations. This enables the algorithm to be applied to a wider range of systems, including complex non - Markovian photosynthesis and photovoltaic systems.

3. **Verification of the Quantum Algorithm**:
- The authors verify the effectiveness of this quantum algorithm by simulating the spin - boson model on IBM quantum computers. The results show that the algorithm performs well in simulating the dynamics of the spin - boson model, especially on the QASM simulator.

### Technical Details
- **Generalized Quantum Master Equation (GQME)**:
\[
\frac{d}{dt}\hat{\sigma}(t)=-\frac{i}{\hbar}\langle L\rangle_0^n\hat{\sigma}(t)-\int_0^t d\tau K(\tau)\hat{\sigma}(t - \tau)
\]
where \(\langle L\rangle_0^n\) is the projected overall system Liouvillian operator and \(K(\tau)\) is the memory kernel super - operator.

- **Calculation of Non - Unitary Propagators**:
\[
\hat{\sigma}(t)=G(t)\hat{\sigma}(0)
\]
By numerically solving the GQME equation, the non - unitary propagator \(G(t)\) can be obtained.

- **Unitary Dilation**:
\[
U_{G'}(t)=\begin{pmatrix}
G'(t)&D_{G'}(t)\\
D_{G'}^\dagger(t)&-G'^\dagger(t)
\end{pmatrix}
\]
where \(D_{G'}(t)=\sqrt{I - G'^\dagger(t)G'(t)}\).

### Experimental Results
- **QASM Simulator**:
- On the QASM simulator, the results of the algorithm are highly consistent with the exact results of classical calculations, verifying the effectiveness of the algorithm.

- **Real Quantum Devices**:
- On IBM quantum computers, although affected by noise, the algorithm can still capture some trends of the exact results, but the quantitative consistency is poor. The authors suggest improving the accuracy by reducing the circuit depth.

### Conclusion
The paper proposes a quantum algorithm based on GQME that can more accurately simulate the dynamics of open quantum systems on NISQ computers. Through the unitary dilation technique, the problem of non - unitary propagators is solved and verified on the QASM simulator. Future work can further improve the performance of the algorithm on real quantum devices by reducing the circuit depth.