Simulating quantum circuits using efficient tensor network contraction algorithms with subexponential upper bound

Thorsten B. Wahl,Sergii Strelchuk
DOI: https://doi.org/10.1103/PhysRevLett.131.180601
2023-08-14
Abstract:We derive a rigorous upper bound on the classical computation time of finite-ranged tensor network contractions in $d \geq 2$ dimensions. Consequently, we show that quantum circuits of single-qubit and finite-ranged two-qubit gates can be classically simulated in subexponential time in the number of gates. Moreover, we present and implement an algorithm guaranteed to meet our bound and which finds contraction orders with vastly lower computational times in practice. In many practically relevant cases this beats standard simulation schemes and, for certain quantum circuits, also a state-of-the-art method. Specifically, our algorithm leads to speedups of several orders of magnitude over naive contraction schemes for two-dimensional quantum circuits on as little as an $8 \times 8$ lattice. We obtain similarly efficient contraction schemes for Google's Sycamore-type quantum circuits, instantaneous quantum polynomial-time circuits, and non-homogeneous (2+1)-dimensional random quantum circuits.
Quantum Physics,Statistical Mechanics,Strongly Correlated Electrons,High Energy Physics - Theory
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an effective method to simulate the classical computing time of quantum circuits, especially for the tensor network contraction problems in two - dimensional and higher - dimensional cases. Specifically, the authors aim to: 1. **Derive strict upper bounds**: Derive a strict upper bound for the time complexity of classical computers when simulating tensor network contractions with a finite range. This includes quantum circuits composed of single - qubit and finite - range two - qubit gates. 2. **Propose efficient algorithms**: Design and implement an algorithm that can meet the above - mentioned upper bounds, and this algorithm can significantly reduce the computing time in practical applications. When dealing with two - dimensional quantum circuits, especially in small systems of 8x8 lattices, this algorithm is several orders of magnitude faster than traditional contraction schemes. 3. **Provide theoretical guarantees**: By using separation theorems (such as the Sphere - Separation Theorem and the Planar - Separation Theorem), the authors provide non - trivial theoretical running - time guarantees for tensor network contraction algorithms for the first time. These methods take into account the specific layout structure of quantum circuits, thereby improving efficiency. ### Specific content - Application of the **Sphere - Separation Theorem (SST)** and the **Planar - Separation Theorem (PST)**: These theorems are used to divide the tensor network into smaller parts, thereby simplifying the contraction process. For example, the SST shows that any set of balls in a d - dimensional space can be divided into two disjoint subsets by removing a part of the balls. \[ |ΓO(S)| \leq c_d k^{1/d} n^{1 - 1/d} \] \[ |ΓE(S)|, |ΓI(S)| \leq \frac{d + 1}{d + 2} n \] - **Classical simulation of quantum circuits**: By compressing the time dimension to zero and representing quantum gates as tensor networks, the authors prove that for any high - dimensional quantum circuit, the upper bound of its classical simulation time can be significantly better than that of the traditional layer - by - layer contraction method. For example, for a two - dimensional quantum circuit of nearest - neighbor two - qubit gates, the tighter upper bound obtained is: \[ t(0) < 2^{N^{1/\log_2(3/2)}} M a'_2 \sqrt{N} \] ### Experimental verification The paper also shows the practical performance of this method on several typical quantum circuits, including short - range instantaneous quantum polynomial - time (IQP) circuits, Google Sycamore - type circuits, and (2 + 1) - dimensional random quantum circuits. The experimental results show that the new method has surpassed the state - of - the - art methods in small - scale systems and also shows better performance in large - scale systems. ### Conclusion By introducing the Sphere - Separation Theorem and the Planar - Separation Theorem, the authors successfully provide strict upper bounds on the classical computing time for high - dimensional tensor network contractions and propose efficient contraction algorithms. These achievements not only help classical computers to benchmark quantum circuits but also help to determine the real boundaries of quantum supremacy.