Arxiv:2107.00677V2 [Quant-Ph] 19 Jul 2021

The fixed angle conjecture for QAOA on regular MaxCut graphs Jonathan Wurtz Department of Physics and Astronomy, Tufts University, Medford, Massachusetts 02155, USA∗ Danylo Lykov Argonne National Laboratory, Lemont, Illinois 60439, USA (Dated: July 20, 2021) The quantum approximate optimization algorithm (QAOA) is a near-term combinatorial optimization algorithm suitable for noisy quantum devices. However, little is known about performance guarantees for p > 2. A recent work [1] computing MaxCut performance guarantees for 3-regular graphs conjectures that any d-regular graph evaluated at particular fixed angles has an approximation ratio greater than some worst-case guarantee. In this work, we provide numerical evidence for this fixed angle conjecture for p < 12. We compute and provide these angles via numerical optimization and tensor networks. These fixed angles serve for an optimization-free version of QAOA, and have universally good performance on any 3 regular graph. Heuristic evidence is presented for the fixed angle conjecture on graph ensembles, which suggests that these fixed angles are \close" to global optimum. Under the fixed angle conjecture, QAOA has a larger performance guarantee than the Goemans Williamson algorithm on 3-regular graphs for p ≥ 11. I. INTRODUCTION ing a tensor network quantum circuit simulator [8]. In this work, we make this observation even more simple by Near-term quantum computers have the potential for demonstrating that there exist a set of fixed angles that advantage in the field of combinatorial optimization. Po- are \universally good". These fixed angles, when eval- tentially, an algorithm running on a small noisy quantum uated on any regular graph of fixed edge weights, will device may soon exhibit quantum advantage by provid- return approximation ratios larger than some guarantee ing better approximate to combinatorial problems than and very close to the exact maximal approximation ra- the best classical approximate solver [2]. One partic- tio. These fixed angles allow QAOA to completely bypass ular class of potential algorithms are variational quan- any variational optimization step, potentially increasing tum algorithms (VQA) [3], which use a hybrid quantum- execution times by a factor of 100 - 1000×. classical loop to optimize ansatz wavefunctions that en- code solutions to combinatorial problems. In this work, we focus on the NP-complete combina- One particular ansatz choice for VQA is the quantum torial optimization problem MaxCut [9]. Given a graph approximate optimization algorithm (QAOA) [4], which G of edges E and vertices V , a MaxCut algorithm strives is generated using p rounds of unitaries alternating be- to find a bipartition of vertices fA; V nAg such that a tween some \mixing" unitary and objective function uni- maximal number of edges E have one vertex in each bi- tary. The ansatz is a function of 2p parameters fγ; βg, partition. In other words, a MaxCut strives to \cut" the which are optimized through repeated query of a classical maximal number of edges via a bipartition. optimizer to a digital quantum device. A graph of N = jV j vertices and M = jEj edges can While it is known that the QAOA converges to the be encoded into an objective function C over N qubits exact result for p ! 1 consistent with the adiabatic and M clauses, as follows theorem [4;5], less is known about the guaranteed per- 1 X formance of the algorithm at low depth. In the original C^ = (1 − σî σ^j ): (1) introduction of Farhi et. al [4], the performance was guar- 2 z z hiji2E arXiv:2107.00677v2 [quant-ph] 19 Jul 2021 anteed to be at least 0:693 for p = 1 and 3 regular graphs, and in a later work [1], the performance was guaranteed The eigenstates of C^ are bipartitions, with eigenval- to be at least 0:7559 for p = 2 and 3 regular graphs. The ues counting the number of cut edges. The partition is same work observed that worst-case graphs have no small assigned by measuring the ansatz state in the Z basis. cycles, and conjectured that the same holds for larger p. The QAOA strives to find optimal parameters which In this work, we provide performance guarantees under maximize the expectation value of the objective function this conjecture for regular graphs for p ≤ 11, and show with respect to some ansatz state. The QAOA ansatz is that for p ≥ 11 the QAOA has a larger performance guar- defined as a function of 2p variables fγ; βg over p rounds antee on 3 regular graphs than the best general purpose of repeated unitaries classical solver. Additionally, it has been observed that optimal QAOA ^ ^ ^ ^ jγ; βi = e−iβpBe−iγpC (··· )e−iβ1Be−iγ1C j+i (2) angles concentrate around typical values [6;7], suggest- ing that QAOA may bypass the quantum-classical vari- where ellipsis represent the p actions of objective function ^ ^ P i ational optimization step, or pre-compute the angles us- C and mixing term B = i σ^x. The initial state j+i is 2 (A) (B) (C) (D) A graph A local subgraph A tensor treegraph decomposition Expectation values FIG. 1. A sketch of how to compute QAOA expectation values. Given some graph G of edges and vertices (A), QAOA will only \see" local structure within p steps of an edge. In this way, a graph can be decomposed into jEj subgraphs S (B), one for each edge, and expectation values be computed independently (D). One method of computing expectation values used in this work are tensor networks (C), which efficiently decompose the local QAOA state into a tensor network for efficient computation. hiji the equally weighted superposition state and maximal where fp is an individual edge contribution to the total eigenstate of B^. cost function. The QAOA ansatz only correlates vertices hiji The QAOA is a hybrid algorithm, in the sense that it within p steps, so each clause fp may be computed on includes both a quantum and a classical device. Through a subgraph Shiji that only includes the edges that are repeated query to a small quantum device to evaluate p G ^ incident from vertices at a distance p from the vertices expectation values Fp (γ; β) = hγ; βjCjγ; βi, a classical i or j. The expectation value thus can be calculated device optimizes the 2p angles fγ; βg to provide optimal as a sum over M independent subgraphs, as sketched in MaxCut solutions to a graph. Fig.1A-B. The performance of the algorithm is characterized by Due to this locality, QAOA cannot take advantage of the approximation ratio, which is the ratio between the G graph structures that are greater than 2p steps away [4]. optimized value Fp (γ; β) and the MaxCut value Additionally, QAOA cannot distinguish between graphs with cycles of size > 2p + 1, which suggests that worst- G case graphs have no small cycles. This intuition was G Fp (γ; β) Cp = MAX: G : (3) made concrete in [1], which proved that worst-case 3 γ,β Cmax regular graphs with the smallest approximation ratio for The approximation ratio ranges between 0 and 1. A p = 1 and 2 are bipartite with no cycles ≤ 3 or 5, re- larger indicates better performance, and a value of 1 in- spectively. A bipartite graph has a MaxCut value of M, dicates the exact solution. Given a class of graphs and cutting every edge. The only subgraphs of a graph with fixed value p, this approximation ratio may be guaran- no small cycles are the tree subgraph, which has no cycles teed to be above some performance guarantee. For 3 (see Fig.1). The optimized value of a worst-case graph G∗ is then regular graphs, C1 ≥ 0:692 [4], and C2 ≥ 0:7559 [1]. For p ! 1 the approximation ratio converges to 1 in accor- dance with the adiabatic theorem [4;5]. X F G∗ (γ; β) = f ij(γ; β) = Mf tree(γ; β): (5) Critically important to the evaluation of the perfor- p p p mance is the fact that QAOA is a local algorithm [10{ hiji2G∗ 13]. Given p steps, a vertex can only be correlated with Thus, the approximation ratio of a worst-case graph other vertices within a distance ≤ p. This is because the is simply the expectation value of the tree subgraph QAOA has a \lightcone" of interaction, which is clear to f tree(γ; β). The expectation value of the tree subgraph see in the Heisenberg picture (See for example, [1] Eq. p evaluated at optimal parameters fγ; βg serves as the 7). The expectation value of the cost function is a sum tree performance guarantee. These worst-case graphs are ex- of clauses ponentially rare but can be found constructively [14]. F G(γ; β) = hγ; βjC^jγ; βi Additionally, it was proven that any graph evaluated p at angles optimum to the tree subgraph have an approx- X 1 = hγ; βj (1 − σî σ^j ) jγ; βi imation ratio larger than the guarantee, for p = 1 and 2. 2 z z hiji2E These facts naturally lead to two conjectures for larger values of p [1] X hiji ≡ fp (γ; β); (4) hiji2E 3 Large Loop conjecture: The worst-case graphs for full state vector of the system and evolving it by apply- fixed p are bipartite and have no cycles less than 2p + 2. ing matrix transformations, the state is represented by a tensor network, and the gates with tensors that have in- Fixed Angle conjecture: Any graph evaluated at fixed put and output indices. The tensor network constructed angles optimal to the tree subgraph will have an approx- in this way can then be contracted in an efficient manner imation ratio larger than the guarantee.

Load more