-
Quantum Computing Inspired Iterative Refinement for Semidefinite Optimization
Authors:
Mohammadhossein Mohammadisiahroudi,
Brandon Augustino,
Pouya Sampourmahani,
Tamás Terlaky
Abstract:
Iterative Refinement (IR) is a classical computing technique for obtaining highly precise solutions to linear systems of equations, as well as linear optimization problems. In this paper, motivated by the limited precision of quantum solvers, we develop the first IR scheme for solving semidefinite optimization (SDO) problems and explore two major impacts of the proposed IR scheme. First, we prove…
▽ More
Iterative Refinement (IR) is a classical computing technique for obtaining highly precise solutions to linear systems of equations, as well as linear optimization problems. In this paper, motivated by the limited precision of quantum solvers, we develop the first IR scheme for solving semidefinite optimization (SDO) problems and explore two major impacts of the proposed IR scheme. First, we prove that the proposed IR scheme exhibits quadratic convergence toward an optimal solution without any assumption on problem characteristics. We also show that using IR with Quantum Interior Point Methods (QIPMs) leads to exponential improvements in the worst-case overall running time of QIPMs, compared to previous best-performing QIPMs. We also discuss how the proposed IR scheme can be used with classical inexact SDO solvers, such as classical inexact IPMs with conjugate gradient methods.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
Quantum Optimization: Potential, Challenges, and the Path Forward
Authors:
Amira Abbas,
Andris Ambainis,
Brandon Augustino,
Andreas Bärtschi,
Harry Buhrman,
Carleton Coffrin,
Giorgio Cortiana,
Vedran Dunjko,
Daniel J. Egger,
Bruce G. Elmegreen,
Nicola Franco,
Filippo Fratini,
Bryce Fuller,
Julien Gacon,
Constantin Gonciulea,
Sander Gribling,
Swati Gupta,
Stuart Hadfield,
Raoul Heese,
Gerhard Kircher,
Thomas Kleinert,
Thorsten Koch,
Georgios Korpas,
Steve Lenk,
Jakub Marecek
, et al. (21 additional authors not shown)
Abstract:
Recent advances in quantum computers are demonstrating the ability to solve problems at a scale beyond brute force classical simulation. As such, a widespread interest in quantum algorithms has developed in many areas, with optimization being one of the most pronounced domains. Across computer science and physics, there are a number of algorithmic approaches, often with little linkage. This is fur…
▽ More
Recent advances in quantum computers are demonstrating the ability to solve problems at a scale beyond brute force classical simulation. As such, a widespread interest in quantum algorithms has developed in many areas, with optimization being one of the most pronounced domains. Across computer science and physics, there are a number of algorithmic approaches, often with little linkage. This is further complicated by the fragmented nature of the field of mathematical optimization, where major classes of optimization problems, such as combinatorial optimization, convex optimization, non-convex optimization, and stochastic extensions, have devoted communities. With these aspects in mind, this work draws on multiple approaches to study quantum optimization. Provably exact versus heuristic settings are first explained using computational complexity theory - highlighting where quantum advantage is possible in each context. Then, the core building blocks for quantum optimization algorithms are outlined to subsequently define prominent problem classes and identify key open questions that, if answered, will advance the field. The effects of scaling relevant problems on noisy quantum devices are also outlined in detail, alongside meaningful benchmarking problems. We underscore the importance of benchmarking by proposing clear metrics to conduct appropriate comparisons with classical optimization techniques. Lastly, we highlight two domains - finance and sustainability - as rich sources of optimization problems that could be used to benchmark, and eventually validate, the potential real-world impact of quantum optimization.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
A quantum central path algorithm for linear optimization
Authors:
Brandon Augustino,
Jiaqi Leng,
Giacomo Nannicini,
Tamás Terlaky,
Xiaodi Wu
Abstract:
We propose a novel quantum algorithm for solving linear optimization problems by quantum-mechanical simulation of the central path. While interior point methods follow the central path with an iterative algorithm that works with successive linearizations of the perturbed KKT conditions, we perform a single simulation working directly with the nonlinear complementarity equations. Combining our appr…
▽ More
We propose a novel quantum algorithm for solving linear optimization problems by quantum-mechanical simulation of the central path. While interior point methods follow the central path with an iterative algorithm that works with successive linearizations of the perturbed KKT conditions, we perform a single simulation working directly with the nonlinear complementarity equations. Combining our approach with iterative refinement techniques, we obtain an exact solution to a linear optimization problem involving $m$ constraints and $n$ variables using at most $\mathcal{O} \left( (m + n) \text{nnz} (A) κ(\mathcal{M}) L \cdot \text{polylog} \left(m, n, κ(\mathcal{M}) \right) \right)$ elementary gates and $\mathcal{O} \left( \text{nnz} (A) L \right)$ classical arithmetic operations, where $ \text{nnz} (A)$ is the total number of non-zero elements found in the constraint matrix, $L$ denotes binary input length of the problem data, and $κ(\mathcal{M})$ is a condition number that depends only on the problem data.
△ Less
Submitted 7 November, 2023;
originally announced November 2023.
-
Improvements to Quantum Interior Point Method for Linear Optimization
Authors:
Mohammadhossein Mohammadisiahroudi,
Zeguan Wu,
Brandon Augustino,
Arriele Carr,
Tamás Terlaky
Abstract:
Quantum linear system algorithms (QLSA) have the potential to speed up Interior Point Methods (IPM). However, a major challenge is that QLSAs are inexact and sensitive to the condition number of the coefficient matrices of linear systems. This sensitivity is exacerbated when the Newton systems arising in IPMs converge to a singular matrix. Recently, an Inexact Feasible Quantum IPM (IF-QIPM) has be…
▽ More
Quantum linear system algorithms (QLSA) have the potential to speed up Interior Point Methods (IPM). However, a major challenge is that QLSAs are inexact and sensitive to the condition number of the coefficient matrices of linear systems. This sensitivity is exacerbated when the Newton systems arising in IPMs converge to a singular matrix. Recently, an Inexact Feasible Quantum IPM (IF-QIPM) has been developed that addresses the inexactness of QLSAs and, in part, the influence of the condition number using iterative refinement. However, this method requires a large number of gates and qubits to be implemented. Here, we propose a new IF-QIPM using the normal equation system, which is more adaptable to near-term quantum devices. To mitigate the sensitivity to the condition number, we use preconditioning coupled with iterative refinement to obtain better gate complexity. Finally, we demonstrate the effectiveness of our approach on IBM Qiskit simulators
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
Generating Linear, Semidefinite, and Second-order Cone Optimization Problems for Numerical Experiments
Authors:
Mohammadhossein Mohammadisiahroudi,
Ramin Fakhimi,
Brandon Augustino,
Tamás Terlaky
Abstract:
The numerical performance of algorithms can be studied using test sets or procedures that generate such problems. This paper proposes various methods for generating linear, semidefinite, and second-order cone optimization problems. Specifically, we are interested in problem instances requiring a known optimal solution, a known optimal partition, a specific interior solution, or all these together.…
▽ More
The numerical performance of algorithms can be studied using test sets or procedures that generate such problems. This paper proposes various methods for generating linear, semidefinite, and second-order cone optimization problems. Specifically, we are interested in problem instances requiring a known optimal solution, a known optimal partition, a specific interior solution, or all these together. In the proposed problem generators, different characteristics of optimization problems, including dimension, size, condition number, degeneracy, optimal partition, and sparsity, can be chosen to facilitate comprehensive computational experiments. We also develop procedures to generate instances with a maximally complementary optimal solution with predetermined optimal partition to generate challenging semidefinite and second-order cone optimization problems. Generated instances enable us to evaluate efficient interior-point methods for conic optimization problems.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
An Inexact Feasible Quantum Interior Point Method for Linearly Constrained Quadratic Optimization
Authors:
Zeguan Wu,
Mohammadhossein Mohammadisiahroudi,
Brandon Augustino,
Xiu Yang,
Tamás Terlaky
Abstract:
Quantum linear system algorithms (QLSAs) have the potential to speed up algorithms that rely on solving linear systems. Interior Point Methods (IPMs) yield a fundamental family of polynomial-time algorithms for solving optimization problems. IPMs solve a Newton linear system at each iteration to find the search direction, and thus QLSAs can potentially speed up IPMs. Due to the noise in contempora…
▽ More
Quantum linear system algorithms (QLSAs) have the potential to speed up algorithms that rely on solving linear systems. Interior Point Methods (IPMs) yield a fundamental family of polynomial-time algorithms for solving optimization problems. IPMs solve a Newton linear system at each iteration to find the search direction, and thus QLSAs can potentially speed up IPMs. Due to the noise in contemporary quantum computers, such quantum-assisted IPM (QIPM) only allows an inexact solution for the Newton linear system. Typically, an inexact search direction leads to an infeasible solution. In our work, we propose an Inexact-Feasible QIPM (IF-QIPM) and show its advantage in solving linearly constrained quadratic optimization problems. We also apply the algorithm to $\ell_1$-norm soft margin support vector machine (SVM) problems and obtain the best complexity regarding dependence on dimension. This complexity bound is better than any existing classical or quantum algorithm that produces a classical solution.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
Solving the semidefinite relaxation of QUBOs in matrix multiplication time, and faster with a quantum computer
Authors:
Brandon Augustino,
Giacomo Nannicini,
Tamás Terlaky,
Luis Zuluaga
Abstract:
Recent works on quantum algorithms for solving semidefinite optimization (SDO) problems have leveraged a quantum-mechanical interpretation of positive semidefinite matrices to develop methods that obtain quantum speedups with respect to the dimension $n$ and number of constraints $m$. While their dependence on other parameters suggests no overall speedup over classical methodologies, some quantum…
▽ More
Recent works on quantum algorithms for solving semidefinite optimization (SDO) problems have leveraged a quantum-mechanical interpretation of positive semidefinite matrices to develop methods that obtain quantum speedups with respect to the dimension $n$ and number of constraints $m$. While their dependence on other parameters suggests no overall speedup over classical methodologies, some quantum SDO solvers provide speedups in the low-precision regime. We exploit this fact to our advantage, and present an iterative refinement scheme for the Hamiltonian Updates algorithm of Brandão et al. (Quantum 6, 625 (2022)) to exponentially improve the dependence of their algorithm on precision. As a result, we obtain a classical algorithm to solve the semidefinite relaxation of Quadratic Unconstrained Binary Optimization problems (QUBOs) in matrix multiplication time. Provided access to a quantum read/classical write random access memory (QRAM), a quantum implementation of our algorithm exhibits a worst case running time of $\mathcal{O} \left(ns + n^{1.5} \cdot \text{polylog} \left(n, \| C \|_F, \frac{1}ε \right) \right)$.
△ Less
Submitted 11 May, 2023; v1 submitted 10 January, 2023;
originally announced January 2023.
-
Quantum Interior Point Methods for Semidefinite Optimization
Authors:
Brandon Augustino,
Giacomo Nannicini,
Tamás Terlaky,
Luis F. Zuluaga
Abstract:
We present two quantum interior point methods for semidefinite optimization problems, building on recent advances in quantum linear system algorithms. The first scheme, more similar to a classical solution algorithm, computes an inexact search direction and is not guaranteed to explore only feasible points; the second scheme uses a nullspace representation of the Newton linear system to ensure fea…
▽ More
We present two quantum interior point methods for semidefinite optimization problems, building on recent advances in quantum linear system algorithms. The first scheme, more similar to a classical solution algorithm, computes an inexact search direction and is not guaranteed to explore only feasible points; the second scheme uses a nullspace representation of the Newton linear system to ensure feasibility even with inexact search directions. The second is a novel scheme that might seem impractical in the classical world, but it is well-suited for a hybrid quantum-classical setting. We show that both schemes converge to an optimal solution of the semidefinite optimization problem under standard assumptions. By comparing the theoretical performance of classical and quantum interior point methods with respect to various input parameters, we show that our second scheme obtains a speedup over classical algorithms in terms of the dimension of the problem $n$, but has worse dependence on other numerical parameters.
△ Less
Submitted 7 September, 2023; v1 submitted 11 December, 2021;
originally announced December 2021.
-
Co-optimizing the Smart Grid and Electric Public Transit Bus System
Authors:
Mertcan Yetkin,
Brandon R. Augustino,
Alberto J. Lamadrid,
Lawrence V. Snyder
Abstract:
As climate change provides impetus for investing in smart cities, with electrified public transit systems, we consider electric public transportation buses in an urban area, which play a role in the power system operations in addition to their typical function of serving public transit demand. Our model considers a social planner, such that the transit authority and the operator of the electricity…
▽ More
As climate change provides impetus for investing in smart cities, with electrified public transit systems, we consider electric public transportation buses in an urban area, which play a role in the power system operations in addition to their typical function of serving public transit demand. Our model considers a social planner, such that the transit authority and the operator of the electricity system co-optimize the power system to minimize the total operational cost of the grid, while satisfying additional transportation constraints on buses. We provide deterministic and stochastic formulations to co-optimize the system. Each stochastic formulation provides a different set of recourse actions to manage the variable renewable energy uncertainty: ramping up/down the conventional generators, or charging/discharging of the transit fleet. We demonstrate the capabilities of the model and the benefit obtained via a coordinated strategy. We compare the efficacies of these recourse actions to provide additional managerial insights. We analyze the effect of different pricing strategies on the co-optimization. We also conduct congestion analysis in the power network, comparing our cooperative approach to a non-cooperative strategy when we assume electrified fleet sizes grow with greater battery capacities. Given the recent momentum towards building smarter cities and electrifying transit systems, our results provide policy directions towards a sustainable future. We test our models using modified MATPOWER case files and verify our results with different sized power networks. This study is motivated by a project with a large transit authority in California.
△ Less
Submitted 25 July, 2021; v1 submitted 14 December, 2020;
originally announced December 2020.