Quadratic Unconstrained Binary Optimization Research Articles

Quantum annealing-based computed tomography using variational approach for a real-number image reconstruction

Objective. Despite recent advancements in quantum computing, the limited number of available qubits has hindered progress in CT reconstruction. This study investigates the feasibility of utilizing quantum annealing-based computed tomography (QACT) with current quantum bit levels. Approach. The QACT algorithm aims to precisely solve quadratic unconstrained binary optimization problems. Furthermore, a novel approach is proposed to reconstruct images by approximating real numbers using the variational method. This approach allows for accurate CT image reconstruction using a small number of qubits. The study examines the impact of projection data quantity and noise on various image sizes ranging from 4 × 4 to 24 × 24 pixels. The reconstructed results are compared against conventional reconstruction algorithms, namely maximum likelihood expectation maximization (MLEM) and filtered back projection (FBP). Main result. By employing the variational approach and utilizing two qubits for each pixel of the image, accurate reconstruction was achieved with an adequate number of projections. Under conditions of abundant projections and lower noise levels, the image quality in QACT algorithm outperformed that of MLEM and FBP algorithms. However, in situations with limited projection data and in the presence of noise, the image quality in QACT was inferior to that in MLEM. Significance. This study developed the QACT reconstruction algorithm using the variational approach for real-number reconstruction. Remarkably, only 2 qubits were required for each pixel representation, demonstrating their sufficiency for accurate reconstruction.

Quantum and quantum-inspired optimization for an in-core fuel management problem

Operation management of nuclear power plants consists of several computationally hard problems. Searching for an in-core fuel loading pattern is among them. The main challenge of this combinatorial optimization problem is the exponential growth of the search space with a number of loading elements. Here we study a reloading problem in a Quadratic Unconstrained Binary Optimization (QUBO) form. Such a form allows us to apply various techniques, including quantum annealing, classical simulated annealing, and quantum-inspired algorithm in order to find fuel reloading patterns for several realistic configurations of nuclear reactors. We present the results of benchmarking the in-core fuel management problem in the QUBO form using the aforementioned computational techniques. This work demonstrates potential applications of quantum computers and quantum-inspired algorithms in the energy industry.

Quantum and quantum-inspired optimization for solving the minimum bin packing problem

Quantum computing devices are believed to be powerful in solving hard computational tasks, in particular, combinatorial optimization problems. In the present work, we consider a particular type of the minimum bin packing problem, which can be used for solving the problem of filling spent nuclear fuel in deep-repository canisters that is relevant for atomic energy industry. We first redefine the aforementioned problem it in terms of quadratic unconstrained binary optimization. Such a representation is natively compatible with existing quantum annealing devices as well as quantum-inspired algorithms. We then present the results of the numerical comparison of quantum and quantum-inspired methods. Results of our study indicate on the possibility to solve industry-relevant problems of atomic energy industry using quantum and quantum-inspired optimization.

Quadratic Unconstrained Binary Optimization Approach for Incorporating Solvency Capital into Portfolio Optimization

In this paper, we consider the inclusion of the solvency capital requirement (SCR) into portfolio optimization by the use of a quadratic proxy model. The Solvency II directive requires insurance companies to calculate their SCR based on the complete loss distribution for the upcoming year. Since this task is, in general, computationally challenging for insurance companies (and therefore, not taken into account during portfolio optimization), employing more feasible proxy models provides a potential solution to this computational difficulty. Here, we present an approach that is also suitable for future applications in quantum computing. We analyze the approximability of the solvency capital ratio in a quadratic form using machine learning techniques. This allows for an easier consideration of the SCR in the classical mean-variance analysis. In addition, it allows the problem to be formulated as a quadratic unconstrained binary optimization (QUBO), which benefits from the potential speedup of quantum computing. We provide a detailed description of our model and the translation into a QUBO. Furthermore, we investigate the performance of our approach through experimental studies.

Multi-satellite cooperative scheduling method for large-scale tasks based on hybrid graph neural network and metaheuristic algorithm

