High-order splitting of non-unitary operators on quantum computers

TL;DR

This paper introduces a high-order operator splitting method using complex coefficients to simulate non-unitary, dissipative quantum dynamics accurately on near-term quantum hardware, overcoming stability issues of traditional methods.

Contribution

It develops a novel complex-coefficient product formula approach for high-order splitting of dissipative dynamics, enabling stable and accurate simulations on quantum computers.

Findings

High-order splitting improves simulation accuracy on noisy hardware.

Complex coefficients allow stable decomposition of dissipative processes.

Demonstrated effective simulation of lossy wave propagation on a trapped-ion processor.

Abstract

Dissipation and irreversibility are central to most physical processes, yet they lead to non-unitary dynamics that are challenging to realise on quantum processors. High-order operator splitting is an attractive approach for simulating unitary dynamics, yet conventional product formulas introduce negative time steps at high orders that are numerically unstable for dissipative dynamics. We show how complex-coefficient product formulas can decompose dissipative dynamics into a sequence of simple Hamiltonian evolutions in real and imaginary time with high-order accuracy. The unitary substages use positive real coefficients, while the dissipative substages use complex coefficients with positive real parts, where the real parts preserve the contractive evolution and the imaginary parts are additional unitary evolutions. We demonstrate the approach by simulating the classical problem of lossy mechanical wave propagation on a trapped-ion quantum processor. A step of order 4 achieves greater accuracy than the steps with low orders 1 and 2, despite the increased circuit depth on noisy hardware. The results suggest that high-order operator splitting is an accurate and practical approach for simulating dissipative dynamics on near-term quantum processors.