Hybrid algorithm simulating non-equilibrium steady states of an open quantum system

Hongyi Zhou,Rui Mao,Xiaoming Sun

2023-09-13

Abstract:Non-equilibrium steady states are a focal point of research in the study of open quantum systems. Previous variational algorithms for searching these steady states have suffered from resource-intensive implementations due to vectorization or purification of the system density matrix, requiring large qubit resources and long-range coupling. In this work, we present a novel variational quantum algorithm that efficiently searches for non-equilibrium steady states by simulating the operator-sum form of the Lindblad equation. By introducing the technique of random measurement, we are able to estimate the nonlinear cost function while reducing the required qubit resources by half compared to previous methods. Additionally, we prove the existence of the parameter shift rule in our variational algorithm, enabling efficient updates of circuit parameters using gradient-based classical algorithms. To demonstrate the performance of our algorithm, we conduct simulations for dissipative quantum transverse Ising and Heisenberg models, achieving highly accurate results. Our approach offers a promising solution for effectively addressing non-equilibrium steady state problems while overcoming computational limitations and implementation challenges.

Quantum Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: efficiently searching for non - equilibrium steady states in open quantum systems. Specifically, the authors aim to overcome the challenges of existing variational algorithms in resource consumption and physical implementation, especially the large qubit resource requirements and long - range coupling problems caused by vectorizing or purifying the system density matrix.

### Problem Background

1. **Open Quantum Systems**: Quantum systems that interact with the environment, and their evolution cannot be described solely by the Schrödinger equation but requires the Lindblad master equation.
2. **Importance of Non - equilibrium Steady States**: These steady states are of great significance in statistical mechanics, biology, and quantum information processing tasks because they possess special transport properties.
3. **Limitations of Existing Methods**: Traditional methods such as matrix product states, tensor network schemes, and real - space renormalization methods, although effective, face computational resource limitations when simulating large - scale quantum systems. Quantum algorithms face difficulties in physical implementation, especially the need to integrate a large number of qubits and mitigate the decoherence effect.

### New Method Proposed in the Paper

The authors propose a new variational quantum algorithm to search for non - equilibrium steady states by directly simulating the operator - sum form of the Lindblad equation. The main innovations include:

- **Reducing Qubit Requirements**: Introducing random measurement techniques to halve the number of required qubits.
- **Efficient Parameter Updates**: Proving the existence of the parameter - shift rule, enabling efficient updating of circuit parameters using gradient - based classical algorithms.
- **High - Precision Simulation**: Simulating the dissipative transverse Ising model and the Heisenberg model, achieving high - precision results.

### Summary of Mathematical Formulas

The Lindblad master equation is described as:
\[ \frac{d\rho}{dt} = L(\rho)= - i[H, \rho]+\sum_{k = 1}^{K}\gamma_k\left(2c_k\rho c_k^{\dagger}-\{c_k c_k^{\dagger}, \rho\}\right), \]
where \( c_k \) is the \( k \) - th transition operator and \( \gamma_k \) is the dissipation intensity.

The steady - state condition is:
\[ L(\rho_{ss}) = 0. \]

The cost function is defined as the square of the Frobenius norm of the Lindblad equation:
\[ C(\vec{\theta}):=\|L(\rho)\|_F^2. \]

Through random measurement and classical shadow estimation, this quadratic function can be effectively estimated, and only half of the qubit resources are required.

### Conclusion

This paper provides an effective solution that can accurately find non - equilibrium steady states while reducing computational resources, thus providing new tools and methods for further research on open quantum systems.