Traditional heuristic optimization algorithms are no longer applicable to the multiple satellites scheduling with large-scale tasks, as they are unable to provide satisfactory performance in terms of convergence and speed when addressing complex constraint conflicts. Therefore, we propose a multi-satellite cooperative scheduling method for large-scale tasks based on a hybrid graph neural network (GNN) and metaheuristic algorithm. By using the representation and extraction capability of GNN for relations in graph, the features of large-scale tasks and their constraint relations were expressed, and the generalized knowledge was extracted. And combined with metaheuristic algorithm, the optimization framework for large-scale satellite task scheduling based knowledge and data was implemented. The method consists of a GNN pre-solver module for static constraint conflicts and a metaheuristic optimization module for dynamic constraint conflicts of satellite missions. In the GNN pre-solver module, the graph sample and aggregate network was used to learn and extract the common knowledge in the static task-conflict graph, and provided effective prior information for the metaheuristic optimization module. The quadratic unconstrained binary optimization loss function was designed with multiple influencing factors, to convert the discrete optimization function into a continuous function; a greedy threshold was also added to improve the task completion rate. Finally, numerical experimental results showed that the proposed method can achieve an efficient solution to the multi-satellite scheduling problem with tens of thousands tasks. Compared with the commonly used multi-satellite scheduling algorithms (the genetic algorithm (GA) and greedy-GA), the proposed method can obtain higher-quality solutions under the same conditions and greatly improve the computational efficiency of large-scale mission planning.

D-Wave Implementation of Quantum Annealing for Optimal Resource Allocation in Disaster Response Operation of Marikina City

Quantum computing shows a positive approach for addressing optimization challenges in NP-hard problems such as the vehicle routing problem (VRP). This study focuses on improving the efficiency of disaster response operations by localizing the application of D-wave quantum annealing in Marikina City. This study uses the Solution Partitioning Solver (SPS) and the Quadratic Unconstrained Binary Optimization (QUBO) formulation to convert the VRP into an equation that can be solved using quantum annealing. The study demonstrates that quantum computing effectively distributes resources during emergency response operations and improves overall operational efficiency. In determining the most effective route for each vehicle, the D-wave Leap API and QUBO representation compute the distances traveled by each vehicle. These findings contribute to the practical applications of quantum computing to revolutionize various fields, including disaster management. Implementing D-wave quantum annealing in Marikina City shows relevance for future advancements in optimizing resource allocation and improving disaster response operations.

Quantum computer-aided job scheduling for storage and retrieval systems

Abstract In this paper, a quantum computer-aided approach to job scheduling for automated storage and retrieval systems is introduced. The approach covers application cases, where various objects need to be transported between storage positions and the order of transport operations can be freely chosen. The objective of job scheduling is to arrange the transport operations in a sequence, where the cumulative costs of the transport operations and empty runs between subsequent transport operations are minimized. The scheduling problem is formulated as an asymmetric quadratic unconstrained binary optimization (QUBO) problem, in which the transport operations are modeled as nodes and empty runs are modeled as edges, with costs assigned to each node and each edge. An Quantum Approximate Optimization Algorithm (QAOA) is used to solve the QUBO. Evaluations of the quantum computer-aided job scheduling approach have been conducted on the IBM Q System One quantum computer in Ehningen. In particular, the running time for the solution of the QUBO has been investigated, as well as the scalability of the approach with respect to the required number of qubits.

Pneumonia detection by binary classification: classical, quantum, and hybrid approaches for support vector machine (SVM)

Early diagnosis of pneumonia is crucial to increase the chances of survival and reduce the recovery time of the patient. Chest X-ray images, the most widely used method in practice, are challenging to classify. Our aim is to develop a machine learning tool that can accurately classify images as belonging to normal or infected individuals. A support vector machine (SVM) is attractive because binary classification can be represented as an optimization problem, in particular as a Quadratic Unconstrained Binary Optimization (QUBO) model, which, in turn, maps naturally to an Ising model, thereby making annealing—classical, quantum, and hybrid—an attractive approach to explore. In this study, we offer a comparison between different methods: (1) a classical state-of-the-art implementation of SVM (LibSVM); (2) solving SVM with a classical solver (Gurobi), with and without decomposition; (3) solving SVM with simulated annealing; (4) solving SVM with quantum annealing (D-Wave); and (5) solving SVM using Graver Augmented Multi-seed Algorithm (GAMA). GAMA is tried with several different numbers of Graver elements and a number of seeds using both simulating annealing and quantum annealing. We found that simulated annealing and GAMA (with simulated annealing) are comparable, provide accurate results quickly, competitive with LibSVM, and superior to Gurobi and quantum annealing.

