Qubit Routing using Graph Neural Network aided Monte Carlo Tree Search
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
PRELIMINARY PREPRINT VERSION: DO NOT CITE The AAAI Digital Library will contain the published version some time after the conference. Qubit Routing using Graph Neural Network aided Monte Carlo Tree Search Animesh Sinha, Utkarsh Azad, Harjinder Singh Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad. Center for Quantum Science and Technology, International Institute of Information Technology, Hyderabad. Abstract existing CNOT gates. This ensures that any non-local quan- Near-term quantum hardware can support two-qubit opera- tum operations are performed only between the qubits that tions only on the qubits that can interact with each other. are nearest-neighbors. This process of circuit transformation Therefore, to execute an arbitrary quantum circuit on the by a compiler routine for the target hardware is known as hardware, compilers have to first perform the task of qubit qubit routing (Cowtan et al. 2019). The output instructions in routing, i.e., to transform the quantum circuit either by insert- the transformed quantum circuit should follow the connec- ing additional SWAP gates or by reversing existing CNOT tivity constraints and essentially result in the same overall gates to satisfy the connectivity constraints of the target topol- unitary evolution as the original circuit (Pozzi et al. 2020). ogy. The depth of the transformed quantum circuits is mini- In the context of NISQ hardware, this procedure is of ex- mized by utilizing the Monte Carlo tree search (MCTS) to treme importance as the transformed circuit will, in general, perform qubit routing by making it both construct each action have higher depth due to the insertion of extra SWAP gates. and search over the space of all actions. It is aided in perform- ing these tasks by a Graph neural network that evaluates the This overhead in the circuit depth becomes more prominent value function and action probabilities for each state. Along due to the high decoherence rates of the qubits and it be- with this, we propose a new method of adding mutex-lock comes essential to find the most optimal and efficient strat- like variables in our state representation which helps factor in egy to minimize it (Cowtan et al. 2019; Herbert and Sen- the parallelization of the scheduled operations, thereby prun- gupta 2018; Pozzi et al. 2020). In this article, we present a ing the depth of the output circuit. Overall, our procedure (re- procedure that we refer to as QRoute. We use Monte Carlo ferred to as QRoute) performs qubit routing in a hardware ag- tree search (MCTS), which is a look-ahead search algorithm nostic manner, and it outperforms other available qubit rout- for finding optimal decisions in the decision space guided by ing implementations on various circuit benchmarks. a heuristic evaluation function (Kocsis and Szepesvari 2006; Munos 2014; Kocsis and Szepesvári 2006). We use it for ex- 1 Introduction plicitly searching the decision space for depth minimization The present-day quantum computers, more generally known and as a stable and performant machine learning setting. It is as Noisy Intermediate-Scale quantum (NISQ) devices aided by a Graph neural network (GNN) (Wang et al. 2018), (Preskill 2018) come in a variety of hardware architec- with an architecture that is used to learn and evaluate the tures (IBM 2021; Arute et al. 2019; Karalekas et al. 2020; heuristic function that will help guide the MCTS. Bourassa et al. 2021), but there exist a few problems per- vading across all of them. These problems constitute the 2 Qubit Routing poor quality of qubits, limited connectivity between qubits, In this section, we begin by defining the problem of qubit and the absence of error-correction for noise-induced errors routing formally and discussing the work done previously in encountered during the execution of gate operations. These the field. place a considerable restriction on the number of instruc- tions that can be executed to perform useful quantum com- 2.1 Describing the Problem putation (Preskill 2018). Collectively these instructions can The topology of quantum hardware can be visualized as a be realized as a sequential series of one or two-qubit gates qubit connectivity graph (Fig. 2). Each node in this graph that can be visualized more easily as a quantum circuit as would correspond to a physical qubit which in turn might shown in Fig. 1a (Childs, Schoute, and Unsal 2019). correspond to a logical qubit. The quantum instruction set, To execute an arbitrarily given quantum circuit on the tar- which is also referred to as quantum circuit (Fig. 1a), is a get quantum hardware, a compiler routine must transform it sequential series of single-qubit and two-qubit gate opera- to satisfy the connectivity constraints of the topology of the tions that act on the logical qubits. The two-qubit gates such hardware (Cowtan et al. 2019). These transformations usu- as CNOT can only be performed between two logical qubits ally include the addition of SWAP gates and the reversal of iff there exists an edge between the nodes that correspond to Copyright © 2022, Association for the Advancement of Artificial the physical qubits, (Herbert and Sengupta 2018). This edge Intelligence (www.aaai.org). All rights reserved. could be either unidirectional or bidirectional, i.e., CNOT
q7 c q7 q8 e q8 q9 b q9 a 1 1 2 f 3 1 2 f 3 4 5 q1 q1 q1 q1 d d q2 q2 q3 q2 q3 q3 q3 q2 e q4 q4 q7 q4 a q5 q5 q5 q5 q6 q6 q9 q6 c c q7 q4 q7 q7 b e q8 q8 q6 q8 b q9 q9 q8 q9 a (a) Quantum circuit (b) 1Decomposed 2 f 3 circuit 4 5 (c) Decomposed transformed circuit 1 2 3 f 4 q1 q1 Figure 1: An q1 example of qubit routing on a quantum d circuit for 3×3 grid architecture (Figure 2a). (a) For simplicity, the original quantum circuit consists d of two-qubit gate only q2 operations. (b) Decompositionq3of the original quantum circuit into series of q2 slices such that all the instructions present inq a slice can be executed in parallel.q The two-qubit gate operations: {d, e} (green) 3 2 comply withq3 the topology of the grid architecture whereas the operations:e {a, b, c, f } (red) do not comply with the topology (and therefore cannot be performed). Note qthat 4 the successive two-qubit gateq7operations on q3 → q4 (blue) are redundant and are notq4considered while routing. (c) Decomposition ofathe transformed quantum circuit we get after qubit routing. Four q5 q5 additional SWAP q5 gates are added that increased the circuit depth to 5, i.e., an overhead circuit depth of 2. The final qubit labels are represented at the end right side of the circuit. q6 The qubits that are not moved q9 (or swapped) are shown in brown ({q1 , q5 }), q6 while the rest of them are shown in blue. c c q7 q4 q7 b e q8 q6 q8 be more than that of the original circuit due to the insertion 1 E2 q9 2 E4 3 b q9 of additional qSWAP 8 gates. This comes from the definition of 1 E1E2 2 E3E4a 3 E5 the term depth in the context of quantum circuits. This can E1 4 E7E3 5 E9E5 1 62 f 3 4 5 be understood by decomposing a quantum circuit into series of individual slices, each of which contains a group of gate 1E7 5 E9 6 4qE6 d E8 E10 q1 operations that have no overlapping qubits, i.e., all the in- E6q27 E11 E8 8 E12 E109 q 3 structions present in a slice can be executed in parallel (Fig. 1b). The depth of the quantum circuit then refers to the min- 7qE11 3 8 E12 9 q 2 imum number of such slices the circuit can be decomposed e (a) 3×3 grid architecture into, i.e., the minimum amount of parallel executions needed q4 with edges (i.e. neighbor- a (b) IBMQX-20q7architecture rep- to execute the circuit. The goal is to minimize the overhead ingq5qubits) labelled resented as a graph q5 depth of the transformed circuit with respect to the original circuit. Figure q6 2: Examples of qubit connectivity q9 graphs for some This goal involves solving two subsequent problems of common quantum architecturesc (i) qubit allocation, which refers to the mapping of program q7 q4 b qubits to logic qubits, and (ii) qubit movement, which refers q8 q6 to routing qubits between different locations such that inter- can be performed either in one direction or in both direc- action can be made possible (Tannu and Qureshi 2019). In tions.q9In this work, we consider only theqbidirectional 8 case, this work, we focus on the latter problem of qubit movement while noting that the direction of a CNOT gate can be re- only and refer to it as qubit routing. However, it should be versed by sandwiching it between a pair of Hadamard gates noted that qubit allocation is also an important problem and acting on both control and target qubits (Garcia-Escartin and it can play an important role in minimizing the effort needed Chamorro-Posada 2011). to perform qubit movement. Given a target hardware topology D and a quantum circuit C, the task of qubit routing refers to transforming this quan- tum circuit by adding a series of SWAP gates such that all 2.2 Related Work its gate operations then satisfy the connectivity constraints The first major attraction for solving the qubit routing prob- of the target topology (Fig. 1c). Formally, for a routing algo- lem was the competition organized by IBM in 2018 to find rithm R, we can represent this process as follows: the best routing algorithm. This competition was won by R(C, D) → C 0 (1) Zulehner, Paler, and Wille (2019), for developing a Com- puter Aided Design-based (CAD) routing strategy. Since 0 Depth of C (transformed quantum circuit) will, in general, then, this problem has been presented in many different
ways. These include graph-based architecture-agnostic so- In every timestep t, QRoute starts by greedily schedul- lution by Cowtan et al. (2019), showing equivalence to the ing all the operations that are both current and local in P. travelling salesman problem by Paler, Zulehner, and Wille To evolve M, it then performs a Monte Carlo tree search (2021), machine learning based methods by Paler et al. (MCTS) to find an optimal set of SWAPs by the evaluation (2020), and heuristic approaches by Venturelli et al. (2017), metrics described in the Section 3.2 such that all operations Baioletti, Rasconi, and Oddi (2021), Chand et al. (2019), etc. in the current timestep put together form a parallelizable set, A reinforcement learning in a combinatorial action space so- i.e., a set of local operations such that no two operations in lution was proposed by Herbert and Sengupta (2018), which the set act on the same qubit. The number of states we can suggested used simulated annealing to search through the encounter in the action space explodes exponentially with combinatorial action space, aided by a Feed-Forward neu- the depth of our search, therefore an explicit search till the ral network to judge the long-term expected depth. This was circuit is done compiling is not possible. Therefore we cut further extended to use Double Deep Q-learning and priori- short our search at some shallow intermediate state, and use tized experience replay by Pozzi et al. (2020). a neural network to get its heuristic evaluation. Recently, Monte Carlo tree search (MCTS), a popular re- The following subsections describe in greater detail the inforcement learning algorithm (Browne et al. 2012) previ- working of the search and the heuristic evaluation. ously proven successful in a variety of domains like playing puzzle games such as Chess and Go (Silver et al. 2016), and 3.1 State and Action Space was used by Zhou, Feng, and Li (2020) to develop a qubit Definition 3.1 (State) It captures entire specification of the routing solution. state of compilation at some timestep t. Abstractly, it is de- scribed as: 2.3 Our Contributions st = (D, Mt , Pt , Lt ) (2) Our work demonstrates the use of MCTS on the task of where, D is the topology of target hardware, and Mt and Qubit Routing and presents state of the art results. Following Pt represents the current values of M and P respectively. are the novelties of this approach: Lt is the set of nodes that are locked by the gate operations • We use an array of mutex locks to represent the current from the previous timestep and therefore cannot be operated state of parallelization, helping to reduce the depth of the in the current timestep. circuits instead of the total quantum volume, in contrast Definition 3.2 (Action) It is a set of SWAP gates (repre- to previous use of MCTS for qubit routing in Zhou, Feng, sented by the pair of qubits it acts on) such that all gates and Li (2020). are local, and its union with the set of operations that were • The actions in each timestep (layer of the output circuit) scheduled in the same timestep forms a parallelizable set. belong to a innumerably large action space. We phrase We are performing a tree search over state-action pairs. the construction of such actions as a Markov decision Since the number of actions that can be taken at any timestep process, making the training stabler and the results better, is exponential in the number of connections on the hardware, particularly at larger circuit sizes, than those obtained by we are forced to build a single action up, step-by-step. performing simulated annealing to search in such action Definition 3.3 (Move) It is a single step in a search proce- spaces (Herbert and Sengupta 2018; Pozzi et al. 2020). dure which either builds up the action or applies it to the Such approach should be applicable to other problems of current state. Moves are of the following two types: a similar nature. 1. SWAP(n1 , n2 ): Inserts a new SWAP on nodes n1 and n2 • Graph neural networks are used as an improved architec- into the action set. Such an insertion is only possible if ture to help guide the tree search. the operation is local and resulting set of operations for Finally, we provide a simple python package containing the timestep form a parallel set. the implementation of QRoute, together with an easy inter- 2. COMMIT: Finishes the construction of the action set for face for trying out different neural net architectures, com- that timestep. It also uses the action formed until now bining algorithms, reward structures, etc. to update the state st (schedules the SWAP gates on the hardware), and resets the action set for the next step. 3 Method In reality, different gate operations take different counts The QRoute algorithm takes in an input circuit and an in- of timesteps for execution. For example, if a hardware re- jective map, M : Q → N , from logical qubits to nodes quires SWAP gate to be broken down into CNOT gates, (physical qubits). Iteratively, over multiple timesteps, it tries then it would take three timesteps for complete execution to schedule the gate operations that are present in the input (Garcia-Escartin and Chamorro-Posada 2011). This means, circuit onto the target hardware. To do so, from the set of operations which are being scheduled must maintain mu- unscheduled gate operations, P, it takes all the current op- tual exclusivity with other other operations over the nodes erations, which are the first unscheduled operation for both which participates in them. This is essential to minimizing the qubits that they act on, and tries to make them into lo- the depth of the circuit since it models paralleizability of op- cal operations, which are those two-qubit operations that in- erations. volve qubits that are mapped to nodes connected on the tar- However, constructing a parallelizable set and represent- get hardware. ing the state of parallelization to our heuristic evaluator is a
1-3 2 3-1 1-3 2 3-1 4-6 5-9 9-5 4-6 5-9 9-5 7-4 8-6 6-8 7-4 8-6 6-8 X X X X 1 X X 1 X X 1 X X 1 X X 1 1 X X 1 1 X X 1 1 X X 1 1 X X 1-3 3-1 2 7-4 5 9 4-6 8 6-4 1-3 2 3-1 X Update 4-6 5-9 9-5 7-4 8-6 6-8 X X 1-3 3-1 2 7-4 5 9 Model 1 X X Update 1 X X 4-6 8 6-4 Evaluation X 1 1 X X Update 1 1 X X (value function) 1-3 3-1 2 Update 7-4 5 9 4-6 8 6-4 X Figure 3: Iteration of a Monte Carlo tree search: (i) select - recursively choosing a node in the search tree for exploration using a selection criteria, (ii) expand - expanding the tree using the available actions if an unexplored leaf node is selected, (iii) rollout - estimating long term reward using a neural network for the action-state pair for the new node added to search tree, and (iv) backup - propagating reward value from the leaf nodes upwards to update the evaluations of its ancestors. challenge. But an analogy can be drawn here to the nodes trated in Fig. 3. In each iteration, it begins traversing down being thought of as “resources” that cannot be shared, and the existing search tree by selecting the node with the maxi- the operations as “consumers” (Dijkstra 2002). This moti- mum UCT value (Eq. 4) at each level. During this traversal, vates us to propose the use of Mutex Locks for this purpose. whenever it encounters a leaf node, it expands the tree by These will lock a node until a scheduled gate operation in- choosing a move m from that leaf node. Then, it estimates volving that node executes completely. Therefore, this al- the scalar evaluation for the new state-action pair and back- lows our framework to naturally handle different types of propagates it up the tree to update evaluations of its ances- operations which take different amounts of time to complete. tors. For every state-action pair, the application of a feasible To build an optimal action set, we would want to select m move m on it will result in a new state-action pair: (s, a) −→ the move m with the maximum true Q-value. But since true (s0 , a0 ). This is a formulation of the problem of search as a Q-values are intractably expensive to compute, we can only Markov Decision Process. Associated with each such state- approximate the Q-values through efficient exploration. We action-move tuple ((s, a), m), we maintain two additional use the Upper Confidence Bound on Trees (UCT) objective values that are used by MCTS: (Kocsis and Szepesvári 2006) to balance exploration and ex- 1. N-value - The number of times we have taken the said ploitation as we traverse through the search tree. Moreover, move m from said state-action pair (s, a). as this problem results in a highly asymmetric tree, since some move block a lot of other moves, while others block 2. Q-value - Given a reward function R, it is the aver- fewer moves, we use the formulation of UCT adapted for age long-term reward expected after taking said move m asymmetric trees (Moerland et al. 2018): over all iterations of the search. (Future rewards are dis- counted by a factor γ) UCT((s, a), m) = Q((s, a), m) + Q((s, a), m) = R((s, a), m) + pP m N ((s, a), m) P 0 0 0 0 0 0 c × p(m|(s, a)) 0 N ((s , a ), m ) · Q((s , a ), m ) N ((s, a), m) γ m P 0 0 0 (4) m0 N ((s , a ), m ) (3) Here, the value p(m|(s, a)) is the prior policy function, which is obtained by adding a Dirichlet noise to the policy 3.2 Monte Carlo Tree Search output of the neural network (Silver et al. 2017). As MCTS Monte Carlo tree search progresses iteratively by executing continues probing the action space, it gets a better estimate its four phases: select, expand, rollout, and backup as illus- of the true values of the actions. This means that it acts as
Figure 4: Graph neural network architecture that approximates the value function and the policy function. a policy enhancement function whose output policy (Eq. 5) 2. Locked Edges - It is a set of edges (pairs of connected can be used to train the neural network’s prior (π), and the nodes) that are still locked due to either of its qubits be- average Q-value computed can be used to train its scalar ing involved in an operation in the current timestep or evaluation (Eq. 6). another longer operation that hasn’t yet terminated from the previous timesteps. π(m|(s, a)) ∝ N ((s, a), m)) (5) 3. Remaining targets - It is a list of the number of gate- P operations that are yet to be scheduled for each logical m Q((s, a), m) qubit. V((s, a)) = P (6) m N ((s, a), m) The SWAP operations each qubit would partake in de- The details of how MCTS progresses have been elabo- pends primarily on its target node, and on those of the nodes rated in the supplementary. Once it gets terminated, i.e., the in its neighborhood that might be competing for the same search gets completed, we go down the tree selecting the resources. It seems reasonable that we can use a Graph Neu- child with the maximum Q-value at each step until a COM- ral Network with the device topology graph for its connec- MIT action is found, we use the action set of the selected tivity since the decision of the optimal SWAP action for state-action pair to schedule SWAPs for the current timestep, some node is largely affected by other nodes in its physical and we re-root the tree at the child node of the COMMIT ac- neighborhood. Therefore, our architecture includes an edge- tion to prepare for the next timestep. convolution block (Wang et al. 2018), followed by some fully-connected layers with Swish (Ramachandran, Zoph, 3.3 Neural Network Architecture and Le 2017) activations for the policy and value heads. The Each iteration of the MCTS requires evaluation of Q-values value function and the policy function computed from this for a newly encountered state-action pair. But these values neural network are returned back to the MCTS. are impossible to be computed exactly since it would in- volve an intractable number of iterations in exploring and 4 Results expanding the complete search tree. Therefore, it is favor- We compare QRoute against the routing algorithms from able to heuristically evaluate the expected long-term reward other state-of-the-art frameworks on various circuit bench- from the state-action pair using a Neural Network, as it acts marks: (i) Qiskit and its three variants (Abraham and et al. as an excellent function approximator that can learn the sym- 2019): (a) basic, (b) stochastic, and (c) sabre, (ii) Deep-Q- metries and general rules inherent to the system. Networks (DQN) from (Pozzi et al. 2020), (iii) Cirq (Cirq So, once the MCTS sends a state-action pair to the eval- Developers 2021), and (iv) t|keti from Cambridge Quan- uator, it begins by committing the action to the state and tum Computing (CQC) (Sivarajah et al. 2020). Qiskit’s tran- getting the resultant state. We then generate the following spiler uses gate commutation rules while perform qubit rout- featurized representation of this state and pass this represen- ing. This strategy is shown to be advantageous in achieving tation through the neural-network architecture as shown in lower circuit depths (Itoko et al. 2019) but was disabled in Fig. 4. our simulations to have a fair comparison. The results for 1. Node Targets - It is a square boolean matrix whose rows DQN shown are adapted from the data provided by the au- and columns correspond to the nodes on a target device. thors Pozzi et al. (2020). An element (i, j) is true iff some logical qubits qx and qy are currently mapped to nodes i and j respectively, 4.1 Random Test Circuits such that (qx , qy ) is the first unscheduled operation that The first benchmark for comparing our performance com- qx partakes in. prises of random circuits. These circuits are generated on
QRoute 4000 2500 175 Cirq Qiskit (basic) 2000 Qiskit (stochastic) 3500 150 Qiskit (sabre) 1500 t|ket Best with QRoute Best w/o QRoute 3000 DQN Output Circuit Depth Output Circuit Depth 125 2500 100 2000 75 1500 50 1000 500 25 0 40 60 80 100 120 140 QRoute DQN t|ket Qiskit (basic) Qiskit (sabre) Qiskit (stochastic) Cirq Number of Input Gates Quantum Compiler Figure 5: Comparative performance of routing algorithms Figure 6: Plots of output circuit depths of routing algorithms on random circuits as a function of the number of two-qubit over small realistic circuits (≤ 100 gates), summed over the operations in the circuit. entire dataset. The inset shows the results on the same data comparing the best performant scheduler excluding and in- cluding QRoute on each circuit respectively. the fly and initialized with the same number of qubits as there are nodes on the device. Then two-qubit gates are put between any pair of qubits chosen at random. In our simu- On this dataset also, closest to QRoute performance is lations, the number of such gates is varied from 30 to 150 shown by Deep-Q-Network guided simulated annealer. To and the results for assessing performance of different frame- compare performances, we look at the average circuit depth works are given in Fig. 5. The experiments were repeated 10 ratio (CDR), which is defined by (Pozzi et al. 2020): times on each circuit size, and final results were aggregated over this repetition. 1 X Output Circuit Depth CDR = (7) Amongst the frameworks compared, QRoute ranks a very #circuits Input Circuit Depth circuits close second only to Deep-Q-Network guided simulated an- nealer (DQN). Nevertheless, QRoute still does consistently The resultant CDR for QRoute is 1.178, where as the re- better than all the other major frameworks: Qiskit, Cirq and ported CDR for the DQN is 1.19. In fact, QRoute outper- t|keti, and it scales well when we increase the number of lay- forms DQN on at least 80% of the circuits. This is signif- ers and the layer density in the input circuit. QRoute shows icant because in contrast to the random circuit benchmark, equivalent performance to DQN on smaller circuits, and on the realistic benchmarks consist of the circuits that are closer the larger circuits it outputs depths which are on average ≤ 4 to the circuits used in useful computation. layers more than those of DQN. Some part of this can be at- tributed to MCTS, in it’s limited depth search, choosing the 4.3 Large Realistic Circuit worse of two moves with very close Q-values, resulting in the scheduling of some unnecessary SWAP operations. For final benchmark, we take eight large circuits ranging from 154 gates to 5960 gates in its input from the IBM- 4.2 Small Realistic Circuits Quantum realistic test dataset (Zulehner, Paler, and Wille Next we test on the set of all circuits which use 100 or less 2018). The results are plotted in Fig. 7. QRoute has the gates from the IBM-Q realistic quantum circuit dataset used best performance of all available routing methods: Qiskit by Zulehner, Paler, and Wille (2018). The comparative per- and t|keti, on every one of these sampled circuits with on formance of all routing frameworks has been shown by plot- an average 13.6% lower circuit depth, and notable increase ting the depths of the output circuits summed over all the cir- in winning difference on the larger circuits. cuits in the test set in Fig. 6. Since the lack of a good initial The results from DQN and Cirq are not available for these qubit allocation becomes a significant problem for all pure benchmarks as they are not designed to scale to such huge routing algorithms on small circuits, we have benchmarked circuits. In case of DQN, the CDR data results were not pro- QRoute on this dataset from three trials with different initial vided for the circuits over 200 gates, mainly because sim- allocations. ulated annealing used in it is computationally expensive. The model presented herein has the best performance on Similarly, for Cirq, it takes several days to compile each of this dataset. We also compare the best result from a pool the near 5000 qubit circuits. In contrast, QRoute is able to of all routers including QRoute against that of another pool compile these circuits in at most 4 hours, and its compila- of the same routers but excluding QRoute. The pool includ- tion process can be sped up by reducing the depth of the ing QRoute gives on average 2.5% lower circuit depth, in- search. Spending more time, however, helps MCTS to bet- dicating that there is a significant number of circuits where ter approximate the Q-values leading to circuits with lower QRoute is the best routing method available. resulting depth.
9000 QRoute t|ket⟩ Qiskit (basic) 8000 Qiskit (stochastic) Qiskit (sabre) 7000 6000 Depth of Output Circuit 5000 4000 3000 2000 1000 0 154 gates 1343 gates 1405 gates 1498 gates 1701 gates 2100 gates 2319 gates 2648 gates 3089 gates 4459 gates 5960 gates Number of Gates in Circuit Figure 7: The results over eight circuits sampled from the large realistic dataset benchmark, the outputs of each routing algorithm are shown for every circuit. 5 Discussion and Conclusion ulated Annealing to build up their actions is that the algo- rithm cannot plan for the gates which are not yet waiting In this article, we have shown that the problem of qubit rout- to be scheduled, those which will come to the head of the ing has a very powerful and elegant formulation in Rein- list once the gates which are currently waiting are executed forcement Learning (RL) which can surpass the results of (Pozzi et al. 2020). QRoute also shares this deficiency, but any classical heuristic algorithm across all sizes of circuits the effect of this issue is mitigated by the explicit tree search and types of architectures. Furthermore, the central idea of which takes into account the rewards that will be accrued building up solutions step-by-step when searching in com- in the longer-term future. There is scope to further improve binatorial action spaces and enforcing constraints using mu- this by feeding the entire list of future targets directly into tex locks, can be adapted for several other combinatorial our neural network by using transformer encoders to han- optimization problems (Mazyavkina et al. 2020; Xing and dle the arbitrary length sequence data. This and other as- Tu 2020; Xu and Lieberherr 2019; Abe et al. 2019; Lat- pects of neural network design will be a primary facet of erre et al. 2018). Our approach is flexible enough to com- future explorations. Another means of improving the perfor- pile circuits of any size onto any device, from small ones mance would be to introduce new actions by incorporating like IBMQX20 with 20 qubits, to much larger hardware like use of BRIDGE gates (Itoko et al. 2019) and gate commu- Google Sycamore (results provided in supplementary) with tation rules (Garcia-Escartin and Chamorro-Posada 2011) 53 qubits (the Circuit Depth Ratio for small realistic circuits alongside currently used SWAP gates. The advantage of for- on Google Sycamore was 1.64). Also, it intrinsically deals mer is that it allows running CNOT gates on non-adjacent with hardware having different primitive instruction set, for qubit without permuting the ordering of the logical qubits; example on hardware where SWAP gates are not a primitive whereas, the latter would allow MCTS to recognize the re- and they get decomposed to 3 operations. QRoute enjoys dundancy in action space, making its exploration and selec- significant tunability; hyperparameters can be changed eas- tion more efficient. ily to alter the tradeoff between time taken and optimality of Finally, we provide an open-sourced access to our soft- decisions, exploration and exploitation, etc. ware library. It will allow researchers and developers to im- QRoute is a reasonably fast method, taking well under 10 plement variants of our methods with minimal effort. We minutes to route a circuit with under 100 operations, and hope that this will aid future research in quantum circuit at most 4 hours for those with upto 5000 operations, when transformations. For review we are providing, the codebase tested on a personal machine with an i3 processor (3.7 GHz) and a multimedia in the supplementary. and no GPU acceleration. Yet more can be desired in terms On the whole, the Monte Carlo Tree Search for building of speed. However, it is hard to achieve any significant im- up solutions in combinatorial action spaces has exceeded the provement without reducing the number of search iterations current state of art methods that perform qubit routing. De- and trading off a bit of performance. More predictive neural spite its success, we note that QRoute is a primitive imple- networks can help squeeze in better speeds. mentation of our ideas, and there is great scope of improve- One of the challenges of methods like DQN, that use Sim- ment in future.
Acknowledgements Herbert, S.; and Sengupta, A. 2018. Using Reinforcement Learning to find Efficient Qubit Routing Policies for De- A. S. and U. A. have contributed equally to the work pre- ployment in Near-term Quantum Computers. arXiv e-prints. sented in this manuscript. A. S. and U. A. would like to ac- knowledge the assistance from their colleague Bhuvanesh IBM. 2021. IBM Quantum. https://quantum-computing. Sridharan for his help in implementing and refining MCTS. ibm.com/. Additionally, they would also like to thank their colleagues Itoko, T.; Raymond, R.; Imamichi, T.; and Matsuo, A. Jai Bardhan, Kalp Shah for their suggestions. 2019. Optimization of Quantum Circuit Mapping using Gate Transformation and Commutation. arXiv e-prints. References Karalekas, P. J.; Tezak, N. A.; Peterson, E. C.; Ryan, C. A.; da Silva, M. P.; and Smith, R. S. 2020. A quantum-classical Abe, K.; Xu, Z.; Sato, I.; and Sugiyama, M. 2019. Solv- cloud platform optimized for variational hybrid algorithms. ing NP-Hard Problems on Graphs with Extended AlphaGo Quantum Sci. Technol., 5: 024003. Zero. arXiv e-prints. Kocsis, L.; and Szepesvari, C. 2006. Bandit Based Monte- Abraham, H.; and et al. 2019. Qiskit: An Open-source Carlo Planning. In ECML. Framework for Quantum Computing. Kocsis, L.; and Szepesvári, C. 2006. Bandit based monte- Arute, F.; Arya, K.; Babbush, R.; and et al. 2019. Quantum carlo planning. In European conference on machine learn- supremacy using a programmable superconducting proces- ing, 282–293. Springer. sor. Nature, 574(7779): 505–510. Laterre, A.; Fu, Y.; Khalil Jabri, M.; Cohen, A.-S.; Kas, D.; Baioletti, M.; Rasconi, R.; and Oddi, A. 2021. A Novel Hajjar, K.; Dahl, T. S.; Kerkeni, A.; and Beguir, K. 2018. Ant Colony Optimization Strategy for the Quantum Circuit Ranked Reward: Enabling Self-Play Reinforcement Learn- Compilation Problem. In EvoCOP, 1–16. ing for Combinatorial Optimization. arXiv e-prints. Bourassa, J. E.; Alexander, R. N.; Vasmer, M.; Patil, A.; Tz- Mazyavkina, N.; Sviridov, S.; Ivanov, S.; and Burnaev, E. itrin, I.; Matsuura, T.; Su, D.; Baragiola, B. Q.; Guha, S.; 2020. Reinforcement Learning for Combinatorial Optimiza- Dauphinais, G.; Sabapathy, K. K.; Menicucci, N. C.; and tion: A Survey. arXiv e-prints. Dhand, I. 2021. Blueprint for a Scalable Photonic Fault- Moerland, T. M.; Broekens, J.; Plaat, A.; and Jonker, C. M. Tolerant Quantum Computer. Quantum, 5: 392. 2018. Monte Carlo Tree Search for Asymmetric Trees. arXiv e-prints. Browne, C. B.; Powley, E.; Whitehouse, D.; Lucas, S. M.; Cowling, P. I.; Rohlfshagen, P.; Tavener, S.; Perez, D.; Munos, R. 2014. From Bandits to Monte-Carlo Tree Search: Samothrakis, S.; and Colton, S. 2012. A Survey of Monte The Optimistic Principle Applied to Optimization and Plan- Carlo Tree Search Methods. IEEE Transactions on Compu- ning. Found. Trends Mach. Learn., 7: 1–129. tational Intelligence and AI in Games, 4(1): 1–43. Paler, A.; Sasu, L. M.; Florea, A.; and Andonie, R. 2020. Machine Learning Optimization of Quantum Circuit Lay- Chand, S.; Singh, H. K.; Ray, T.; and Ryan, M. 2019. Rollout outs. arXiv:2007.14608. based heuristics for the quantum circuit compilation prob- lem. In 2019 IEEE Congress on Evolutionary Computation Paler, A.; Zulehner, A.; and Wille, R. 2021. NISQ circuit (CEC), 974–981. IEEE. compilation is the travelling salesman problem on a torus. Quantum Science and Technology, 6(2): 025016. Childs, A. M.; Schoute, E.; and Unsal, C. M. 2019. Circuit Pozzi, M. G.; Herbert, S. J.; Sengupta, A.; and Mullins, R. D. Transformations for Quantum Architectures. In 14th Con- 2020. Using Reinforcement Learning to Perform Qubit ference on the Theory of Quantum Computation, Communi- Routing in Quantum Compilers. arXiv e-prints. cation and Cryptography (TQC 2019), volume 135 of Leib- niz International Proceedings in Informatics (LIPIcs), 3:1– Preskill, J. 2018. Quantum computing in the NISQ era and 3:24. Schloss Dagstuhl–Leibniz-Zentrum fuer Informatik. beyond. Quantum, 2: 79. ISBN 978-3-95977-112-2. Ramachandran, P.; Zoph, B.; and Le, Q. V. 2017. Searching for Activation Functions. arXiv e-prints. Cirq Developers. 2021. Cirq. Silver, D.; Huang, A.; Maddison, C. J.; Guez, A.; Sifre, L.; Cowtan, A.; Dilkes, S.; Duncan, R.; Krajenbrink, A.; Sim- van den Driessche, G.; and et al. 2016. Mastering the game mons, W.; and Sivarajah, S. 2019. On the Qubit Rout- of Go with deep neural networks and tree search. Nature, ing Problem. In 14th Conference on the Theory of Quan- 529(7587): 484–489. tum Computation, Communication and Cryptography (TQC Silver, D.; Hubert, T.; Schrittwieser, J.; Antonoglou, I.; Lai, 2019), volume 135 of Leibniz International Proceedings in M.; Guez, A.; Lanctot, M.; Sifre, L.; Kumaran, D.; Graepel, Informatics (LIPIcs), 5:1–5:32. Schloss Dagstuhl–Leibniz- T.; Lillicrap, T. P.; Simonyan, K.; and Hassabis, D. 2017. Zentrum fuer Informatik. ISBN 978-3-95977-112-2. Mastering Chess and Shogi by Self-Play with a General Re- Dijkstra, E. W. 2002. Cooperating Sequential Processes, inforcement Learning Algorithm. arXiv e-prints. 65–138. New York, NY: Springer New York. ISBN 978-1- Sivarajah, S.; Dilkes, S.; Cowtan, A.; Simmons, W.; Edging- 4757-3472-0. ton, A.; and Duncan, R. 2020. t|keti: A retargetable compiler Garcia-Escartin, J. C.; and Chamorro-Posada, P. 2011. for NISQ devices. Quantum Science and Technology, 6(1): Equivalent Quantum Circuits. arXiv e-prints. 014003.
Tannu, S. S.; and Qureshi, M. K. 2019. Not All Qubits Are Created Equal: A Case for Variability-Aware Policies for NISQ-Era Quantum Computers. In Proceedings of the Twenty-Fourth International Conference on Architectural Support for Programming Languages and Operating Sys- tems, ASPLOS ’19, 987–999. New York, NY, USA: Associ- ation for Computing Machinery. ISBN 9781450362405. Venturelli, D.; Do, M.; Rieffel, E. G.; and Frank, J. 2017. Temporal Planning for Compilation of Quantum Approxi- mate Optimization Circuits. In IJCAI, 4440–4446. Wang, Y.; Sun, Y.; Liu, Z.; Sarma, S. E.; Bronstein, M. M.; and Solomon, J. M. 2018. Dynamic Graph CNN for Learn- ing on Point Clouds. arXiv e-prints. Xing, Z.; and Tu, S. 2020. A Graph Neural Network Assisted Monte Carlo Tree Search Approach to Traveling Salesman Problem. IEEE Access, 8: 108418–108428. Xu, R.; and Lieberherr, K. 2019. Learning Self-Game-Play Agents for Combinatorial Optimization Problems. arXiv e- prints. Zhou, X.; Feng, Y.; and Li, S. 2020. A Monte Carlo Tree Search Framework for Quantum Circuit Transforma- tion. In Proceedings of the 39th International Confer- ence on Computer-Aided Design, ICCAD ’20. New York, NY, USA: Association for Computing Machinery. ISBN 9781450380263. Zulehner, A.; Paler, A.; and Wille, R. 2018. IBM Qiskit Developer Challenge. https://www.ibm.com/blogs/research/ 2018/08/winners-qiskit-developer-challenge/. Accessed on: 2021-03-23. Zulehner, A.; Paler, A.; and Wille, R. 2019. An Efficient Methodology for Mapping Quantum Circuits to the IBM QX Architectures. IEEE Transactions on Computer-Aided De- sign of Integrated Circuits and Systems, 38(7): 1226–1236.
You can also read