Hybrid Classical-Quantum Branch-and-Bound Algorithm for Solving Integer Linear Problems

Quantum annealers are suited to solve several logistic optimization problems expressed in the QUBO formulation. However, the solutions proposed by the quantum annealers are generally not optimal, as thermal noise and other disturbing effects arise when the number of qubits involved in the calculation is too large. In order to deal with this issue, we propose the use of the classical branch-and-bound algorithm, that divides the problem into sub-problems which are described by a lower number of qubits. We analyze the performance of this method on two problems, the knapsack problem and the traveling salesman problem. Our results show the advantages of this method, that balances the number of steps that the algorithm has to make with the amount of error in the solution found by the quantum hardware that the user is willing to risk. The results are obtained using the commercially available quantum hardware D-Wave Advantage, and they outline the strategy for a practical application of the quantum annealers.

Haplotype-resolved assembly of diploid and polyploid genomes using quantum computing

Precision medicine's emphasis on individual genetic variants highlights the importance of haplotype-resolved assembly, a computational challenge in bioinformatics given its combinatorial nature. While classical algorithms have made strides in addressing this issue, the potential of quantum computing remains largely untapped. Here, we present the vehicle routing problem (VRP) assembler: an approach that transforms this task into a vehicle routing problem, an optimization formulation solvable on a quantum computer. We demonstrate its potential and feasibility through a proof of concept on short synthetic diploid and triploid genomes using a D-Wave quantum annealer. To tackle larger-scale assembly problems, we integrate the VRP assembler with Google's OR-Tools, achieving a haplotype-resolved local assembly across the human major histocompatibility complex (MHC) region. Our results show encouraging performance compared to Hifiasm with phasing accuracy approaching the theoretical limit, underscoring the promising future of quantum computing in bioinformatics.

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.

Exploring the Limitations of Hybrid Adiabatic Quantum Computing for Emission Tomography Reconstruction

Our study explores the feasibility of quantum computing in emission tomography reconstruction, addressing a noisy ill-conditioned inverse problem. In current clinical practice, this is typically solved by iterative methods minimizing a L2 norm. After reviewing quantum computing principles, we propose the use of a commercially available quantum annealer and employ corresponding hybrid solvers, which combine quantum and classical computing to handle more significant problems. We demonstrate how to frame image reconstruction as a combinatorial optimization problem suited for these quantum annealers and hybrid systems. Using a toy problem, we analyze reconstructions of binary and integer-valued images with respect to their image size and compare them to conventional methods. Additionally, we test our method's performance under noise and data underdetermination. In summary, our method demonstrates competitive performance with traditional algorithms for binary images up to an image size of 32×32 on the toy problem, even under noisy and underdetermined conditions. However, scalability challenges emerge as image size and pixel bit range increase, restricting hybrid quantum computing as a practical tool for emission tomography reconstruction until significant advancements are made to address this issue.

Quantum annealing-based clustering of single cell RNA-seq data

Cluster analysis is a crucial stage in the analysis and interpretation of single-cell gene expression (scRNA-seq) data. It is an inherently ill-posed problem whose solutions depend heavily on hyper-parameter and algorithmic choice. The popular approach of K-means clustering, for example, depends heavily on the choice of K and the convergence of the expectation-maximization algorithm to local minima of the objective. Exhaustive search of the space for multiple good quality solutions is known to be a complex problem. Here, we show that quantum computing offers a solution to exploring the cost function of clustering by quantum annealing, implemented on a quantum computing facility offered by D-Wave [1]. Out formulation extracts minimum vertex cover of an affinity graph to sub-sample the cell population and quantum annealing to optimise the cost function. A distribution of low-energy solutions can thus be extracted, offering alternate hypotheses about how genes group together in their space of expressions.

Quantum advantage in variational Bayes inference

Variational Bayes (VB) inference algorithm is used widely to estimate both the parameters and the unobserved hidden variables in generative statistical models. The algorithm-inspired by variational methods used in computational physics-is iterative and can get easily stuck in local minima, even when classical techniques, such as deterministic annealing (DA), are used. We study a VB inference algorithm based on a nontraditional quantum annealing approach-referred to as quantum annealing variational Bayes (QAVB) inference-and show that there is indeed a quantum advantage to QAVB over its classical counterparts. In particular, we show that such better performance is rooted in key quantum mechanics concepts: i) The ground state of the Hamiltonian of a quantum system-defined from the given data-corresponds to an optimal solution for the minimization problem of the variational free energy at very low temperatures; ii) such a ground state can be achieved by a technique paralleling the quantum annealing process; and iii) starting from this ground state, the optimal solution to the VB problem can be achieved by increasing the heat bath temperature to unity, and thereby avoiding local minima introduced by spontaneous symmetry breaking observed in classical physics based VB algorithms. We also show that the update equations of QAVB can be potentially implemented using ⌈logK⌉ qubits and 𝒪(K) operations per step, where K is the number of values hidden categorical variables can take. Thus, QAVB can match the time complexity of existing VB algorithms, while delivering higher performance.