Decoding of Polar Codes Using Quadratic Unconstrained Binary Optimization

Decoding of Polar Codes Using Quadratic Unconstrained Binary Optimization

Novel real number representations in Ising machines and performance evaluation: Combinatorial random number sum and constant division.

Quantum annealing machines are next-generation computers for solving combinatorial optimization problems. Although physical simulations are one of the most promising applications of quantum annealing machines, a method how to embed the target problem into the machines has not been developed except for certain simple examples. In this study, we focus on a method of representing real numbers using binary variables, or quantum bits. One of the most important problems for conducting physical simulation by quantum annealing machines is how to represent the real number with quantum bits. The variables in physical simulations are often represented by real numbers but real numbers must be represented by a combination of binary variables in quantum annealing, such as quadratic unconstrained binary optimization (QUBO). Conventionally, real numbers have been represented by assigning each digit of their binary number representation to a binary variable. Considering the classical annealing point of view, we noticed that when real numbers are represented in binary numbers, there are numbers that can only be reached by inverting several bits simultaneously under the restriction of not increasing a given Hamiltonian, which makes the optimization very difficult. In this work, we propose three new types of real number representation and compared these representations under the problem of solving linear equations. As a result, we found experimentally that the accuracy of the solution varies significantly depending on how the real numbers are represented. We also found that the most appropriate representation depends on the size and difficulty of the problem to be solved and that these differences show a consistent trend for two annealing solvers. Finally, we explain the reasons for these differences using simple models, the minimum required number of simultaneous bit flips, one-way probabilistic bit-flip energy minimization, and simulation of ideal quantum annealing machine.

Quantum‐Based Combinatorial Optimization for Optimal Sensor Placement in Civil Structures

Over the last decade, concepts such as industry 4.0 and the Internet of Things (IoT) have contributed to the increase in the availability and affordability of sensing technology. In this context, structural health monitoring (SHM) arises as an especially interesting field to integrate and develop these new sensing capabilities, given the criticality of structural application for human life and the elevated costs of manual monitoring. Due to the scale of structural systems, one of the main challenges when designing a modern monitoring system is the optimal sensor placement (OSP) problem. The OSP problem is combinatorial in nature, making its exact solution infeasible in most practical cases, usually requiring the use of metaheuristic optimization techniques. While approaches such as genetic algorithms (GAs) have been able to produce significant results in many practical case studies, their ability to scale up to more complex structures is still an area of open research. This study proposes a novel quantum‐based combinatorial optimization approach to solve the OSP problem approximately, within the context of SHM. For this purpose, a quadratic unconstrained binary optimization (QUBO) model formulation is developed, taking as a starting point of the modal strain energy (MSE) of the structure. The framework is tested using numerical simulations of Warren truss bridges of varying scales. The results obtained using the proposed framework are compared against exhaustive search approaches to verify their performance. More importantly, a detailed discussion of the current limitations of the technology and the future paths of research in the area is presented to the reader.

A Collaborative Neurodynamic Algorithm for Quadratic Unconstrained Binary Optimization

A Collaborative Neurodynamic Algorithm for Quadratic Unconstrained Binary Optimization

Oscillator-based optimization: design, emulation, and implementation

