<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <issn pub-type="ppub">1613-0073</issn>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Does the structure of the QUBO problem afect the efectiveness of quantum annealing? An empirical perspective</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Politecnico di Milano</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Italy</string-name>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Riccardo Pellini</string-name>
          <email>riccardo.pellini@polimi.it</email>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Maurizio Ferrari Dacrema</string-name>
          <email>maurizio.ferrari@polimi.it</email>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Paolo Cremonesi</string-name>
          <email>paolo.cremonesi@polimi.it</email>
        </contrib>
      </contrib-group>
      <pub-date>
        <year>2023</year>
      </pub-date>
      <abstract>
        <p>In recent years there has been a significant interest in exploring the potential of Quantum Annealers (QA) as heuristic solvers of Quadratic Unconstrained Binary Optimization (QUBO) problems. Some problems are more dificult to solve on QA and understanding why is not straightforward, because an analytical study of the underlying physical system is intractable for large QUBO problems. This work consists in an empirical analysis of the features making a QUBO problem dificult to solve on QA, based on clusters of QUBO instances identified with Hierarchical Clustering. The analysis reveals correlations between specific values of the features and the ability of QA to solve efectively the instances. These initial results open new research opportunities to inform the development of new AI methods supporting quantum computation (e.g., for minor embedding or error mitigation) that are better tailored to the characteristics of the problem, as well as to develop better QUBO formulations for known problems in order to improve the quality of the solutions found by QA. 1 htp:/ceur-ws.org CEUR Workshop Proceedings (CEUR-WS.org) ISN1613-073</p>
      </abstract>
      <kwd-group>
        <kwd>quantum computing</kwd>
        <kwd>quantum annealing</kwd>
        <kwd>optimization</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>CEUR
ceur-ws.org</p>
    </sec>
    <sec id="sec-2">
      <title>1. Introduction</title>
      <p>Quantum Annealers (QA) are heuristic solvers of Quadratic Unconstrained Binary Optimization
(QUBO) problems. It is known that some problems are more dificult to solve efectively with
QA compared to other ones. However, it’s challenging to study analytically the Hamiltonian
for large QUBO problems and it is often not clear how to use these findings to develop new
general QUBO formulations that are easier to solve on QA. Furthermore, AI techniques based
on the characteristics of QUBO problems able to support QA, such as for minor embedding or
error mitigation, are lacking, since the characteristics which represent the dificulty of a QUBO
problem are unknown.</p>
      <p>
        In this work, we study with an empirical perspective the characteristics making a QUBO
problem dificult to solve on QA, in particular when it requires too many qubits for the analytical
study of the Hamiltonian, trying to answer to the question: does the structure of a QUBO problem
afects the efectiveness of QA?
1These results have been presented at the 12th Adiabatic Quantum Computing Conference (AQC 2023), June
CEUR
Workshop
Proceedings
In this study we consider instances of the Maximum Cut, Minimum Vertex Cover, Graph
Coloring, Set Partitioning, Number Partitioning problems [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. Each instance has 30-32 QUBO
variables, corresponding to 100-150 qubits on the physical embedding of QA. For all instances
we compute several features based on the energy distribution of their solutions and on the
spectral representation of the QUBO problem as a graph. In particular, the Spectral Flatness
measures the uniformity of the spectrum of the distribution of the solutions of a QUBO instance
[
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]; the Graph Spectral Flatness measures the uniformity of the spectrum of a signal on a
graph [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ], when both the signal and the graph are related to the same QUBO instance.
      </p>
      <p>
        All instances are solved with the D-Wave Advantage QA and also with Simulated Annealing
and Tabu Search, considered as baseline methods. Instances are then clustered on the base of
their features and the clusters are validated with Silhouette Coeficient. In order to corroborate
the obtained results, we also include four test instances related to the Feature Selection problem
[
        <xref ref-type="bibr" rid="ref4">4</xref>
        ] and check to which cluster they are closest to, according to the computed features.
      </p>
    </sec>
    <sec id="sec-3">
      <title>3. Results and Discussion</title>
      <p>The analysis reveals correlations between the clusters and the ability of QA to solve the instances
efectively. In particular, we see that instances characterized by high levels of Spectral Flatness
are solved optimally by QA. Furthermore, we see also that instances which have low levels of
Graph Spectral Flatness are solved optimally too by the QA.</p>
      <p>All these findings are corroborated by the test instances we have considered, confirming the
relationship between the quality of the solution found by the QA and the structure of the QUBO
problem. These initial results open new research opportunities to inform the development of
new AI methods supporting quantum computation that are better tailored to the characteristics
of the problem, as well as to develop better QUBO formulations for known problems in order to
improve the quality of the solutions found by QA.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <surname>Glover</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kochenberger</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Du</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          ,
          <article-title>Quantum bridge analytics i: a tutorial on formulating and using qubo models</article-title>
          ., 4OR
          <volume>17</volume>
          (
          <issue>4</issue>
          ) (
          <year>2019</year>
          )
          <fpage>335</fpage>
          -
          <lpage>371</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>B.</given-names>
            <surname>Boashash</surname>
          </string-name>
          , G. Azemi,
          <string-name>
            <given-names>N. Ali</given-names>
            <surname>Khan</surname>
          </string-name>
          ,
          <article-title>Principles of time-frequency feature extraction for change detection in non-stationary signals: Applications to newborn eeg abnormality detection</article-title>
          ,
          <source>Pattern Recognition</source>
          <volume>48</volume>
          (
          <year>2015</year>
          )
          <fpage>616</fpage>
          -
          <lpage>627</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>A.</given-names>
            <surname>Ortega</surname>
          </string-name>
          , Introduction to Graph Signal Processing, Cambridge University Press,
          <year>2022</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>M.</given-names>
            <surname>Ferrari Dacrema</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Moroni</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Nembrini</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Ferro</surname>
          </string-name>
          , G. Faggioli,
          <string-name>
            <given-names>P.</given-names>
            <surname>Cremonesi</surname>
          </string-name>
          ,
          <article-title>Towards feature selection for ranking and classification exploiting quantum annealers</article-title>
          ,
          <source>in: SIGIR '22: The 45th International ACM SIGIR Conference on Research and Development in Information Retrieval</source>
          , Madrid, Spain,
          <source>July 11 - 15</source>
          ,
          <year>2022</year>
          , ACM,
          <year>2022</year>
          , pp.
          <fpage>2814</fpage>
          -
          <lpage>2824</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>