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

arXiv:2511.19659v4 Announce Type: replace Abstract: Dissipation and irreversibility are central to many physical systems, 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 ill-conditioned for dissipative dynamics. We show how block encodings of complex-coefficient product formulas can be constructed by 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; the real parts preserve the contractive evolution, and the imaginary parts are additional unitary evolutions. We use known formulas of orders 4 and 6 and give a new symmetric formula of order 8 in this restricted class. We apply the approach to the classical problem of lossy mechanical wave propagation in statevector simulations and on a trapped-ion quantum processor. Orders 4 and 6 require comparable CNOT gates to orders 1 and 2 at large error tolerances, but quickly become more efficient as the tolerance is reduced. On the quantum hardware, a single step of order 4 achieves lower error than order 1 and comparable error to order 2, while order 6 is noise-limited. Our results show that operator splitting can combine practical circuit constructions with high-order convergence for non-unitary dissipative dynamics.

Abstract:Dissipation and irreversibility are central to many physical systems, 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 ill-conditioned for dissipative dynamics. We show how block encodings of complex-coefficient product formulas can be constructed by 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; the real parts preserve the contractive evolution, and the imaginary parts are additional unitary evolutions. We use known formulas of orders 4 and 6 and give a new symmetric formula of order 8 in this restricted class. We apply the approach to the classical problem of lossy mechanical wave propagation in statevector simulations and on a trapped-ion quantum processor. Orders 4 and 6 require comparable CNOT gates to orders 1 and 2 at large error tolerances, but quickly become more efficient as the tolerance is reduced. On the quantum hardware, a single step of order 4 achieves lower error than order 1 and comparable error to order 2, while order 6 is noise-limited. Our results show that operator splitting can combine practical circuit constructions with high-order convergence for non-unitary dissipative dynamics.

arXiv:2511.19659 [quant-ph]

(or arXiv:2511.19659v4 [quant-ph] for this version)

https://doi.org/10.48550/arXiv.2511.19659 Focus to learn more arXiv-issued DOI via DataCite

arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.

Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.

Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.

analiza aiqbench

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

Subiecte conexe