Quantum Annealing Designs Nonhemolytic Antimicrobial Peptides in a Discrete Latent Space

Increasing the variety of antimicrobial peptides is crucial in meeting the global challenge of multi-drug-resistant bacterial pathogens. While several deep-learning-based peptide design pipelines are reported, they may not be optimal in data efficiency. High efficiency requires a well-compressed latent space, where optimization is likely to fail due to numerous local minima. We present a multi-objective peptide design pipeline based on a discrete latent space and D-Wave quantum annealer with the aim of solving the local minima problem. To achieve multi-objective optimization, multiple peptide properties are encoded into a score using non-dominated sorting. Our pipeline is applied to design therapeutic peptides that are antimicrobial and non-hemolytic at the same time. From 200 000 peptides designed by our pipeline, four peptides proceeded to wet-lab validation. Three of them showed high anti-microbial activity, and two are non-hemolytic. Our results demonstrate how quantum-based optimizers can be taken advantage of in real-world medical studies.

Quantum Annealing Boosts Prediction of Multimolecular Adsorption on Solid Surfaces Avoiding Combinatorial Explosion

Quantum annealing has been used to predict molecular adsorption on solid surfaces. Evaluation of adsorption, which takes place in all solid surface reactions, is a crucially important subject for study in various fields. However, predicting the most stable coordination by theoretical calculations is challenging for multimolecular adsorption because there are numerous candidates. This report presents a novel method for quick adsorption coordination searches using the quantum annealing principle without combinatorial explosion. This method exhibited much faster search and more stable molecular arrangement findings than conventional methods did, particularly in a high coverage region. We were able to complete a configurational prediction of the adsorption of 16 molecules in 2286 s (including 2154 s for preparation, only required once), whereas previously it has taken 38 601 s. This approach accelerates the tuning of adsorption behavior, especially in composite materials and large-scale modeling, which possess more combinations of molecular configurations.

Dynamic Asset Allocation with Expected Shortfall via Quantum Annealing

Recent advances in quantum hardware offer new approaches to solve various optimization problems that can be computationally expensive when classical algorithms are employed. We propose a hybrid quantum-classical algorithm to solve a dynamic asset allocation problem where a target return and a target risk metric (expected shortfall) are specified. We propose an iterative algorithm that treats the target return as a constraint in a Markowitz portfolio optimization model, and dynamically adjusts the target return to satisfy the targeted expected shortfall. The Markowitz optimization is formulated as a Quadratic Unconstrained Binary Optimization (QUBO) problem. The use of the expected shortfall risk metric enables the modeling of extreme market events. We compare the results from D-Wave's 2000Q and Advantage quantum annealers using real-world financial data. Both quantum annealers are able to generate portfolios with more than 80% of the return of the classical optimal solutions, while satisfying the expected shortfall. We observe that experiments on assets with higher correlations tend to perform better, which may help to design practical quantum applications in the near term.

Development of research network on Quantum Annealing Computation and Information using Google Scholar data

We build and analyse the network of 100 top-cited nodes (research papers and books from Google Scholar; the strength or citation of the nodes range from about 44 000 up to 100) starting in early 1980 until last year. These searched publications (papers and books) are based on Quantum Annealing Computation and Information categorized into four different sets: (A) Quantum/Transverse Field Spin Glass Model, (B) Quantum Annealing, (C) Quantum Adiabatic Computation and (D) Quantum Computation Information in the title or abstract of the searched publications. We fitted the growth in the annual number of publication ([Formula: see text]) in each of these four categories, A-D, to the form [Formula: see text] where [Formula: see text] denotes the time in years. We found the scaling time [Formula: see text] to be of the order of about 10 years for categories A and C, whereas [Formula: see text] is of the order of about 5 years for categories B and D. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Role of quantum fluctuation in inducing ergodicity in the spin glass phase and its effect in quantum annealing