The emergence of digital computers has profoundly reshaped our interactions with technology and the processing of information. Despite excelling in data processing and arithmetics, these computers face limitations in tackling complex nondeterministic-polynomial (NP) problems. In response, researchers have started searching for new computational paradigms that possess the natural tendency of solving these problems. Oscillator-based optimizers are one such paradigm, where the idea is to exploit the parallelism of oscillators networks in order to efficiently solve NP problems. This involves a process of mapping a given optimization task to a quadratic unconstrained binary optimization program and then mapping the resulting program onto an inter-oscillator coupling circuit encoding its coefficients. This paper presents a comprehensive approach to constructing oscillator-based optimizers, offering both the rationale for employing oscillator networks and formulas for linking optimization coefficients to inter-oscillator coupling. Here, we cover most aspects of oscillator-based optimization starting from the design of the network up to its technical implementation. Moreover, we provide a platform-independent wave digital algorithm, which allows for emulating our network’s behavior in a highly parallel fashion.Graphical

New advances for quantum‐inspired optimization

AbstractAdvances in quantum computing with applications in combinatorial optimization have evolved at an increasing rate in recent years. The quadratic unconstrained binary optimization (QUBO) model is at the center of these developments and has become recognized as an effective alternative method for representing a wide variety of combinatorial optimization problems. Additional momentum has resulted from the arrival of quantum computers and their ability to solve the Ising spin glass problem, another form of the QUBO model. This paper highlights advances in solving QUBO models and extensions to more general polynomial unconstrained binary optimization (PUBO) models as important alternatives to traditional approaches. Computational experience is provided that compares the performance of unique quantum‐inspired metaheuristic solvers—the Next Generation Quantum (NGQ) solver for QUBO models and the NGQ‐PUBO solver for PUBO models—with the performance of CPLEX and the Dwave quantum advantage solver. Extensive results, including experiments with a set of large set partitioning problems representing the largest QUBO models reported in the literature to date, along with maximum diversity and max cut problem sets, disclose that our solvers outperform both CPLEX and Dwave by a wide margin in terms of both computational time and solution quality.

Quafu-Qcover: Explore combinatorial optimization problems on cloud-based quantum computers

We introduce Quafu-Qcover, an open-source cloud-based software package developed for solving combinatorial optimization problems using quantum simulators and hardware backends. Quafu-Qcover provides a standardized and comprehensive workflow that utilizes the quantum approximate optimization algorithm (QAOA). It facilitates the automatic conversion of the original problem into a quadratic unconstrained binary optimization (QUBO) model and its corresponding Ising model, which can be subsequently transformed into a weight graph. The core of Qcover relies on a graph decomposition-based classical algorithm, which efficiently derives the optimal parameters for the shallow QAOA circuit. Quafu-Qcover incorporates a dedicated compiler capable of translating QAOA circuits into physical quantum circuits that can be executed on Quafu cloud quantum computers. Compared to a general-purpose compiler, our compiler demonstrates the ability to generate shorter circuit depths, while also exhibiting superior speed performance. Additionally, the Qcover compiler has the capability to dynamically create a library of qubits coupling substructures in real-time, utilizing the most recent calibration data from the superconducting quantum devices. This ensures that computational tasks can be assigned to connected physical qubits with the highest fidelity. The Quafu-Qcover allows us to retrieve quantum computing sampling results using a task ID at any time, enabling asynchronous processing. Moreover, it incorporates modules for results preprocessing and visualization, facilitating an intuitive display of solutions for combinatorial optimization problems. We hope that Quafu-Qcover can serve as an instructive illustration for how to explore application problems on the Quafu cloud quantum computers.

Encoding Molecular Docking for Quantum Computers.

Molecular docking is important in drug discovery but is burdensome for classical computers. Here, we introduce Grid Point Matching (GPM) and Feature Atom Matching (FAM) to accelerate pose sampling in molecular docking by encoding the problem into quadratic unconstrained binary optimization (QUBO) models so that it could be solved by quantum computers like the coherent Ising machine (CIM). As a result, GPM shows a sampling power close to that of Glide SP, a method performing an extensive search. Moreover, it is estimated to be 1000 times faster on the CIM than on classical computers. Our methods could boost virtual drug screening of small molecules and peptides in future.

Exploration of new chemical materials using black-box optimization with the D-wave quantum annealer

