Conic programs arising in physics, quantum information, machine learning, and engineering are often defined over sparse graphs. Although such problems can be solved in polynomial time using classical interior-point solvers, the computational complexity scales unfavorably with graph size. We propose a variational quantum paradigm for solving conic programs, including quadratically constrained quadratic programs and semidefinite programs. We encode primal variables via the state of a parameterized quantum circuit (PQC) and dual variables via the probability mass function associated with a second PQC. The Lagrangian function can thus be expressed as scaled expectations of quantum observables. We pursue approximately stationary points of the Lagrangian by minimizing/maximizing the Lagrangian over the parameters of the first/second PQC. This is accomplished in a hybrid fashion: gradients are estimated using the two PQCs, while their parameters are updated classically using a primal-dual method. We propose permuting primal variables so that related observables have a banded form, enabling efficient measurement. We provide a complexity analysis that is useful to determine which problem types may enjoy quantum advantage. The framework is applied to the AC OPF problem, a large-scale optimization problem central to electric power system operation. Numerical tests on the IEEE 57-node system using PennyLane's simulator show that the proposed doubly variational quantum framework can find high-quality OPF solutions. While this demonstration does not yield a quantum speedup, the results serve as a proof-of-concept and highlight challenges toward practical quantum advantage. Although showcased for OPF, the framework has broader scope, including conic programs with many variables and constraints, problems defined over sparse graphs, and training quantum machine learning models to satisfy constraints.
Conic programs arising in physics, quantum information, machine learning, and engineering are often defined over sparse graphs. Although such problems can be solved in polynomial time using classical interior-point solvers, the computational complexity scales unfavorably with graph size. We propose a variational quantum paradigm for solving conic programs, including quadratically constrained quadratic programs and semidefinite programs. We encode primal variables via the state of a parameterized quantum circuit (PQC) and dual variables via the probability mass function associated with a second PQC. The Lagrangian function can thus be expressed as scaled expectations of quantum observables. We pursue approximately stationary points of the Lagrangian by minimizing/maximizing the Lagrangian over the parameters of the first/second PQC. This is accomplished in a hybrid fashion: gradients are estimated using the two PQCs, while their parameters are updated classically using a primal-dual method. We propose permuting primal variables so that related observables have a banded form, enabling efficient measurement. We provide a complexity analysis that is useful to determine which problem types may enjoy quantum advantage. The framework is applied to the AC OPF problem, a large-scale optimization problem central to electric power system operation. Numerical tests on the IEEE 57-node system using PennyLane's simulator show that the proposed doubly variational quantum framework can find high-quality OPF solutions. While this demonstration does not yield a quantum speedup, the results serve as a proof-of-concept and highlight challenges toward practical quantum advantage. Although showcased for OPF, the framework has broader scope, including conic programs with many variables and constraints, problems defined over sparse graphs, and training quantum machine learning models to satisfy constraints.
Quantum computing provides alternative algorithmic approaches to solve computationally challenging problems and has shown potential advantages compared to classical methods for specific problem classes. This thesis explores variational quantum algorithms (VQAs) for solving the alternating current optimal power flow (AC OPF) problem and the AC power flow (PF) problem. Both tasks are fundamental to power system analysis, operation, and planning, but they become increasingly challenging as grid size and complexity grow. The first main contribution of the thesis introduces a variational quantum framework for solving the OPF problem, formulated as a nonconvex quadratically constrained quadratic program. The proposed method encodes primal and dual optimization variables using two parameterized quantum circuits (PQCs), allowing the OPF Lagrangian to be expressed as expectations over quantum observables. An OPF solution can be found by minimizing the Lagrangian over the parameters of the first PQC and maximizing it over the parameters of the second PQC. The proposed framework seeks saddle points of the Lagrangian in a hybrid fashion. The gradients of the Lagrangian are estimated using the two PQCs, and PQC parameters are updated classically using a primal-dual method. A key innovation is to permute the primal variables of the original problem, so the corresponding quantum observables adhere to a banded form that facilitates a novel qubit-efficient quantum measurement protocol. The second main contribution presents a quantum machine learning (QML) approach for learning AC PF solutions using a data-embedded PQC. The method first poses a single PF instance as a nonlinear least-squares problem over PQC parameters. Second, the PF specifications are embedded as features into another PQC, and the resultant QML model is trained to learn the PF mapping. Once presented with a new problem instance, the trained QML model can perform the prediction step efficiently in real-time operation. Third, a novel protocol is devised to efficiently measure PF observables by leveraging the graph structure of a power network. The proposed quantum AC PF/OPF framework sets the foundations for addressing more elaborate grid tasks via quantum computing. More broadly, the proposed quantum paradigms can be used to achieve safe QML as PQCs are trained to comply with data-dependent constraints.
Variational quantum approaches have shown great promise in finding near-optimal solutions to computationally challenging tasks, including solving optimization problems. Nonetheless, optimization problems with constraints may not have been handled in a disciplined fashion thus far. To address this gap, this work proposes a hybrid quantum-classical algorithmic paradigm termed the variational quantum eigensolver with constraints (VQEC) that extends the celebrated VQE to handle optimization with constraints. As with the standard VQE, the vector of optimization variables is captured by the state of a variational quantum circuit (VQC). To deal with constraints, VQEC optimizes a Lagrangian function classically over both the VQC parameters as well as the dual variables associated with constraints. To comply with the quantum setup, variables are updated via a perturbed primal-dual method leveraging the parameter shift rule. Among a wide gamut of potential applications, we showcase how VQEC can approximately solve quadratically constrained binary optimization problems, find stochastic binary policies satisfying quadratic constraints on the average and in probability, and solve large-scale linear programs over the probability simplex. Under an assumption on the error for the VQC to approximate an arbitrary probability mass function, we provide bounds on the optimality gap attained by a VQC. Numerical tests on a quantum simulator investigate the effect of various parameters and corroborate that VQEC can generate high-quality solutions.
Variational quantum approaches have shown great promise in finding near-optimal solutions to computationally challenging tasks. Nonetheless, enforcing constraints in a disciplined fashion has been largely unexplored. To address this gap, this work proposes a hybrid quantum-classical algorithmic paradigm termed VQEC that extends the celebrated VQE to handle optimization with constraints. As with the standard VQE, the vector of optimization variables is captured by the state of a variational quantum circuit (VQC). To deal with constraints, VQEC optimizes a Lagrangian function classically over both the VQC parameters as well as the dual variables associated with constraints. To comply with the quantum setup, variables are updated via a perturbed primal-dual method leveraging the parameter shift rule. Among a wide gamut of potential applications, we showcase how VQEC can approximately solve quadratically-constrained binary optimization (QCBO) problems, find stochastic binary policies satisfying quadratic constraints on the average and in probability, and solve large-scale linear programs (LP) over the probability simplex. Under an assumption on the error for the VQC to approximate an arbitrary probability mass function (PMF), we provide bounds on the optimality gap attained by a VQC. Numerical tests on a quantum simulator investigate the effect of various parameters and corroborate that VQEC can generate high-quality solutions.
Alternating Current Optimal Power Flow Unit Commitment (AC-OPF-UC) is a difficult mixed-integer nonlinear optimization problem that combines binary generator commitment decisions with nonconvex continuous AC power-flow constraints. In this work, we investigate whether hybrid quantum-classical variational algorithms can improve the solution of single-period AC-OPF-UC relative to classical approaches. To the best of our knowledge, this is the first study to directly evaluate quantum or hybrid quantum-classical algorithms for the full AC-OPF-UC problem. We consider two candidate algorithms for improving AC-OPF-UC solution quality relative to purely classical methods on ideal quantum hardware. The first applies QAOA directly to a fully discretized formulation of the problem, with equality and inequality constraints incorporated through penalty terms and slack variables. Although conceptually straightforward, this approach requires a prohibitively large number of qubits even for small instances. The second, qubit-efficient approach encodes only the binary generator status variables on a quantum computer, while optimizing the continuous power-flow variables classically for each sampled bitstring. We benchmark this method on randomly generated AC-OPF-UC instances with 5 to 13 generators and compare it against SCIP, SMAC, and uniform random sampling. Our simulations show that the qubit-efficient hybrid method does not outperform uniform sampling. These results suggest that in order to establish potential advantage of the variational hybrid strategy considered here over the best classical algorithms, if any, much larger system sizes (25+ generators) need to be tested, which is beyond our computational capacity. Alternatively, different approaches, such as quantum versions of branch-and-bound methods, may be more promising.
Optimal Power Flow (OPF) is a fundamental optimisation problem in power systems, aiming to minimize the generation cost dispatch while satisfying operational constraints. Traditional optimisation methods, especially the AC-OPF formulation, struggle to converge or find optimal solutions due tothenon-linear andnon-convexnatureoftheproblem. Thisstudyaddressesthe AC-OPF using a hybrid iterative algorithm that combines classical and quantum computing, where Quantum Annealing assists in optimizing decision variables. ThemethodologysplitstheAC-OPFintotwointerrelatedproblems. Thefirstsub-problemisthe classical OPFformulation, whichincorporatesallelectricalnetworkcomponentsandminimizes the objective function under technical and operational constraints. The second sub-problem, also known as the master problem, is a binary problem that assumes activation states of the generators, while minimizingcostsofgenerationandstart-upcostsfordifferenttimesteps. This forms aBinary Integer Linear Programming (BILP) problem, which is well-suited for quantum computing due to its discrete structure and combinatorial complexity. Once the BILP is defined, it is converted into a Hamiltonian and evolved using Quantum Annealing throughtheQiboframework. Theresultingoptimalbinarycombinationisthenapplied to the network model and solved classically. If the solution is feasible and optimal, the process ends. If not, the combination is excluded from the BILP and added as a constraint to guide the search for the next best feasible solution. This methodology is tested on two IEEE networks of 9 to evaluate computational speedup and assess the limitations of actual quantum hardware. Once this is done, future work will explore the use of Quantum Machine Learning to generate 24-hour load and generation profiles, while in this work, an LSTM is used, enabling the AC-OPF solutions.
Abstract Quadratically Constrained Quadratic Programs (QCQPs) are an important class of optimization problems with diverse real-world applications. In this work, we propose a variational quantum algorithm for general QCQPs. By encoding the variables in the amplitude of a quantum state, the requirement for the qubit number scales logarithmically with the dimension of the variables, which makes our algorithm suitable for current quantum devices. Using the primal-dual interior-point method in classical optimization, we can deal with general quadratic constraints. Our numerical experiments on typical QCQP problems, including Max-Cut and optimal power flow problems, demonstrate better performance of our hybrid algorithm over classical counterparts.