We review the numerical studies on the critical behaviour of the quantum Sherrington-Kirkpatrick (SK) spin glass model, which indicate that a quantum critical behaviour is observed up to a low but non-zero value of temperature. We revisit the numerical investigations on the spin glass order parameter distributions, which identify a low temperature along with high transverse field spin glass phase where the order parameter distribution becomes a delta function in the thermodynamic limit indicating the restoration of replica symmetry and ergodic nature of the system. In the remaining spin glass phase associated with high temperature and low transverse field, the observed distribution is broad akin to the Parisi order parameter distribution. This essentially indicates the non-ergodic behaviour of the system. We further discuss the annealing dynamics studies on the quantum SK model. Such investigations reveal the system size independence of annealing time when the annealing paths go through the ergodic spin glass region. Interestingly, when such dynamics are performed in the non-ergodic spin glass phase the annealing time becomes an increasing function of the system size. Spin autocorrelation shows faster relaxation in the ergodic spin glass region compared with that found in the non-ergodic spin glass region. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Statistics of the number of defects after quantum annealing in a thermal environment

We study the statistics of the kink number generated by quantum annealing in a one-dimensional transverse Ising model coupled to a bosonic thermal bath. Using the freezing ansatz for quantum annealing in the thermal environment, we show the relation between the ratio of the second to the first cumulant of the kink number distribution and the average kink density. The theoretical result is confirmed thoroughly by numerical simulation using the non-Markovian infinite time-evolving block decimation which we proposed recently. The simulation using D-Wave's quantum annealer is also discussed. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Greedy parameter optimization for diabatic quantum annealing

A shorter processing time is desirable for quantum computation to minimize the effects of noise. We propose a simple procedure to variationally determine a set of parameters in the transverse-field Ising model for quantum annealing (QA) appended with a field along the [Formula: see text]-axis. The method consists of greedy optimization of the signs of coefficients of the [Formula: see text]-field term based on the outputs of short annealing processes. We test the idea in the ferromagnetic system with all-to-all couplings and spin-glass problems, and find that the method outperforms the traditional form of QA and simulated annealing in terms of the success probability and the time to solution, in particular, in the case of shorter annealing times, achieving the goal of improved performance while avoiding noise. The non-stoquastic [Formula: see text] term can be eliminated by a rotation in the spin space, resulting in a non-trivial diabatic control of the coefficients in the stoquastic transverse-field Ising model, which may be feasible for experimental realization. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Threshold theorem in isolated quantum dynamics with stochastic control errors

We investigate the effect of stochastic control errors in the time-dependent Hamiltonian on isolated quantum dynamics. The control errors are formulated as time-dependent stochastic noise in the Schrödinger equation. For a class of stochastic control errors, we establish a threshold theorem that provides a sufficient condition to obtain the target state, which should be determined in noiseless isolated quantum dynamics, as a relation between the number of measurements and noise strength. The theorem guarantees that if the sum of the noise strengths is less than the inverse of computational time, the target state can be obtained through a constant-order number of measurements. If the opposite is true, the number of measurements to guarantee obtaining the target state increases exponentially with computational time. Our threshold theorem can be applied to any isolated quantum dynamics such as quantum annealing and adiabatic quantum computation. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Quantum adiabatic theorem for unbounded Hamiltonians with a cutoff and its application to superconducting circuits

We present a new quantum adiabatic theorem that allows one to rigorously bound the adiabatic timescale for a variety of systems, including those described by originally unbounded Hamiltonians that are made finite-dimensional by a cutoff. Our bound is geared towards the qubit approximation of superconducting circuits and presents a sufficient condition for remaining within the [Formula: see text]-dimensional qubit subspace of a circuit model of [Formula: see text] qubits. The novelty of this adiabatic theorem is that, unlike previous rigorous results, it does not contain [Formula: see text] as a factor in the adiabatic timescale, and it allows one to obtain an expression for the adiabatic timescale independent of the cutoff of the infinite-dimensional Hilbert space of the circuit Hamiltonian. As an application, we present an explicit dependence of this timescale on circuit parameters for a superconducting flux qubit and demonstrate that leakage out of the qubit subspace is inevitable as the tunnelling barrier is raised towards the end of a quantum anneal. We also discuss a method of obtaining a [Formula: see text] effective Hamiltonian that best approximates the true dynamics induced by slowly changing circuit control parameters. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Understanding domain-wall encoding theoretically and experimentally

