We investigate the performance of different annealers for power flow analysis. The annealers are D-Wave’s simulated annealer Neal, D-Wave’s quantum-classical hybrid annealer, D-Wave’s Advantage™ system (QA), Fujitsu’s simulated annealer, and Fujitsu’s digital annealer V3 (DA). We implement Quadratic Unconstrained Binary Optimization (QUBO) and Ising model formulations, with the latter offering finer control over complex voltage adjustments. We also propose a new method for the iterative and adaptive updating of complex voltages for the QUBO and Ising model formulations. Six different test systems are experimented with. QA and DA show superior performance over classical annealers. DA effectively manages the 15-bus system, whereas QA encounters difficulties in embedding the problem graph onto the hardware graph because of the limited qubit connectivity. This constraint confines QA to the 14-bus system with the QUBO formulation and the 4-bus system with the Ising model formulation. The best performance is associated with different annealers across different test systems, which suggests that adjusting the threshold can improve precision if the compiler and annealer are capable of handling the variables involved.
We investigate the performance of different annealers for power flow analysis. The annealers are D-Wave’s simulated annealer Neal, D-Wave’s quantum-classical hybrid annealer, D-Wave’s Advantage™ system (QA), Fujitsu’s simulated annealer, and Fujitsu’s digital annealer V3 (DA). We implement Quadratic Unconstrained Binary Optimization (QUBO) and Ising model formulations, with the latter offering finer control over complex voltage adjustments. We also propose a new method for the iterative and adaptive updating of complex voltages for the QUBO and Ising model formulations. Six different test systems are experimented with. QA and DA show superior performance over classical annealers. DA effectively manages the 15-bus system, whereas QA encounters difficulties in embedding the problem graph onto the hardware graph because of the limited qubit connectivity. This constraint confines QA to the 14-bus system with the QUBO formulation and the 4-bus system with the Ising model formulation. The best performance is associated with different annealers across different test systems, which suggests that adjusting the threshold can improve precision if the compiler and annealer are capable of handling the variables involved.
We investigate the performance of different annealers for power flow analysis using adiabatic computing. The annealers include D-Wave's simulated annealer Neal, D-Wave's quantum-classical hybrid annealer, D-Wave's Advantages system (QA), Fujitsu's classical simulated annealer, and Fujitsu's digital annealer V3 (DA). We implement Quadratic Unconstrained Binary Optimization (QUBO) and Ising model formulations, with the latter offering finer control over complex voltage adjustments. Different test systems are experimented with to systematically evaluate the annealers. The evaluation is based on the accuracy, the annealer's capability to handle the decision variables, and the computational time needed. QA and DA show superior performance over classical annealers for our application. DA effectively manages larger test systems, whereas QA encounters difficulties embedding the problem graph onto the hardware graph because of the limited qubit connectivity. This constraint confines QA to the 14-bus system with the QUBO formulation and the 4-bus system with the Ising model formulation. The best performance is associated with different annealers across different test systems, which suggests that adjusting the threshold can improve precision if the compiler and annealer are capable of handling the number of variables involved.
Power flow (PF) analysis is a foundational computational method to study the flow of power in an electrical network. This analysis involves solving a set of non-linear and non-convex differential-algebraic equations. State-of-the-art solvers for PF analysis, therefore, face challenges with scalability and convergence, specifically for large-scale and/or ill-conditioned cases characterized by high penetration of renewable energy sources, among others. The adiabatic quantum computing paradigm has been proven to efficiently find solutions for combinatorial problems in the noisy intermediate-scale quantum (NISQ) era, and it can potentially address the limitations posed by state-of-the-art PF solvers. For the first time, we propose a novel adiabatic quantum computing approach for efficient PF analysis. Our key contributions are (i) a combinatorial PF algorithm and a modified version that aligns with the principles of PF analysis, termed the adiabatic quantum PF algorithm (AQPF), both of which use Quadratic Unconstrained Binary Optimization (QUBO) and Ising model formulations; (ii) a scalability study of the AQPF algorithm; and (iii) an extension of the AQPF algorithm to handle larger problem sizes using a partitioned approach. Numerical experiments are conducted using different test system sizes on D-Wave's Advantage™ quantum annealer, Fujitsu's digital annealer V3, D-Wave's quantum-classical hybrid annealer, and two simulated annealers running on classical computer hardware. The reported results demonstrate the effectiveness and high accuracy of the proposed AQPF algorithm and its potential to speed up the PF analysis process while handling ill-conditioned cases using quantum and quantum-inspired algorithms.
This study further explores reformulating power flow (PF) analysis as a discrete combinatorial optimization problem, proposed in our earlier study using the Adiabatic Quantum Power Flow (AQPF) algorithm, which can be executed on Ising machines, including quantum and quantum-inspired hardware. This approach provides a new representation of the underlying equations, analogous to how neural networks approximate complex functions using simple operations. While the resulting combinatorial optimization problem is NP-hard, it is compatible with emerging quantum hardware designed to address such complexity. We introduce the Adiabatic Quantum Optimal Power Flow (AQOPF) algorithm, which transforms the classical optimal power flow (OPF) equations into quadratic unconstrained binary optimization (QUBO) models. Furthermore, the AQPF and AQOPF algorithms are evaluated on standard test cases ranging from 4- to 1354-bus systems using D-Wave's Advantage\texttrademark\ system (QA), its hybrid quantum-classical solver (HA), and Fujitsu's third-generation Digital Annealer (DAv3) and Quantum-Inspired Integrated Optimization (QIIO) platform. Both full and partitioned formulations are investigated, with particular attention to scalability and robustness in ill-conditioned scenarios. The results demonstrate that the algorithms can reproduce feasible PF and OPF solutions and exhibit promising computational scalability when supported by scalable hardware.
Abstract Quadratic unconstrained binary optimization (QUBO) formulations of combinatorial optimization problems allow for solving them using various quantum heuristics. While large-scale quantum computations are currently still out of reach, we can already numerically test such QUBO formulations on a perhaps surprisingly large scale. In this work, we benchmark Fujitsu’s Digital Annealer (DA) on the Max-Cut problem, which captures the main complexity of the QUBO problem. We make a comprehensive benchmark against leading other heuristic algorithms on graphs with up to 53,000 variables by focusing on the wall-clock time. Moreover, we compare the DA performance against published performance results of the D-Wave hybrid quantum–classical annealer and the recently proposed QIS3 heuristic. Based on performance statistics for over 2000 graphs from the MQLib, we find that the DA yields competitive results. We hope that this benchmark demonstrates the extent to which large QUBO instances can be heuristically solved today, yielding consistent results across different solvers.
Combinatorial optimization problems are ubiquitous in various disciplines and applications. Many heuristic algorithms have been devoted to solve these types of problems. In order to increase the efficiency for finding the optimal solutions, an application-specific hardware, called digital annealer (DA) has been developed for solving combinatorial optimization problems using quadratic unconstrained binary optimization (QUBO) formulations. In this study, we formulated the number partitioning problem and the graph partitioning problem into QUBO forms and solved such problems with the DA developed by Fujitsu Ltd. The QUBO formulation of the number partitioning problem is fully connected. The DA found the overall runtime for the optimal solution to be less than 30 seconds for 6500 binary variables. For the graph partitioning problem, we adopted modularity as the metric for determining the quality of the partitions. For Zachary's Karate Club graph, the modularity obtained was 0.445, a 6% increase against D-wave Quantum Annealer and Simulated Annealing. Moreover, to explore the DA's potential applications to real-world problems, we used the search for communities or virtual microgrids in a power distribution network as an example. The problem was formulated into graph partitioning. It is shown that the DA effectively identified community structures in the IEEE 33-bus and IEEE 118-bus network.
QUBO formulations of combinatorial optimization problems allow for solving them using various quantum heuristics. While large-scale quantum computations are currently still out of reach, we can already numerically test such QUBO formulations on a perhaps surprisingly large scale. In this work, we benchmark Fujitsu’s Digital Annealer (DA) on the Max-Cut problem, which captures the main complexity of the QUBO problem. We make a comprehensive benchmark against leading other heuristic algorithms on graphs with up to 53,000 variables by focusing on the wall-clock time. Moreover, we compare the DA performance against published performance results of the D-Wave hybrid quantum-classical annealer and the recently proposed QIS3 heuristic. Based on performance statistics for over 2,000 graphs from the MQLib, we find that the DA yields competitive results. We hope that this benchmark demonstrates the extent to which large QUBO instances can be heuristically solved today, yielding consistent results across different solvers.