The quantum approximation optimisation algorithm (QAOA) is one of the leading
candidates for testing the applicability of gate-model quantum resources at solving
optimisation problems on small-sized quantum hardware. A combinatorial problem that
has been understudied with QAOA is correlation clustering: given a weighted (+1,-1)
graph, where the edge weights indicate whether two nodes are similar (positive edge
weight) or different (negative edge weight), the task in correlation clustering is to find a
clustering that either maximises agreements, minimises disagreements or a combination
of both. In this thesis, we design Hamiltonian formulations that encode correlation
clustering problems such that they can be solved with QAOA. For all Hamiltonian
formulations we propose circuit implementations and study their complexities. To
benchmark the performances of a basic QAOA algorithm using these formulations, we
use numerical simulations on complete graph data sets. For one of the formulations,
which uses a multi-level approach naturally suitable to qudit systems, we investigate
the performances of several optimisers and introduce heuristic strategies to further
improve its performance. On all instances in our data-sets, which include complete as
well as Erdős–Renyi graphs, the improved algorithm shows competitive performances
for QAOA depth p ≥ 2. We also show that for this algorithm at p = 1 parameters exists
such that it has a performance guarantee of 0.670 on 3-regular graphs.
The quantum approximation optimisation algorithm (QAOA) is one of the leading
candidates for testing the applicability of gate-model quantum resources at solving
optimisation problems on small-sized quantum hardware. A combinatorial problem that
has been understudied with QAOA is correlation clustering: given a weighted (+1,-1)
graph, where the edge weights indicate whether two nodes are similar (positive edge
weight) or different (negative edge weight), the task in correlation clustering is to find a
clustering that either maximises agreements, minimises disagreements or a combination
of both. In this thesis, we design Hamiltonian formulations that encode correlation
clustering problems such that they can be solved with QAOA. For all Hamiltonian
formulations we propose circuit implementations and study their complexities. To
benchmark the performances of a basic QAOA algorithm using these formulations, we
use numerical simulations on complete graph data sets. For one of the formulations,
which uses a multi-level approach naturally suitable to qudit systems, we investigate
the performances of several optimisers and introduce heuristic strategies to further
improve its performance. On all instances in our data-sets, which include complete as
well as Erdős–Renyi graphs, the improved algorithm shows competitive performances
for QAOA depth p ≥ 2. We also show that for this algorithm at p = 1 parameters exists
such that it has a performance guarantee of 0.670 on 3-regular graphs.
We study the correlation clustering problem using the quantum approximate optimization algorithm (QAOA) and qudits, which constitute a natural platform for such non-binary problems. Specifically, we consider a neutral atom quantum computer and propose a full stack approach for correlation clustering, including Hamiltonian formulation of the algorithm, analysis of its performance, identification of a suitable level structure for 87Sr and specific gate design. We show the qudit implementation is superior to the qubit encoding as quantified by the gate count. For single layer QAOA, we also prove (conjecture) a lower bound of 0.6367 ( 0.6699 ) for the approximation ratio on 3-regular graphs. Our numerical studies evaluate the algorithm's performance by considering complete and Erdős-Rényi graphs of up to 7 vertices and clusters. We find that in all cases the QAOA surpasses the Swamy bound 0.7666 for the approximation ratio for QAOA depths p≥2 . Finally, by analysing the effect of errors when solving complete graphs we find that their inclusion severely limits the algorithm's performance.
Quantum computing promises solutions to classically difficult and new-found problems through controlling the subtleties of quantum computing. The Quantum Approximate Optimisation Algorithm (QAOA) is a recently proposed quantum algorithm designed to tackle difficult combinatorial optimisation problems utilising both quantum and classical computation. The hybrid nature, generality and typically low gate-depth make it a strong candidate for near-term implementation in quantum computing. Finding the practical limits of the algorithm is currently an open problem. Until now, no tools to facilitate the design and validation of probabilistic quantum optimisation algorithms such as the QAOA on a non-trivial scale exist. Graph similarity is a long standing classically difficult problem withstanding decades of research from academia and industry. Determining the maximal edge overlap between all possible node label permutations is an NP-Complete task and provides an apt measure of graph similarity. We introduce a novel quantum optimisation simulation package facilitating investigation of all constituent components of the QAOA from desktop to cluster scale using graph similarity as an example. Our simulation provides flexibility and performance. We investigate eight classical optimisation methods each at six levels of decomposition. Moreover an encoding for permutation based problems such as graph similarity through edge overlap to the QAOA allows for significant quantum memory savings at the cost of additional operations. This compromise extends into the classical portion of the algorithm as the inclusion of infeasible solutions creates a challenging cost-function landscape. We present performance analysis of our simulation and of the QAOA setting a precedent for investigating and validating numerous other difficult problems to the QAOA as we move towards realising practical quantum computation.
The Quantum Approximate Optimization Algorithm (QAOA) is one of the promising near-term algorithms designed to find approximate solutions for combinatorial optimization problems. The algorithm prepares a parametrized state that is aimed to maximize the expectation value of the objective function of the problem. The circuit for QAOA consists of p layers, and it depends on 2p parameters, the determination of which can be carried out using a variational quantum eigensolver (VQE) subroutine. This is a variational approach for finding optimal parameters, where the quantum state is prepared with a quantum computer, and a classical computer is used for calculating the corresponding objective function value and performing an outer loop optimization routine in order to find optimal parameters. Earlier work on the algorithm has shown that QAOA is closely related to the Quantum Adiabatic Algorithm, which results in noticeable patterns in the optimal parameters. To this purpose, two methods, INTERP and FOURIER, were proposed to exploit this structure by Zhou et al. (2018). In this thesis, I extend upon these results by examining the effectiveness of one of their methods, the INTERP method, on the Max-Cut problem from graph theory. The performance of QAOA using this method was studied on several classes of graphs and benchmarked against the currently best-known classical algorithm Goemans-Williamson, which achieves an approximation ratio of ρ ≈ 0.878 on arbitrary graphs. The classes of graphs examined include cyclic, weighted and unweighted 3-regular, and Erdős–Rényi graphs with edge probabilities of 0.5 and 0.75 and graph sizes ranging from 4 to 16. Moreover, it was investigated if the INTERP method offers an advantage compared to random initialization of parameters and whether or not the method is polynomial in p. It was found that different graph classes lead to similar but different parameter patterns. Furthermore, the performance of QAOA using the INTERP method was dependent on the graph class considered. For the small graphs investigated it was found that QAOA outperforms Goemans-Williamson on ER-0.50, ER-0.75 and weighted 3-regular graphs from p=7 and beyond with respect to how well the expectation of the objective function approximates the optimal value, while for unweighted 3-regular graphs p=6 suffices. However, the calculation of the parameters is costly and requires a number of function evaluations that is on the order of the number of possible bipartitions of small graph with sizes n ≤ 16. Moreover, the numerical results suggested that indeed the method is polynomial in p advocating that the algorithm might be a strong alternative to Goemans-Williamson for large graph sizes, once hardware allows it.
Abstract The quantum approximate optimisation algorithm (QAOA) is a hybrid quantum–classical algorithm used to approximately solve combinatorial optimisation problems. It involves multiple iterations of a parameterised ansatz that consists of a problem and mixer Hamiltonian, with the parameters being classically optimised. While QAOA can be implemented on near-term quantum hardware, physical limitations such as gate noise, restricted qubit connectivity, and state-preparation-and-measurement (SPAM) errors can limit circuit depth and decrease performance. To address these limitations, this work introduces the eXpressive QAOA (XQAOA), an overparameterised variant of QAOA that assigns more classical parameters to the ansatz to improve its performance at low depths. XQAOA also introduces an additional Pauli-Y component in the mixer Hamiltonian, allowing the mixer to implement arbitrary unitary transformations on each qubit. To benchmark the performance of XQAOA at unit depth, we derive its closed-form expression for the MaxCut problem and compare it to QAOA, Multi-Angle QAOA (MA-QAOA) (Herrman et al 2022 Sci. Rep. 12 6781), a classical-relaxed algorithm, and the state-of-the-art Goemans–Williamson algorithm on a set of unweighted regular graphs with 128 and 256 nodes for degrees ranging from 3 to 10. Our results indicate that at unit depth, XQAOA has benign loss landscapes with local minima concentrated near the global optimum, allowing it to consistently outperform QAOA, MA-QAOA, and the classical-relaxed algorithm on all graph instances and the Goemans–Williamson algorithm on graph instances with degrees greater than 4. Small-scale simulations also reveal that unit-depth XQAOA invariably surpasses both QAOA and MA-QAOA on all tested depths up to five. Additionally, we find an infinite family of graphs for which XQAOA solves MaxCut exactly and analytically show that for some graphs in this family, special cases of XQAOA are capable of achieving a much larger approximation ratio than QAOA. Overall, XQAOA is a more viable choice for variational quantum optimisation on near-term quantum devices, offering competitive performance at low depths.
The transition from theoretical quantum computing to practical real-world applications heavily relies onhybrid quantum-classical architectures designed for NoisyIntermediate-Scale Quantum (NISQ) devices. Among these approaches, the Quantum Approximate Optimization Algorithm(QAOA) is a prominent strategy for addressing complex combinatorial optimization challenges, such as the NP-hard MaxCut problem. Instead of performing the entire computationon a quantum processor, QAOA utilizes a quantum circuit toevaluate a parameterized cost function, whereas a classical optimizer refines the parameters iteratively. This study presentsa rigorous empirical evaluation comparing the performanceof QAOA against established classical solvers on small-scaleMax-Cut instances. All experiments were implemented andexecuted using Python (Qiskit Aer) on a state-vector simulator.Using Erdős–Rényi random graphs and regular graphs with4 to 20 nodes, we benchmark QAOA at depths p = 2 andp = 3 against simulated annealing, the Goemans–WilliamsonSDP-based algorithm, and a linear-time greedy heuristic. Ourexperimental results demonstrate that low-depth QAOA implementations exhibit severe performance degradation, prohibitive computational overhead (O(2^n) scaling), and low convergencerates compared to mature classical algorithms. At 20 nodes,QAOA (p = 2) achieves only a 0.467 approximation ratio versus0.856 for Goemans–Williamson, with a 210× runtime slowdownover greedy heuristics
This thesis is devoted to four modern quantum optimisation algorithms and their applications. In the first chapter, we provide a complete Quantum Approximate Optimisation Algorithm (QAOA) pipeline for the discrete global minimum-variance portfolio (DGMVP) problem. We introduce an ansatz with efficient initial-state preparation, new hard mixing operators satisfying the constraints of the DGMVP problem, and scalable cost operators utilising binary encoding of the qubits. Benchmarking the method with the most effective optimisation strategies---dual annealing and unfrozen layerwise optimisation---indicates a scaling advantage over classical methods in finding the global optimum of the DGMVP problem. In the second chapter, we show an end-to-end quantum/quantum-inspired Goemans-Williamson semidefinite programming (SDP) method for relaxing Quadratic Unconstrained Binary Optimisation (QUBO) problems and extracting their rounding energy. Our method identifies a class of QUBO instances that is effectively solvable in polylogarithmic time relative to the problem dimension, under the assumptions of Pauli sparsification of the QUBO cost matrix, reducibility of the diagonal constraints, and efficient evaluation of the corresponding Gibbs-state expectation values. Further, we identify the commuting 1D Hamiltonian and the Kronecker graph as practically important instances. Our numerical experiments show that our method solves the SDP and QUBO with up to 250 variables for 1D Hamiltonian instances and achieves results competitive with the best classical solver for Kronecker graphs, giving strong evidence of the power of our algorithm. In the third chapter, we improve the current QAOA method with an iterative warm-start initial-state preparation method. We implement the t-th order statistic state as the warm-start initial state generated from each QAOA optimisation loop. We achieve an improved approximation ratio when solving the 3-regular MaxCut problem, and find that the worst-case approximation ratio converges to that of the best classical algorithm for p=1 QAOA. Additionally, our method achieves more favourable scaling than standalone QAOA, single-string warm-started QAOA, and a classical constrained-sampling approach when solving the DGMVP problem. In the fourth chapter, we provide a hard-constrained Grover adaptive search (GAS) method for fixed-sum quadratic binary optimisation problems, and an extension of the alternating direction method of multipliers (ADMM) framework for fixed-sum quartic binary optimisation. Our method guarantees at least a quadratic speed-up over the classical method, and exponentially reduces the gate count for solving fixed-sum binary optimisation when compared with the soft-constrained GAS in the fault-tolerant quantum computing regime. In summary, our contributions include efficient algorithms and novel frameworks that substantially improve performance and broaden the applicability of these quantum optimisation methods. The proposed methods exhibit favourable scaling behaviour and indicate potential quantum advantages in solving challenging optimisation problems, thereby motivating further research at the interface of quantum computing and optimisation.