We analyse the method of encoding pairwise interactions of higher-than-binary discrete variables (these models are sometimes referred to as discrete quadratic models) into binary variables based on domain walls on one-dimensional Ising chains. We discuss how this is relevant to quantum annealing, but also many gate model algorithms such as VQE and QAOA. We theoretically show that for problems of practical interest for quantum computing and assuming only quadratic interactions are available between the binary variables, it is not possible to have a more efficient general encoding in terms of number of binary variables per discrete variable. We furthermore use a D-Wave Advantage 1.1 flux qubit quantum annealing computer to show that the dynamics effectively freeze later for a domain-wall encoding compared with a traditional one-hot encoding. This second result could help explain the dramatic performance improvement of domain wall over one-hot, which has been seen in a recent experiment on D-Wave hardware. This is an important result because usually problem encoding and the underlying physics are considered separately, our work suggests that considering them together may be a more useful paradigm. We argue that this experimental result is also likely to carry over to a number of other settings, we discuss how this has implications for gate-model and quantum-inspired algorithms. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Quantum Annealing in the NISQ Era: Railway Conflict Management

We are in the noisy intermediate-scale quantum (NISQ) devices' era, in which quantum hardware has become available for application in real-world problems. However, demonstrations of the usefulness of such NISQ devices are still rare. In this work, we consider a practical railway dispatching problem: delay and conflict management on single-track railway lines. We examine the train dispatching consequences of the arrival of an already delayed train to a given network segment. This problem is computationally hard and needs to be solved almost in real time. We introduce a quadratic unconstrained binary optimization (QUBO) model of this problem, which is compatible with the emerging quantum annealing technology. The model's instances can be executed on present-day quantum annealers. As a proof-of-concept, we solve selected real-life problems from the Polish railway network using D-Wave quantum annealers. As a reference, we also provide solutions calculated with classical methods, including the conventional solution of a linear integer version of the model as well as the solution of the QUBO model using a tensor network-based algorithm. Our preliminary results illustrate the degree of difficulty of real-life railway instances for the current quantum annealing technology. Moreover, our analysis shows that the new generation of quantum annealers (the advantage system) does not perform well on those instances, either.

Quantum annealing for industry applications: introduction and review

Quantum annealing (QA) is a heuristic quantum optimization algorithm that can be used to solve combinatorial optimization problems. In recent years, advances in quantum technologies have enabled the development of small- and intermediate-scale quantum processors that implement the QA algorithm for programmable use. Specifically, QA processors produced by D-Wave systems have been studied and tested extensively in both research and industrial settings across different disciplines. In this paper we provide a literature review of the theoretical motivations for QA as a heuristic quantum optimization algorithm, the software and hardware that is required to use such quantum processors, and the state-of-the-art applications and proofs-of-concepts that have been demonstrated using them. The goal of our review is to provide a centralized and condensed source regarding applications of QA technology. We identify the advantages, limitations, and potential of QA for both researchers and practitioners from various fields.

Feature Selection for Recommender Systems with Quantum Computing

The promise of quantum computing to open new unexplored possibilities in several scientific fields has been long discussed, but until recently the lack of a functional quantum computer has confined this discussion mostly to theoretical algorithmic papers. It was only in the last few years that small but functional quantum computers have become available to the broader research community. One paradigm in particular, quantum annealing, can be used to sample optimal solutions for a number of NP-hard optimization problems represented with classical operations research tools, providing an easy access to the potential of this emerging technology. One of the tasks that most naturally fits in this mathematical formulation is feature selection. In this paper, we investigate how to design a hybrid feature selection algorithm for recommender systems that leverages the domain knowledge and behavior hidden in the user interactions data. We represent the feature selection as an optimization problem and solve it on a real quantum computer, provided by D-Wave. The results indicate that the proposed approach is effective in selecting a limited set of important features and that quantum computers are becoming powerful enough to enter the wider realm of applied science.

Dynamic Topology Reconfiguration of Boltzmann Machines on Quantum Annealers

Boltzmann machines have useful roles in deep learning applications, such as generative data modeling, initializing weights for other types of networks, or extracting efficient representations from high-dimensional data. Most Boltzmann machines use restricted topologies that exclude looping connectivity, as such connectivity creates complex distributions that are difficult to sample. We have used an open-system quantum annealer to sample from complex distributions and implement Boltzmann machines with looping connectivity. Further, we have created policies mapping Boltzmann machine variables to the quantum bits of an annealer. These policies, based on correlation and entropy metrics, dynamically reconfigure the topology of Boltzmann machines during training and improve performance.