In materials informatics, searching for chemical materials with desired properties is challenging due to the vastness of the chemical space. Moreover, the high cost of evaluating properties necessitates a search with a few clues. In practice, there is also a demand for proposing compositions that are easily synthesizable. In the real world, such as in the exploration of chemical materials, it is common to encounter problems targeting black-box objective functions where formalizing the objective function in explicit form is challenging, and the evaluation cost is high. In recent research, a Bayesian optimization method has been proposed to formulate the quadratic unconstrained binary optimization (QUBO) problem as a surrogate model for black-box objective functions with discrete variables. Regarding this method, studies have been conducted using the D-Wave quantum annealer to optimize the acquisition function, which is based on the surrogate model and determines the next exploration point for the black-box objective function. In this paper, we address optimizing a black-box objective function containing discrete variables in the context of actual chemical material exploration. In this optimization problem, we demonstrate results obtaining parameters of the acquisition function by sampling from a probability distribution with variance can explore the solution space more extensively than in the case of no variance. As a result, we found combinations of substituents in compositions with the desired properties, which could only be discovered when we set an appropriate variance.

Individual subject evaluated difficulty of adjustable mazes generated using quantum annealing

In this study, the maze generation using quantum annealing is proposed. We reformulate a standard algorithm to generate a maze into a specific form of a quadratic unconstrained binary optimization problem suitable for the input of the quantum annealer. To generate more difficulty mazes, we introduce an additional cost function Qupdate to increase the difficulty. The difficulty of the mazes was evaluated by the time to solve the maze of 12 human subjects. To check the efficiency of our scheme to create the maze, we investigated the time-to-solution of a quantum processing unit, classical computer, and hybrid solver. The results show that Qupdate generates difficult mazes tailored to the individual. Furthermore, it show that the quantum processing unit is more efficient at generating mazes than other solvers. Finally, we also present applications how our results could be used in the future.

An Efficient Closed-Form Formula for Evaluating r-Flip Moves in Quadratic Unconstrained Binary Optimization

Quadratic unconstrained binary optimization (QUBO) is a classic NP-hard problem with an enormous number of applications. Local search strategy (LSS) is one of the most fundamental algorithmic concepts and has been successfully applied to a wide range of hard combinatorial optimization problems. One LSS that has gained the attention of researchers is the r-flip (also known as r-Opt) strategy. Given a binary solution with n variables, the r-flip strategy “flips” r binary variables to obtain a new solution if the changes improve the objective function. The main purpose of this paper is to develop several results for the implementation of r-flip moves in QUBO, including a necessary and sufficient condition that when a 1-flip search reaches local optimality, the number of candidates for implementation of the r-flip moves can be reduced significantly. The results of the substantial computational experiments are reported to compare an r-flip strategy-embedded algorithm and a multiple start tabu search algorithm on a set of benchmark instances and three very-large-scale QUBO instances. The r-flip strategy implemented within the algorithm makes the algorithm very efficient, leading to very high-quality solutions within a short CPU time.

Cutting-Edge Trajectory Optimization through Quantum Annealing

This paper introduces an innovative approach to explore the capabilities of Quantum Annealing (QA) for trajectory optimization in dynamic systems. The proposed method involves transforming trajectory optimization problems into equivalent binary optimization problems using Quadratic Unconstrained Binary Optimization (QUBO) representation. The procedure is general and adaptable, making it applicable to a wide range of optimal control problems that entail the satisfaction of dynamic, boundary, and path constraints. Specifically, the trajectory is discretized and approximated using polynomials. In contrast to the conventional approach of determining the polynomial degree (n) solely based on the number of boundary conditions, a specific factor is introduced in our method to augment the polynomial degree. As a result, the ultimate polynomial degree is calculated as a composite of two components: n = l + (m−1), where m denotes the count of boundary conditions and l signifies the number of independent variables. By leveraging inverse dynamics, the control required to follow the approximated trajectory can be determined as a linear function of independent variables l. As a result, the optimization function, which is represented by the integral of the square of the control, can be formulated as a QUBO problem and the QA is employed to find the optimal binary solutions.