A QUBO Formulation of the Stereo Matching Problem for D-Wave Quantum Annealers

In this paper, we propose a methodology to solve the stereo matching problem through quantum annealing optimization. Our proposal takes advantage of the existing Min-Cut/Max-Flow network formulation of computer vision problems. Based on this network formulation, we construct a quadratic pseudo-Boolean function and then optimize it through the use of the D-Wave quantum annealing technology. Experimental validation using two kinds of stereo pair of images, random dot stereograms and gray-scale, shows that our methodology is effective.

Finding a Hadamard Matrix by Simulated Quantum Annealing

Hard problems have recently become an important issue in computing. Various methods, including a heuristic approach that is inspired by physical phenomena, are being explored. In this paper, we propose the use of simulated quantum annealing (SQA) to find a Hadamard matrix, which is itself a hard problem. We reformulate the problem as an energy minimization of spin vectors connected by a complete graph. The computation is conducted based on a path-integral Monte-Carlo (PIMC) SQA of the spin vector system, with an applied transverse magnetic field whose strength is decreased over time. In the numerical experiments, the proposed method is employed to find low-order Hadamard matrices, including the ones that cannot be constructed trivially by the Sylvester method. The scaling property of the method and the measurement of residual energy after a sufficiently large number of iterations show that SQA outperforms simulated annealing (SA) in solving this hard problem.

Efficiency of quantum vs. classical annealing in nonconvex learning problems

Quantum annealers aim at solving nonconvex optimization problems by exploiting cooperative tunneling effects to escape local minima. The underlying idea consists of designing a classical energy function whose ground states are the sought optimal solutions of the original optimization problem and add a controllable quantum transverse field to generate tunneling processes. A key challenge is to identify classes of nonconvex optimization problems for which quantum annealing remains efficient while thermal annealing fails. We show that this happens for a wide class of problems which are central to machine learning. Their energy landscapes are dominated by local minima that cause exponential slowdown of classical thermal annealers while simulated quantum annealing converges efficiently to rare dense regions of optimal solutions.

Ergodic and localized regions in quantum spin glasses on the Bethe lattice

By considering the quantum dynamics of a transverse-field Ising spin glass on the Bethe lattice, we find the existence of a many-body localized (MBL) region at small transverse field and low temperature. The region is located within the thermodynamic spin glass phase. Accordingly, we conjecture that quantum dynamics inside the glassy region is split into a small MBL region and a large delocalized (but not necessarily ergodic) region. This has implications for the analysis of the performance of quantum adiabatic algorithms.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.

Social Milieu Oriented Routing: A New Dimension to Enhance Network Security in WSNs

In large-scale wireless sensor networks (WSNs), in order to enhance network security, it is crucial for a trustor node to perform social milieu oriented routing to a target a trustee node to carry out trust evaluation. This challenging social milieu oriented routing with more than one end-to-end Quality of Trust (QoT) constraint has proved to be NP-complete. Heuristic algorithms with polynomial and pseudo-polynomial-time complexities are often used to deal with this challenging problem. However, existing solutions cannot guarantee the efficiency of searching; that is, they can hardly avoid obtaining partial optimal solutions during a searching process. Quantum annealing (QA) uses delocalization and tunneling to avoid falling into local minima without sacrificing execution time. This has been proven a promising way to many optimization problems in recently published literatures. In this paper, for the first time, with the help of a novel approach, that is, configuration path-integral Monte Carlo (CPIMC) simulations, a QA-based optimal social trust path (QA_OSTP) selection algorithm is applied to the extraction of the optimal social trust path in large-scale WSNs. Extensive experiments have been conducted, and the experiment results demonstrate that QA_OSTP outperforms its heuristic opponents.

Using thermal boundary conditions to engineer the quantum state of a bulk magnet

The degree of contact between a system and the external environment can alter dramatically its proclivity to quantum mechanical modes of relaxation. We show that controlling the thermal coupling of cubic-centimeter-sized crystals of the Ising magnet LiHo(x)Y(1-x)F4 to a heat bath can be used to tune the system between a glassy state dominated by thermal excitations over energy barriers and a state with the hallmarks of a quantum spin liquid. Application of a magnetic field transverse to the Ising axis introduces both random magnetic fields and quantum fluctuations, which can retard and speed the annealing process, respectively, thereby providing a mechanism for continuous tuning between the destination states. The nonlinear response of the system explicitly demonstrates quantum interference between internal and external relaxation pathways.