JP2018529142A - 量子処理装置及び方法 - Google Patents

Abstract

複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法を提供する。前記方法は、計算問題を量子システムの問題ハミルトニアンにコード化するステップを含む。前記問題ハミルトニアンは、複数の調整可能なパラメータを有する単体ハミルトニアンであり、前記コード化するステップは、前記計算問題から前記複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む。前記方法は、前記量子システムを初期量子状態から前記量子システムの最終ハミルトニアンの基底状態に向けて移行させるステップを更に含む。前記最終ハミルトニアンは前記問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンの複数の調整可能なパラメータは前記問題コード化設定に存在し、前記短距離ハミルトニアンはｄが前記計算問題から独立しているｄ体ハミルトニアンである。前記方法は、前記複数のキュービットの少なくとも一部分を測定することで、前記量子システムの読み出しを取得するステップを更に含む。前記方法は、前記読み出しから前記計算問題の解を決定するステップを更に含む。
【選択図】図１

Description

本発明は、量子システム、具体的には、複数の量子ビット（キュービット）を有する量子システムを用いて、計算問題の解を計算する装置及び方法に関する。

古典的な情報処理に基づく計算装置（つまり、量子力学的効果を利用しない計算装置）は、かつて、特定の計算しか行うことができないハードウェアに組み込まれた計算機として始まった。完全にプログラム可能なコンピュータへの移行が当該技術分野に革命を起こし、情報化時代が始まった。現在、量子計算装置（つまり、計算問題を解くために、場合によっては古典的な情報処理を利用することに加え、量子力学的効果を利用する計算装置）は、ある意味、ハードウェアに組み込まれた計算機の段階にあり、特別に設計された、すなわちハードウェアに組み込まれた計算問題に対してのみ対処可能である。具体的には、全てのプラットフォーム及び規律（disciplines）にまたがって現存する全ての量子計算装置は、完全にプログラム可能で拡張性のあるものにまで未だ達していない。

例えば、Ｄ−ＷＡＶＥ ＳＹＳＴＥＭＳ社製の超電導キュービットに基づく量子計算装置は、キュービットの２つのグループを形成する。この２つのグループにおいて、異なるグループのキュービット間では任意の相互作用が生じるが、同じグループのキュービット間では相互作用は生じない。新たなキュービットのグループを追加しても良いが、この場合には同じ制約を受ける。したがって、この量子計算装置はハードウェアに組み込まれた制約により制限的であるので、必要な相互作用に対し、完全にプログラム可能で且つ拡張性のあるアーキテクチャを実現することが不可能である。

したがって、量子システムを用いて計算問題を解くためには、改善された方法及び装置が必要とされる。

実施形態によれば、複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法が提供される。前記方法は、前記計算問題を量子システムの問題ハミルトニアンにコード化するステップを含む。問題ハミルトニアンは複数の調整可能なパラメータを有する単体ハミルトニアンであり、前記コード化するステップは、前記複数の調整可能なパラメータに対する計算問題から問題コード化設定を決定するステップを含む。前記方法は、前記量子システムを初期量子状態から前記量子システムの最終ハミルトニアンの基底状態に向けて移行させるステップを更に含む。前記最終ハミルトニアンは前記問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンの前記複数の調整可能なパラメータは前記問題コード化設定に存在し、前記短距離ハミルトニアンはｄが前記計算問題から独立しているｄ体ハミルトニアンである。前記方法は、前記複数のキュービットの少なくとも一部分を測定することで、前記量子システムの読み出しを取得するステップを更に含む。前記方法は、前記読み出しから前記計算問題の解を決定するステップを更に含む。

更なる実施形態によれば、複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法が提供される。前記方法は、前記計算問題を量子システムの問題ハミルトニアンにコード化するステップを含む。前記問題ハミルトニアンは複数の調整可能なパラメータを有する単体ハミルトニアンであり、前記コード化するステップは、前記計算問題から前記複数の調整可能なパラメータに対し問題コード化設定を決定するステップを含む。前記方法は、初期量子状態の前記量子システムを初期化するステップを更に含む。前記方法は、量子アニーリングを実行することにより前記量子システムを前記初期量子状態から最終量子状態に移行させるステップを更に含む。量子アニーリングを実行するステップは、前記量子システムの初期ハミルトニアンから前記量子システムの最終ハミルトニアンに移行するステップを含む。前記最終ハミルトニアンは前記問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンの前記複数の調整可能なパラメータは前記問題コード化設定に存在し、前記短距離ハミルトニアンはｄ体ハミルトニアンであり、このｄは前記計算問題から独立している。前記方法は、前記複数のキュービットの少なくとも一部分を測定することで前記最終量子状態の読み出しを取得するステップを更に含む。前記方法は、前記読み出しから前記計算問題の解を決定するステップを更に含む。

更なる実施形態によれば、計算問題の解を計算する装置が提供される。前記装置は、複数のキュービットを有する量子システムを備える。前記装置は、前記量子システムを前記量子システムの基底状態に向けて冷却するのに適した冷却ユニットを更に備える。前記装置は、量子アニーリングによって前記量子システムの初期ハミルトニアンを前記量子システムの最終ハミルトニアンに移行させるのに適したプログラム可能な量子アニーリングユニットを更に備える。前記最終ハミルトニアンは問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンは複数の調整可能なパラメータを有する単体ハミルトニアンである。前記装置は、前記複数のキュービットの少なくとも一部分を測定するのに適した測定装置を更に備える。前記装置は、前記プログラム可能な量子アニーリングユニット及び前記測定装置に接続された古典的計算システムを更に備える。前記古典的計算システムは、計算問題を入力として受信するステップと、前記計算問題を前記問題ハミルトニアンにコード化するステップであって、前記計算問題から前記問題ハミルトニアンの前記複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む、前記コード化するステップと、前記問題コード化設定を前記量子アニーリングユニットに伝達するステップと、を実行するように構成される。前記プログラム可能な量子アニーリングユニットは、前記古典的計算システムから前記問題コード化設定を受信するステップと、量子アニーリングによって前記初期ハミルトニアンを前記最終ハミルトニアンに移行させるステップであって、前記問題ハミルトニアンの前記複数の調整可能なパラメータが前記問題コード化設定に存在する、前記移行させるステップと、を実行するように構成される。前記古典的計算システムは、前記測定装置から前記量子システムの読み出しを受信するステップと、前記読み出しから前記計算問題の解を決定するステップと、を更に実行するように構成される。

更なる実施形態によれば、計算問題の解を計算するためのプログラム可能な量子アニーリング装置が提供される。前記プログラム可能な量子アニーリング装置は、２次元格子に沿って配置された複数の超電導キュービットを有する量子システムを備える。前記プログラム可能な量子アニーリング装置は、複数の調整可能な磁束を生成するように構成された複数の磁束バイアスユニットを有する磁束バイアスアセンブリであって、各調整可能な磁束が前記複数の超電導キュービットのうちの単一の超伝導キュービットに作用する、磁束バイアスアセンブリを更に備える。前記プログラム可能な量子アニーリング装置は、プラケットハミルトニアンに従い、前記複数の超電導キュービットを結合するように構成された少なくとも１つの超電導量子干渉装置を有する結合ユニットを更に備える。前記プログラム可能な量子アニーリング装置は、磁束バイアスユニット及び結合ユニットに接続された制御装置を更に備える。前記制御装置は、前記量子システムの問題ハミルトニアンの複数の調整可能なパラメータに対する問題コード化設定を受信するステップであって、前記問題ハミルトニアンが単体ハミルトニアンであり、前記問題コード化設定が計算問題をコード化する、前記受信するステップを実行するように構成されている。前記制御装置は、前記磁束バイアスアセンブリ及び前記結合ユニットを制御して、量子アニーリングによって前記量子システムの初期ハミルトニアンを前記量子システムの最終ハミルトニアンに移行させるステップを更に実行するように構成されている。前記最終ハミルトニアンは前記プラケットハミルトニアン及び前記問題ハミルトニアンの総和であり、前記問題ハミルトニアンの前記複数の調整可能なパラメータは前記問題コード化設定に存在する。

実施形態は、前記開示されたシステム及び装置の操作方法や、実施形態に係る方法を実行するための前記開示されたシステムの使用も対象としている。

実施形態に組み合わせ可能な更なる利点、特徴、側面及び詳細については、従属する請求項、明細書及び図面から明白である。

実施形態に係る、量子システムを用いて計算問題の解を計算する装置を示す。 実施形態に係る複数のキュービットの配置例を示す。 実施形態に係る複数のキュービットの配置例を示す。 実施形態に係る複数のキュービットの配置例を示す。 実施形態に係る単体ハミルトニアンの概念を示す。 実施形態に係る短距離ハミルトニアンの概念を示す。 実施形態に係る短距離ハミルトニアンの概念を示す。 複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 実施形態に従った、問題ハミルトニアン及び対応する最終ハミルトニアンへの計算問題の具体的なコード化を示す。 誤差に対する量子システムのロバスト性に関する実施形態の利点を示す。

以下、様々な例の実施形態を詳しく参照する。そのうちの１以上の例については、各図に図示している。各例は、説明のために提示しているのであって、本発明はその例に限定されるものではない。例えば、１つの実施形態の一部として図示されている又は記載されている特徴は、他の実施形態にも使用できる、又は他の実施形態と共に使用することで新たな実施形態を作り出すことが可能である。本明細書の開示は、そのような改良やバリエーションを含むことを意図している。

下記の図面に関する記載において、同一の参照符号は同一の構成要素を意味する。概して、各実施形態については相違点のみを述べる。図面に示す構造は、必ずしも実際の縮尺通りに描写しているわけではなく、実施形態を理解し易くするために誇張された描写を含む場合がある。

実施形態は、複数のキュービットを有する量子システムに関する。本明細書に記載のキュービットは、量子力学２準位システム（quantum mechanical two-level system）とも称される。キュービットは、起こり得るキュービットの量子状態を表した２つの量子基本状態（quantum basis states）｜０＞及び｜１＞を有する。量子力学の重ね合わせの原理によれば、ａ｜０＞＋ｂ｜１＞の式である全ての重ね合わせがキュービットに起こり得る量子状態である。上記のａ及びｂは複素数である。数学的には、キュービットは２次元のベクトル空間で表すことができる。複数のキュービットは、複数のキュービットのそれぞれが量子状態｜０＞又は量子状態｜１＞の何れかである形態に対応する量子基本状態を有する。例えば、キュービットが５つあるとすれば、５つのキュービットの量子基本状態は｜００１０１＞となり得る。量子状態｜００１０１＞は、第１、第２及び第４キュービットが量子状態｜０＞であり、第３及び第５キュービットが量子状態｜１＞である形態を表している。キュービットが複数のｍ個ある場合、２^ｍ通りの量子基本状態が存在する。重ね合わせの原理からして、複数のキュービットに対し２通りの量子状態を考えると、量子基本状態の重ね合わせは複数のキュービットの量子状態でもある。例えば、ａ、ｂ及びｃが複素数である式ａ｜００１０１＞＋ｂ｜１１１１０＞＋ｃ｜１１１１１＞の重ね合わせは、複数のキュービットの量子状態である。数学的には、複数のｍ個のキュービットからなる量子システムは、２^ｍ次元のベクトル空間で表すことができる。

複数のキュービットは、複数の超電導キュービット（例えば、トランスモン（transmon）キュービット又は磁束キュービット）を含んでいてもよいし、又は複数の超電導キュービットから構成されていてもよい。超電導キュービットは第１次及び第２次超電導ループを含んでいてもよい。第１次超電導ループで時計回り及び反時計回りにそれぞれ伝搬する超電導電流は、超電導キュービットの量子基本状態｜１＞及び｜０＞を形成可能である。さらに、第２次超電導ループを通る磁束バイアスは、量子基本状態｜０＞及び｜１＞を結合可能である。

別の態様として、量子システムは、トラップイオンのシステムを利用して実現してもよい。この場合、キュービットの量子基本状態｜０＞及び｜１＞は、ゼーマン若しくは超微細マニホールドの２段階で形成されるか、又はＣａ４０＋等のアルカリ土類イオン若しくはアルカリ土類のように正電荷を帯びたイオンの禁制光学遷移（forbidden optical transition）のいたるところで形成される。

更に別の態様として、量子システムは、レーザフィールドから光格子又は間隔の広い格子に捕獲された極低温の中性アルカリ原子等の極低温原子を利用して実現してもよい。上記の原子は、レーザ冷却により基底状態に向けて進行可能である。キュービットの量子基本状態は、原子の基底状態及び高位リュードベリ状態（high-lying Rydberg state）によって形成される。キュービットはレーザ光によって扱うことができる。

更に別の態様として、量子システムは、量子ドットを用いて実現してもよい。量子ドットキュービットは、ＧａＡｓ／ＡｌＧａＡｓヘテロ構造で組み立てられてもよい。キュービットは、単一井戸型から二重井戸型へポテンシャルを断熱的に調整することによって準備可能なスピン状態でコード化される。

更に別の態様として、量子システムは、ダイヤモンド結晶の点欠陥であるＮＶセンター等の固体結晶の不純物を用いて実現してもよい。他の不純物、例えば、クロム不純物に関係する着色中心、固体結晶の希土類イオン、又は炭化ケイ素の欠陥中心については研究中である。ＮＶセンターは、２つの不対電子を有し、場合によっては周囲の核スピンと共に、キュービットを実現するために使用可能な長寿命の２つの著しい欠陥レベルの特定を可能にするスピン−１基底状態（spin-1 ground state）を与える。

実施形態によれば、量子システムは、１つ以上又は複数の個別のｑ段階量子システムを有していてもよい。ｑは定数である。例えば、ｑは２〜８の範囲内（例えば、３、４、５又は６）である。個別のｑ段階量子システムは、ｑ状態｜０＞，｜１＞，・・・｜ｑ−１＞からなる基底を有する。個別の３段階量子システムは、「キュートリット」と称される。

量子システムのハミルトニアンは、量子システムの単一の相互作用又は複数の相互作用を表すことができる。ハミルトニアンは、量子システムに作用する演算子である。ハミルトニアンの固有値は、量子システムのエネルギースペクトルに対応する。ハミルトニアンの基底状態は、最少エネルギーを有する量子システムの量子状態である。ハミルトニアンの基底状態は、零度の量子状態であってもよい。

本明細書に記載の古典的計算システムは、古典的ビットで作動する計算システムを意味する。古典的計算システムは、古典的ビットで情報を処理するための中央処理装置（ＣＰＵ）、及び／又は古典的ビットで情報を記憶する記憶装置を有していてもよい。古典的計算システムは、パーソナルコンピュータ（ＰＣ）等の１つ以上の従来型コンピュータ、及び／又は従来型コンピュータのネットワークを有していてもよい。

実施形態の詳細を述べる前に、本開示事項のいくつかの側面について、図１を参照して説明する。図１は、実施形態に係る計算問題の解を計算する装置４００の一例を示している。

図１に示す装置４００は、量子システム４２０を利用して計算問題の解を計算するのに適している。量子システム４２０は、図１において黒色点で表された複数のキュービット１００を有する。図１に示す実施形態によれば、複数のキュービット１００は２次元格子１２０、具体的には２次元正方格子に沿って配置されている。

図１は、量子システム４２０を冷却するように構成された冷却ユニット４１０を更に示す。冷却ユニット４１０は、量子システム４２０を動作温度まで冷却することができる。

図１は、古典的計算システム４５０を更に示す。古典的計算システム４５０は、解くべき計算問題４５２を入力として受信するように構成されている。計算問題４５２は、巡回セールスマン問題やイジングスピンモデル問題等のＮＰ困難問題であってもよい。ここでのＮＰは、「非決定性多項式時間（nondeterministic polynomial time）」の略称である。

古典的計算システム４５０は、更に、計算問題４５２を量子システム４２０の問題ハミルトニアン（problem Hamiltonian）４７２にコード化するように構成されている。図１に例示した実施形態によれば、問題ハミルトニアン４７２は、式Ｈ^ｐｒｏｂ＝Σ_ｋＪ_ｋσ_ｚ ^（ｋ）を有する。σ_ｚ ^（ｋ）は、複数のキュービット１００のうちｋ番目のキュービットに作用するパウリ演算子である。各Ｊ_ｋは、１つ以上の外部実体（external entities）によって決定される調整可能なパラメータであり、キュービットｋ毎に個別に調整可能な磁場を例示できる。例えば、Ｊ_ｋは、ｋ番目のキュービットに影響する調整可能な磁場の強度であってもよい。磁場等の複数の調整可能な外部実体が与えられてもよい。この場合、各調整可能な外部実体が複数のキュービットのうちの１つのキュービットに影響する。外部実体を調整することにより、計算問題４５２に応じてパラメータＪ_ｋを調整可能である。

古典的計算システム４５０によって実行される問題ハミルトニアン４７２での計算問題４５２をコード化するステップは、計算問題４５２から、複数の調整可能なパラメータＪ_ｋに対する問題コード化の設定（problem-encoding configuration）を決定するステップを含む。各調整可能なパラメータＪ_ｋに対して、パラメータの値は計算問題４５２に応じて決定される。したがって、問題コード化設定は、計算問題に依存する。

図１は、プログラム可能な量子アニーリングユニット４３０を更に示す。量子アニーリングユニット４３０は、量子システム４２０の初期ハミルトニアンから量子システム４２０の最終ハミルトニアンに移行することで量子アニーリングを実行するのに適している。

図１に関して述べた実施形態によれば、初期ハミルトニアンは式Ｈ^ｉｎｉｔ＝Σ_ｋａ_ｋσ_ｘ ^（ｋ）を有する。ａ_ｋは係数であり、σ_ｘ ^（ｋ）は複数のキュービット１００のｋ番目のキュービットに作用するパウリ演算子である。パウリ演算子σ_ｚ ^（ｋ）及びσ_ｘ ^（ｋ）は非可換のパウリ演算子、特に反可換のパウリ演算子でもよい。初期ハミルトニアンＨ^ｉｎｉｔは計算問題４５２から独立していてもよい。

最終ハミルトニアンは、問題ハミルトニアン４７２及び短距離ハミルトニアン（short-range Hamiltonian）４７４の総和である。例示した実施形態によれば、短距離ハミルトニアン４７４は、プラケットに対応する複数のキュービットからなるグループ間の相互作用を表すプラケットハミルトニアン（plaquette Hamiltonian）である。プラケットは、例えば、キュービットが配置される２次元正方格子の複数の単純な正方形であってもよい。後で詳しく説明する図７は、実施形態に係る２次元格子のプラケット３７０の例を示す。短距離ハミルトニアンは、例えば、４体プラケットハミルトニアン等のｄ体ハミルトニアンであり、ｄは計算問題４５２から独立している。実施形態によれば、短距離ハミルトニアン４７４は、計算問題４５２から独立していてもよい。

上記のように、計算問題４５２は、問題ハミルトニアン４７２でコード化される。特に、調整可能なパラメータＪ_ｋの問題コード化設定でコード化される。実施形態によれば、問題ハミルトニアン４７２及び短距離ハミルトニアン４７４の総和である最終ハミルトニアン４７０が、計算問題４５２の解に関する情報を含む基底状態を有するようにコード化が行われる。したがって、量子システム４２０が最終ハミルトニアン４７０の基底状態にある場合、量子システム４２０を測定することにより、計算問題に関する情報が明らかになり得る。

実施形態によれば、図１の矢符４９９で示すように、量子システム４２０は、最終ハミルトニアン４７０の基底状態に向けて進行する。このとき、問題ハミルトニアン４７２の複数の調整可能なパラメータは、問題コード化設定に存在する。図１に示す実施形態によれば、量子アニーリングユニット４３０が量子アニーリングを実行することにより、量子システム４２０は、最終ハミルトニアン４７０の基底状態に向けて進行する。量子アニーリングの実行は、初期ハミルトニアンから最終ハミルトニアン４７０に移行するステップを含む。

例示した実施形態によれば、量子システム４２０は、量子システム４２０を初期ハミルトニアンの基底状態に向けて冷却することにより、初期量子状態で初期化される。さらに、プログラム可能な量子アニーリングユニット４３０は、量子アニーリングを実行することにより、量子システム４２０を初期時間での初期量子状態から最終時間での最終量子状態に進行させるのに適している。量子アニーリングは、量子システム４２０を初期量子状態から最終量子状態に進行させるために、初期時間での初期ハミルトニアンから最終時間での最終ハミルトニアン４７０に移行するステップを含む。量子アニーリングは、量子システム４２０が冷却ユニット４１０によって実質的な動作温度に維持されている間に実行してもよい。

量子アニーリングは、徐々に、例えば断熱的に、初期ハミルトニアンＨ^ｉｎｉｔから最終ハミルトニアンＨ^{ｆｉｎａｌ}＝Ｈ^ｐｒｏｂ＋Ｈ^ＳＲに移行するステップを含む。Ｈ^ＳＲは、補間ハミルトニアンＨ（ｔ）を経由した短距離ハミルトニアンである。図１に例示した実施形態によれば、補間ハミルトニアンは、式Ｈ（ｔ）＝Ａ（ｔ）Ｈ^ｉｎｉｔ＋Ｂ（ｔ）Ｈ^ｐｒｏｂ＋Ｃ（ｔ）Ｈ^ＳＲを有する。Ｈ^ｉｎｉｔは初期ハミルトニアン、Ｈ^{ｆｉｎａｌ}は最終ハミルトニアンを表し、ｔは時間パラメータで、Ａ（ｔ）、Ｂ（ｔ）及びＣ（ｔ）は時間パラメータｔに応じた補間係数である。ｔが初期時間ｔ_０の場合、補間係数Ａ（ｔ_０）は初期値１に等しくてもよく、補間係数Ｂ（ｔ_０）は初期値０に等しくてもよい。或いは、ｔが初期時間ｔ_０の場合、補間係数Ａ（ｔ_０）は補間係数Ｂ（ｔ_０）よりはるかに大きくてもよい。ｔが最終時間ｔ_ｆｉｎの場合、補間係数Ａ（ｔ_ｆｉｎ）は最終値０に等しくてもよく、Ｂ（ｔ_ｆｉｎ）及びＣ（ｔ_ｆｉｎ）はそれぞれ最終値１に等しくてもよい。これにより、補間ハミルトニアンＨ（ｔ_ｆｉｎ）がＨ^{ｆｉｎａｌ}に等しくなる。或いは、ｔが最終時間ｔ_ｆｉｎの場合、補間係数Ａ（ｔ_ｆｉｎ）は補間係数Ｂ（ｔ_ｆｉｎ）及びＣ（ｔ_ｆｉｎ）よりもはるかに小さい。量子アニーリングの実行は、徐々に、例えば断熱的に、補間係数Ａ（ｔ）、Ｂ（ｔ）及びＣ（ｔ）を初期時間での初期値から最終時間での最終値へ変換するステップを含んでいてもよい。したがって、補間ハミルトニアンは、初期時間での初期ハミルトニアンから最終時間での最終ハミルトニアンに徐々に変換される。特に、上述の量子アニーリングの手順は、常にＣ（ｔ）＝Ｂ（ｔ）を満たすように行われる。

特定の理論に縛られることなく量子力学の断熱定理等を考慮すると、初期ハミルトニアンから最終ハミルトニアン４７０への移行が十分に時間をかけて実行される場合、初期時間から最終時間までの範囲の時間パラメータｔの全ての値に対し、量子システム４２０の量子状態は、基底状態になるか、又は少なくとも補間ハミルトニアンＨ（ｔ）の基底状態にうまく近似される。したがって、量子アニーリングは、初期時間での初期量子状態を最終時間での最終量子状態に進行させる。最終量子状態は、最終ハミルトニアンの基底状態であるか、又は少なくとも最終ハミルトニアン４７０の基底状態にうまく近似される。

図１は、量子システム４２０の測定に適した測定装置４４０を更に示す。図示のように、測定装置４４０は、複数のキュービット１００のうちの一部分４２５を測定するのに適するものであってもよい。測定装置４４０を使用して前記一部分４２５を測定することで、最終量子状態の読み出しを取得することができる。最終ハミルトニアンの基底状態にうまく近似される最終量子状態は、計算問題４５２の解に関する情報を含む。最終量子状態の読み出しが、解に関する情報を明らかにすることができる。図１に示す実施形態によれば、図１に矢符４４５で示すように、読み出しは測定装置４４０から古典的計算システム４５０に提供される。古典的計算システム４５０は、その読み出しから計算問題の解４９０を決定することができる。古典的計算システム４５０は、少なくとも計算問題の試行解（trial solution）を決定し、その試行解が実際に計算問題の解になっているのかを検証する。ＮＰ問題の場合、その検証は多項式時間で実行可能な計算であり、通常容易に計算できる。計算問題の解が見つからないことが分かると、計算問題の解が見つかるまでそのプロセスが繰り返される。

上記の観点から、一つの実施形態によれば、複数のキュービットを含む量子システムを用いて計算問題の解を計算する方法を提供する。前記方法は、図８にボックス５１０で示すように、計算問題を量子システムの問題ハミルトニアンにコード化するステップを含む。問題ハミルトニアンは、複数の調整可能なパラメータを含む単体ハミルトニアンであり、コード化ステップは、計算問題から複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む。前記方法は、図８にボックス５２０で示すように、量子システムを初期量子状態から量子システムの最終ハミルトニアンの基底状態に向けて移行させるステップを更に含む。最終ハミルトニアンは、問題ハミルトニアン及び短距離ハミルトニアンの総和であり、問題ハミルトニアンの複数の調整可能なパラメータは、問題コード化設定に存在する。いくつかの実施形態において、短距離ハミルトニアンは、ｄが計算問題から独立しているｄ体ハミルトニアンである。前記方法は、図８にボックス５３０で示すように、複数のキュービットの少なくとも一部分を測定することで、量子システムの読み出しを取得するステップを更に含む。前記方法は、図８にボックス５４０で示すように、その読み出しから計算問題の解を決定するステップを更に含む。

したがって、実施形態によれば、量子システムを用いて、ＮＰ困難問題等の計算問題の解を決定することが可能になる。古典的計算システムのみを用いて計算問題の解を決定する、つまり量子システムを用いずに計算問題の解を決定する場合と比較すると、実施形態では、計算問題を解くのに必要な計算時間が短縮される。換言すれば、古典的計算システムと比較すると、実施形態では、計算問題をより速く解くことが可能になり、又は古典的計算システムで計算するには時間がかかり過ぎるような解であっても見つけることが可能となるかもしれない。

更なる利点は、問題ハミルトニアンが単体ハミルトニアンである側面に関する。他の種類の問題ハミルトニアン、特にキュービットの大きなグループ間での相互作用を伴う、又は遠く離れたキュービット間での相互作用（長距離の相互作用）を伴う問題ハミルトニアンは、実行不可能であるか、又は少なくとも量子システム及び量子計算を行う構成要素の非常に複雑なセットアップを必要とする。一方、上述の単体の問題ハミルトニアンは、もっと単純なセットアップ、つまりもっと単純な量子処理装置を用いることにより実現可能である。さらに、実施形態の調整可能なパラメータを有する問題ハミルトニアンは、幅広い計算問題をコード化することができる完全にプログラム可能なシステムを提供する。したがって、実施形態に係る装置及び方法は、ＮＰ困難問題等の幅広い計算問題の解を計算することができる。問題ハミルトニアンが必要とする特定の相互作用がシステムのハードウェアに組み込まれているため、限られた数の問題しかコード化できないシステムと比較すると、柔軟性が増した非常に強力な装置及び方法が提供される。

更なる利点は、最終ハミルトニアンが問題ハミルトニアン及び短距離ハミルトニアンの総和である側面に関する。短距離ハミルトニアンは、被加数ハミルトニアン（summand Hamiltonian）の総和であってもよい。本明細書に記載のように、被加数ハミルトニアンは制約ハミルトニアン（constraint Hamiltonian）であってもよい。短距離ハミルトニアンを有することで、遠く離れたキュービット間の相互作用を操作する必要が無いという利点が得られる。これは、量子システムで実現するには実行不可能な、又は少なくとも非常に複雑な量子処理装置のセットアップを必要とするような長距離の相互作用を必要とするハミルトニアンと大いに異なる。

短距離ｄ体ハミルトニアンのパラメータｄが計算問題から独立している場合、これは、どの計算問題をコード化するかに関わらず、同じ量子処理装置で計算可能であることを意味する。短距離ハミルトニアンが計算問題から独立している場合、異なる計算問題に対し、短距離ハミルトニアンによって決定されるキュービット間の相互作用を変更する必要が無いという更なる利点が得られる。

実施形態は、計算問題の解の計算に対し拡張可能なアーキテクチャを提供する。所定の量子システムでは、特定の最大サイズの多種多様な計算問題の解を計算することができる。最大サイズは、量子システムのキュービットの数により決定される。この最大サイズを超える計算問題の解を計算するために、より大きな量子システム、つまり数多くのキュービットを有する量子システムを、実施形態の対応する問題ハミルトニアン、短距離ハミルトニアン及び最終ハミルトニアンと共に提供することで、より大きなサイズの計算問題を処理することができる。適度に多くのキュービットを有する量子システムを選択することにより、どのような所望サイズの計算問題でもその解を計算できる。実施形態によれば、量子システムのキュービットの数に関わらず、問題ハミルトニアンは単体ハミルトニアンであり、最終ハミルトニアンは問題ハミルトニアン及び短距離ハミルトニアンの総和である。したがって、計算問題の解を計算するために拡張可能なアーキテクチャが提供される。

いくつかの実施形態によれば、計算問題は決定問題であってもよい。決定問題とは、イエスかノーかを問う質問が定式化された計算問題のことをいう。決定問題の解は、「イエス」か「ノー」かのいずれかである。或いは、決定問題の解は、１つの古典ビット、つまり０か１のいずれかであってもよい。他の実施形態によれば、計算問題は決定問題と異なる態様で定式化されてもよい。

計算問題は、例えば、コンピュータ科学、物理、化学又は工学の分野で考えられる種々の計算問題のいずれか１つであってもよい。範囲を限定する意図ではなく、説明目的で、以下の３つの計算問題の例について検討する。以下に検討する３つの例は決定問題の例である。

実施形態の計算問題の第１の例は「巡回セールスマン問題」である。巡回セールスマン問題は、都市の第１リスト、及びその第１リストに挙げられた各２都市間の距離の第２リストを伴う。巡回セールスマン問題では、次の問題について質問される。「第１リスト、第２リスト及び定数Ｋが与えられた場合に、セールスマンが（i）第１リストに挙げられた各都市を必ず１回ずつ巡回し、（ii）巡回開始時の都市に戻るという条件で、最大Ｋの巡回距離は存在するのか？」

実施形態の計算問題の第２の例は、数学的グラフの彩色に関する「３彩色問題」である。数学的グラフは頂点のセットと２つの頂点を結ぶ辺のセットとを含む。数学的グラフの３彩色では、数学的グラフの各頂点が３色（例えば「赤」、「緑」、「青」）のうちの１色に割り当てられる。この際、辺で結ばれるいずれの２つの頂点にも異なる色が割り当てられる。いくつかの数学的グラフでは、３彩色は存在しないかもしれない。３彩色問題では、次の問題について質問される。「数学的グラフが与えられた場合に、３彩色は存在するのか？」

実施形態の計算問題の第３の例は、イジングスピンモデルに関する。イジングスピンモデルは、複数のスピンｓ_１，ｓ_２，・・・，ｓ_ｎ間の相互作用を表す物理モデルである。各スピンｓ_ｉは１か−１のどちらかの値を有する変数であり、ｉは１からｎまでの値である。複数のスピンについては、イジングエネルギー関数Ｈ（ｓ_１，ｓ_２，・・・ｓ_ｎ）が考えられ、このイジングエネルギー関数は次の式を有する。
Ｈ（ｓ_１，ｓ_２，・・・ｓ_ｎ）＝Σ_ｉｊｃ_ｉｊｓ_ｉｓ_ｊ＋Σ_ｉｃ_ｉｓ_ｉ

各ｃ_ｉｊは結合係数で、各ｃ_ｉはフィールド係数である。イジングエネルギー関数は、対相互作用（pair-wise interactions）を伴う。スピンｓ_ｉ及びｓ_ｊ間の対相互作用はイジングエネルギー関数中の項ｃ_ｉｊｓ_ｉｓ_ｊで表される。結合係数ｃ_ｉｊの絶対値は、スピンｓ_ｉ及びｓ_ｊ間の対相互作用の強さを反映する。結合係数ｃ_ｉｊの符号は、例えば強磁性の又は反強磁性の相互作用等の対相互作用の特性を反映する。イジングスピンモデルは、長距離イジングスピンモデルであってもよい。長距離イジングスピンモデルは、距離測度に応じて互いに離れた一対のスピン同士の相互作用を含んでもよい。長距離イジングスピンモデルは、少なくとも２つのスピン間の最大距離の対数である距離を隔てて互いに離れたスピン対間の相互作用を含んでもよい。全対全イジングスピンモデル（all-to-all Ising spin model）等のいくつかの長距離イジングスピンモデルは、全てのスピン対間での相互作用を伴ってもよい。例えば、各結合係数ｃ_ｉｊが零でないイジングスピンモデルが長距離イジングスピンモデルであると考えられる。

イジングエネルギー関数は、スピンｓ_ｉ及びスピンｓ_ｉに影響するが他のスピンには影響しない外部フィールドとの相互作用を表す項ｃ_ｉｓ_ｉを更に含む。スピンｓ_ｉに影響するフィールドの強さと方向は、それぞれフィールド係数ｃ_ｉの絶対値と符号で表される。イジングスピンモデルに関係する計算問題（ここでは、イジングスピンモデル問題という）は、次のように定式化できる。「結合係数ｃ_ｉｊのセット、フィールド係数ｃ_ｉのセット及び定数Ｋを与えた場合に、Ｈ（ｓ_１，ｓ_２，・・・ｓ_ｎ）がＫより小さくなるようなスピンの形態（ｓ_１，ｓ_２，・・・ｓ_ｎ）は存在するのか？」

実施形態によれば、計算問題は複数の入力変数を含んでいてもよい。その複数の入力変数は、解くべき計算問題に関する情報を表してもよい。例えば、上述の３つの計算問題の例を参照すると、複数の入力変数は、巡回セールスマン問題の場合、都市の第１リスト及び距離の第２リストであり、３彩色問題の場合、グラフの頂点及び辺のセットであり、イジングスピンモデル問題の場合、結合係数ｃ_ｉｊとフィールド係数ｃ_ｉのセットである。

実施形態によれば、計算問題の解の計算は、計算問題に対する試行解の計算を含んでいてもよい。試行解は、計算問題の実際の解かもしれないし、そうでないかもしれない。複雑性クラスＮＰ（complexity class NP）に属する計算問題に係る実施形態の場合、以下に記載するように、計算問題の解の計算が証拠変数のセットの計算を含んでいてもよい。

実施形態によれば、計算問題の解は、複数のキュービットを有する量子システムを用いて計算される。複数のキュービットは、少なくとも８つのキュービットであり、特に少なくとも３つのキュービットであってもよい。更に、又は或いは、複数のキュービットは、Ｎ個のキュービットを含んでいてもよい。Ｎは、１００〜１００００で、好ましくは、これよりも大きい。図示されている複数のキュービット１００は、図で表現及び説明する目的で示しているだけであって、キュービットの実際の数はそれとは異なることが理解できるであろう。

量子システムのキュービットは、平面又は曲面を含む２次元表面又は３次元表面に配置されてもよい。図２〜４は、実施形態の複数のキュービット１００の異なる空間的配置を示す。これらの空間的配置は、例えばキュービット及び／又は他の異なる量子システム（キュートリットのようなｑ段階システム）が統合された量子チップ等の量子計算装置のレイアウトであってもよい。図２に示すように、複数のキュービット１００は、図２に破線で示す２次元の平面的な表面１１０に沿って配置されてもよい。図２に示す２次元表面１１０は、複数のキュービットの２次元の空間的配置を視覚的に表現することが目的であって、２次元表面１１０は複数のキュービット１００が配置される物理的に触れることのできる表面である必要がないことを理解できるであろう。以下に記載する２次元格子又は３次元格子に沿って配置された複数のキュービットの実施形態にも同様の考えが適用される。

図３に示す更なる実施形態によれば、複数のキュービット１００は、破線で示すような２次元格子１２０に沿って配置される。２次元格子又は３次元格子等の格子は、規則的なグリッドに沿って空間的に配置された複数のノードを含んでいてもよい。図３では、複数の黒色点で表わされた複数のキュービット１００が、２次元格子１２０のノードに対応する。図３のように、複数のキュービット１００の各キュービットは、２次元格子１２０のノードに配置されてもよい。図３に示す実施形態の例では、２次元格子１２０は２次元正方格子である。別の実施形態によれば、２次元格子１２０は、例えば六方格子若しくは三角格子、又は他の如何なる種類の２次元格子であってもよい。

実施形態によれば、複数のキュービットは３次元格子に沿って配置されてもよい。図３に関して上述したことと同様に、複数のキュービットは３次元格子のノードに対応していてもよい。複数のキュービットの各キュービットは、３次元格子のノードに配置されてもよい。３次元格子は３次元正方格子であってもよい。２次元格子の場合と同様に、他の種類の３次元格子も考えることができる。

２次元格子は平面的な構造であるので、３次元格子や他の不規則な空間的配置等に比べて、キュービットのより単純な空間的配置を提供できる。

実施形態によれば、複数のキュービットは２次元格子の一部分、又は３次元格子の一部分に沿って配置される。図４は、２次元格子の三角形部分１２１に沿って配置された複数のキュービット１００の実施形態の例を表している。図４は、三角形部分１２１の上面図である。三角形部分は、実施形態の方法を実行するように構成された実施形態に係る量子計算装置のレイアウトに対応する。異なる形状を有する格子の一部分も考えることができる。

問題ハミルトニアンは、複数の調整可能なパラメータを有する単体ハミルトニアンである。実施形態の量子システムの単体ハミルトニアンは、２つ以上のキュービットからなるグループ間で相互作用が生じないハミルトニアンであってもよい。単体ハミルトニアンは、複数の被加数ハミルトニアンの総和であってもよい。各被加数ハミルトニアンは、複数のキュービットのうちの１つのキュービットに作用する。単体ハミルトニアンは、式Ｈ＝Σ_ｉＨ_ｉを有していてもよい。各Ｈ_ｉはｉ番目のキュービットのみに作用する被加数ハミルトニアンである。単体ハミルトニアンは、複数のキュービットと磁場又は電場等の外部実体との間の相互作用を表してもよい。この場合、各キュービットは外部実体に個別に作用する。

図５は、実施形態の単体ハミルトニアンの概略図である。範囲を限定することなく具体性を持たせるため、図５に示す複数のキュービットは、図４に示すものと同様に、三角形を成す２次元正方格子の一部分に配置された１０個のキュービット、すなわちキュービット２０１〜２１０を含む。図５を参照して説明する単体ハミルトニアンは、１０個の被加数ハミルトニアン２２１〜２３０の総和である。図５において、各被加数ハミルトニアン２２１〜２３０は、各被加数ハミルトニアンが１つのキュービットに作用することを示すため、１つのキュービットを囲む正方形として概略的に図示している。例えば、被加数ハミルトニアン２２１がキュービット２０１には作用するが、残りのキュービット２０２〜２１０には作用しないことを示すため、被加数ハミルトニアン２２１は、キュービット２０１を囲む正方形及びキュービット２０１で表わされる。

単体ハミルトニアンである問題ハミルトニアンは、上述の被加数ハミルトニアンの総和であってもよい。問題ハミルトニアンの複数の調整可能なパラメータは、被加数ハミルトニアンの複数の調整可能なパラメータを含んでいてもよい。単体ハミルトニアンの１つ以上の被加数ハミルトニアン、特に各被加数ハミルトニアンは、１つ以上の調整可能なパラメータを含んでいてもよい。

実施形態の問題ハミルトニアンの調整可能なパラメータは、複数のキュービットのうちの１つのキュービットと外部実体との相互作用の強さ及び／又は方向を表すパラメータであってもよい。外部実体は、１つ以上の磁場、１つ以上の電場、及び／又は１つ以上のレーザフィールド、マイクロ波、機械的変形からの位相シフトのうちの少なくとも１つを含んでいてもよい。外部実体を調整することにより、及び／又はキュービットと外部実体との相互作用の強さ及び／又は種類を調整することにより、問題ハミルトニアンの調整可能なパラメータを調整することを実現できる。したがって、調整可能なパラメータとは、量子システムのハードウェアに組み込まれていない相互作用等の調整可能な相互作用を表してもよい。

他の実施形態と組み合わせ可能な実施形態によれば、問題ハミルトニアンの複数の調整可能なパラメータは、複数のキュービットに作用する単体フィールドの複数のフィールド強度及び／又は複数のフィールド方向を含んでもよい。複数のキュービットに作用する前記フィールドは、例えば、超電導キュービットに関する実施形態では、１つ以上の磁場及び／又は１つ以上の電場を含んでいてもよい。

単体フィールドとは、複数のキュービットのうちの１つのキュービットに影響するフィールドである。実施形態によれば、複数の単体フィールドは、異なり得るフィールド強度及び／又は異なり得るフィールド方向に従った対応するキュービットに影響を与える異なる単体フィールドを含んでいてもよい。例えば、第１の単体フィールドと第２の単体フィールドは、それぞれ複数のキュービットのうちの第１のキュービットと第２のキュービットに影響する。この場合、第１の単体フィールドと第２の単体フィールドは、両方とも磁場等であり、異なるフィールド強度及び／又はフィールド方向を有していてもよい。

他の実施形態と組み合わせ可能な実施形態によれば、単体ハミルトニアンは式Σ_ｋＪ_ｋσ_ｚ ^（ｋ）を有する。σ_ｚ ^（ｋ）は複数のキュービットのｋ番目のキュービットのパウリ演算子であり、各Ｊ_ｋは係数であり、係数Ｊ_ｋは単体ハミルトニアンの複数の調整可能なパラメータを形成する。いくつかの実施形態によれば、パウリ演算子σ_ｚ ^（ｋ）は第１の空間方向に関連するパウリ演算子であってもよい。

複数の超電導キュービットを有する量子システムの場合、問題ハミルトニアン等の単体ハミルトニアンは、複数の超電導キュービットと相互に作用する複数の磁束により実現できる。磁束又は磁束バイアスは、超電導キュービットの一次超電導ループを通過し、二次超電導ループも通過して広がる。問題ハミルトニアンの複数の調整可能なパラメータは、複数の磁束又は磁束バイアスを調整することにより調整可能である。

トラップイオンで実現される量子システムの場合、各イオンを空間的に分離するか、エネルギーで分離することにより対処可能である。空間的分離の場合は、音響光学偏向器、音響光学変調器、マイクロミラーデバイス等を通過した及び／又はこれらから反射したレーザビームの使用を伴う。エネルギーでの分離は、内部遷移周波数を変化させる磁場勾配の使用を伴い、そうすることで、印加フィールドの離調（detuning）等のエネルギーの相違を通じて選択が可能になる。単体ハミルトニアンは、内部遷移とのレーザフィールド若しくはマイクロ波の共鳴又は無共鳴によって、又は空間的な磁場の相違によって実現できる。

量子ドットで実現される量子システムの場合、単体ハミルトニアンは電場で実現できる。

ＮＶセンターで実現される量子システムの場合、マイクロ波パルスを適用した磁気共鳴を用いて、キュービットの状態をナノ秒の時間尺度でコヒーレントに操作可能である。キュービットの状態の選択的な操作は、すぐ近くの核スピンの状態次第で達成できる。

計算問題は、問題コード化設定上にマッピングされてもよい。問題コード化設定は、計算問題に関する情報に依存していてもよいし、及び／又は、その情報を含んでいてもよい。問題コード化設定を決定する動作は、複数の調整可能なパラメータのそれぞれの値を決定する及び／又は計算することを含んでいてもよい。それぞれの値は、計算問題から決定及び／又は計算される。

実施形態によれば、異なる計算問題は、対応する異なる問題コード化設定を決定することにより問題ハミルトニアンにコード化される。例えば、第１の計算問題及び第２の計算問題は、問題ハミルトニアンにコード化され、これらは、複数の調整可能なパラメータに対する第１の問題コード化設定と第２の問題コード化設定となる。第２の計算問題が第１の計算問題と異なる場合、調整可能なパラメータの第２の問題コード化設定は第１の問題コード化設定と異なってもよい。

実施形態によれば、上記方法は計算問題又は少なくともその計算問題に関する情報を、図１に示す古典的計算システム４５０等の古典的計算システムに提供することを含んでもよい。例えば、上述の計算問題の複数の入力変数は、古典的計算システムに提供されてもよい。実施形態によれば、計算問題は、古典的計算システムによって問題コード化設定にコード化されてもよい。古典的計算システムは、計算問題、例えば、計算問題の複数の入力変数から問題コード化設定を計算するように構成されていてもよい。

本明細書で使用する用語「短距離ハミルトニアン」は、複数のキュービットの相互作用を表すハミルトニアンのことである。短距離ハミルトニアンでは、相互作用遮断距離よりも長い距離で互いに離れているキュービット間で相互作用が生じない。相互作用遮断距離は一定の距離であってもよい。相互作用遮断距離は、複数のキュービットの中のキュービット間の最大キュービット距離と比べて非常に短くてもよい。例えば、相互作用遮断距離は、最大キュービット距離の３０％以下、具体的には２０％以下、より具体的には１０％以下であってもよい。格子に沿って配置された複数のキュービットの場合、短距離ハミルトニアンはｒ範囲のハミルトニアンであってもよい。ｒ範囲のハミルトニアンでは、格子の基本距離（格子定数）のｒ倍よりも長い距離で互いに離れているキュービット間で相互作用が生じない。ｒは１〜５であり、例えばｒ＝２^１／２，２，３，４又は５である。実施形態の格子の基本距離についての考え方は、図６及び７等を参照して以下に説明する。

量子システムのキュービットの数に関わらず、量子システムのプラケットハミルトニアン及び最近接対ハミルトニアン（pairwise nearest-neighbor Hamiltonian）は、短距離ハミルトニアンと見なされる。

短距離ハミルトニアンの例としては、単体ハミルトニアンが挙げられる。単体ハミルトニアンの場合、２つ以上のキュービットからなるグループ間では相互作用が生じず、各キュービットと磁場や電場等の外部実体との間でしか相互作用が生じないため、相互作用遮断距離は零であると考えられる。

図６及び７は、複数のキュービット１００が、２次元正方格子１２０に沿って配置され、２次元正方格子の三角形部分を成す２次元正方格子のノードの位置に位置する実施形態について、短距離ハミルトニアンの更なる例を示す。範囲を限定することなく具体性を持たせるため、図６及び７で例示される２次元正方格子１２０は、１０行１０列からなる１０×１０正方格子内に三角形状に配置された５５個のキュービットを含む。例えば、点線で示す行３９１のように、Ｘ方向３１０に沿った２次元正方格子１２０のキュービットのいずれかの行を検討した場合、その行において連続するキュービットは、互いに基本距離Ｄ分だけ離れて配置される。これをＸ方向の格子定数ともいう。基本距離Ｄは、符号３５０で示される。同様に、例えば、点線で示す列３９２のように、Ｙ方向３２０に沿った２次元正方格子１２０のキュービットのいずれかの列を検討した場合、その列において連続するキュービットは、互いに基本距離分だけ離れて配置される。これをＹ方向の格子定数ともいう。図６及び７において、正方格子の格子並びにＸ方向及びＹ方向の基本距離（格子定数）は同じである。しかしながら、Ｘ方向及びＹ方向の格子定数は異なることもある。図示のように、Ｘ方向３１０はＹ方向３２０に対して直交している。図６及び図７に示す複数のキュービット１００の最大キュービット距離は、キュービット３０１及び３０２の間の距離である。最大キュービット距離は、（９√２）Ｄに等しい。

図６を参照して説明する短距離ハミルトニアンの例は、最近接対ハミルトニアンである。最近接対ハミルトニアンは、２次元格子１２０上の対の近接するキュービット間での相互作用のみを伴う。一対の近接するキュービットとは、基本距離Ｄ分だけ互いに離れた一対のキュービットのことをいう。図６に示すキュービット３６２及び３６４は、一対の近接するキュービットの例になっている。最近接対ハミルトニアンは、複数の被加数ハミルトニアンの総和であってもよく、この場合、各被加数ハミルトニアンは、一対の近接するキュービット間の相互作用を表す。図６を参照して説明する最近接対ハミルトニアンの場合、相互作用遮断距離は基本距離Ｄに等しい。したがって、相互作用遮断距離は、最大キュービット距離に比べると非常に短い。すなわち、相互作用遮断距離Ｄは、最大キュービット距離の１０％未満である。

図７を参照して説明する短距離ハミルトニアンの例は、プラケットハミルトニアンである。図７において、黒色円で示す５５個のキュービットが再び２次元正方格子１２０に配置され、三角形を形成する。図７において符号３７０で示すように、２次元正方格子１２０のプラケットは、２次元正方格子１２０の基本の正方形である。プラケット３７０は、キュービット３７１、３７２、３７３及び３７４を有する。キュービット３７１は、キュービット３７２及び３７４から基本距離Ｄ分だけ離れて配置され、キュービット３７３もまた、キュービット３７２及び３７４から基本距離Ｄ分だけ離れて配置されている。さらに、キュービットのプラケットを完成させるために、黒色の矩形で示す予備キュービット（auxiliary qubit）が新たな線上に追加される。例えば、予備キュービット３０５は、キュービット３０２、３０３及び３０４のプラケットを完成させる。予備キュービットは、特定の量子状態、例えば｜１＞で準備できる。この格子の形状では、プラケットハミルトニアンは、２次元正方格子１２０のプラケットに対応する４つのキュービットからなるグループ間、又は３つのキュービット及び１つの予備キュービットからなるグループ間での相互作用のみを伴う。プラケットハミルトニアンは、複数の被加数ハミルトニアンの総和であってもよい。各被加数ハミルトニアンは、格子上のキュービットのプラケットに対応する相互作用、又はキュービット及び予備キュービットのプラケットに対応する相互作用を表してもよい。或いは、予備キュービットを１つも使用せず、プラケットハミルトニアンが３つのキュービットのみの間での相互作用を示す被加数ハミルトニアンを含む。図７を参照して説明するプラケットハミルトニアンでは、プラケットの２つのキュービット間の最大距離が√２Ｄであるため、相互作用遮断距離は√２Ｄである。例えば、キュービット３７１及び３７３の間の距離は√２Ｄである。したがって、相互作用遮断距離は、最大キュービット距離に比べると非常に短い。すなわち、相互作用遮断距離√２Ｄは、最大キュービット距離の１２％未満である。

複数の超電導キュービットを有する量子システムの場合、プラケットハミルトニアンは、複数の補助キュービット（ancillary qubit）を用いることで実現できる。補助キュービットは、例えば各プラケットの中央等、各プラケットの内側に配置される。式Ｋ_ｋｍσ_ｚ ^（ｋ）σ_ｚ ^（ｍ）のキュービット間の相互作用は、誘導結合ユニット等の結合ユニットにより実現できる。結合ユニットは、超電導量子干渉装置を含む。調整可能な磁束バイアスを超電導量子干渉装置に適用することで、係数Ｋ_ｋｍを調整することが可能になる。そして、プラケットハミルトニアンの被加数ハミルトニアンが、量子基本状態の｜０＞及び｜１＞間の強制エネルギーの差異に対応する式σ_ｚ ^（ｋ）σ_ｚ ^（ｍ）と単体のσ_ｚ ^（１）の対相互作用のみを含むＨ_ｓｒ，ｐ＝Ｃ（σ_ｚ ^（１）＋σ_ｚ ^（２）＋σ_ｚ ^（３）＋σ_ｚ ^（４）−２σ_ｚ ^（ｐ）−１）^２により実現可能である。σ_ｚ ^（ｐ）は補助キュービットを表す。短距離ハミルトニアンは、被加数ハミルトニアンＨ_ｓｒ，ｐの総和である。補助キュービットを伴う実施形態では、複数の補助キュービットに対し、式ｈΣ_ｐσ_ｘ ^（ｐ）の単体ハミルトニアンが初期ハミルトニアンに追加される。

或いは、プラケットハミルトニアンは、補助キュービットを用いることなく、例えばトランズモンキュービットとして３島型超伝導装置（three-island superconducting devices）を用いて実現可能である。２つの追加の超伝導量子干渉装置を結合ユニットに統合し、プラケットの４つのキュービットをコプレーナ共振器に容量的に結合することにより、式−Ｃσ_ｚ ^（１）σ_ｚ ^（２）σ_ｚ ^（３）σ_ｚ ^（４）の被加数ハミルトニアンが実現可能である。結合係数Ｃは、２つの追加の超伝導量子干渉装置を通じて、時間依存性の磁束バイアスにより調整可能である。

トラップイオンにより実現される量子システムの場合、２つのイオン間の相互作用はフォノンバスで送信される。この場合、レーザ又はマイクロ波が使用され、フォノンのブルーサイドバンド遷移及び／又はレッドサイドバンド遷移で離調される。レーザの強度及び離調によって、相互作用強度が調整可能になる。リュードベリ励起を通じた直接相互作用も利用可能である。

低温原子により実現される量子システムの場合、キュービット間の相互作用は、ｄ原子を励起するレーザの離調によって制御可能である。この場合、ハミルトニアンはｄ体ハミルトニアンである。プラケットハミルトニアンは、ｄ体相互作用又は２体相互作用を有する補助キュービットのいずれから実行されてもよい。

量子ドットにより実現される量子システムの場合、２つのキュービット間の相互作用は、電場勾配及び磁場により調整される。短距離ハミルトニアンは、パルスシーケンス及び磁場により実現可能である。プラケットハミルトニアンは、追加の補助キュービットを、プラケットの全ての対に作用する短距離ハミルトニアンと共に用いることで実現される。

ＮＶセンターにより実現される量子システムの場合、ＮＶセンター間の相互作用は、これらを光照射野（light field）に結合することにより伝送される。

他の実施形態と組み合わせ可能な実施形態によれば、複数のキュービットは、２次元格子に沿って配置される。短距離ハミルトニアンは、２次元格子のプラケットに対応する４つのキュービットのグループ間での相互作用を伴う。実施形態によれば、短距離ハミルトニアンはプラケットハミルトニアンであってもよい。

他の実施形態と組み合わせ可能ないくつかの実施形態によれば、短距離ハミルトニアンはｄ体ハミルトニアンである。ｄは２、３、４、５、６、７又は８である。ｄ体ハミルトニアンは、複数のキュービットの相互作用を表すハミルトニアンであってもよい。この場合、ｄ＋１以上のキュービットからなるグループ間での共同相互作用（joint interactions）は生じない。ｄ体ハミルトニアンは、ｄ以下のキュービットからなるグループ間での相互作用を伴ってもよい。ｄ体ハミルトニアンは、複数の被加数ハミルトニアンの総和であってもよい。この場合、各被加数ハミルトニアンは、ｄ以下のキュービットからなるグループ間での共同相互作用を表す。

例えば、単体ハミルトニアンは、ｄ＝１であるｄ体ハミルトニアンと考えることができる。更なる例として、最近接対ハミルトニアンは、ｄ＝２であるｄ体ハミルトニアンと見なされてもよい。更なる例として、プラケットハミルトニアンは、ｄ＝４であるｄ体ハミルトニアンと見なされてもよい。他の実施形態と組み合わせ可能な実施形態によれば、短距離ハミルトニアンは、ｄ＝４であるｄ体ハミルトニアンであってもよい。ｄの値は格子の形状に依存してもよい。例えば、六方格子の場合、プラケットは６つのキュービットを伴い、プラケットハミルトニアンは６体ハミルトニアンとなり得る。

ｄ＝４等の小さな値のｄを有するｄ体ハミルトニアンである短距離ハミルトニアンを有することは有利である。より大きな値のｄを有するｄ体ハミルトニアンに比べて、対応するキュービット間の相互作用をより容易に設計できるからである。

計算問題のサイズは、計算問題を特定するのに必要とされる古典的情報単位の数の大きさであってもよい。計算問題のサイズは、計算問題の入力変数の数に依存してもよい。計算問題のサイズは、入力変数の数が増えるに従い増加してもよい。計算問題のサイズは、入力変数の数と等しくてもよい。例えば、巡回セールスマン問題の場合、前記サイズは、第１リスト及び第２リストの長さの総和である。更なる例として、イジングスピンモデル問題の場合、前記サイズは、スピンＳ_ｉの数ｎである。

第１のサイズを有する第１の計算問題の場合、対応する最終ハミルトニアンは、第１の問題ハミルトニアン及び第１の短距離ハミルトニアンの総和であってもよい。第２のサイズを有する第２の計算問題の場合、対応する最終ハミルトニアンは、第２の問題ハミルトニアン及び第２の短距離ハミルトニアンの総和であってもよい。第２のサイズが第１のサイズと同じである場合、第２の短距離ハミルトニアンは第１の短距離ハミルトニアンと同じであってもよい。第２のサイズが第１のサイズと異なる場合、第２の短距離ハミルトニアンは第１の短距離ハミルトニアンと異なってもよい。例えば、イジングスピンモデルに関し、第１の計算問題は、結合係数及びフィールド係数の第１セットを有するＮスピンに対する第１のイジングスピンモデル問題であり、第２の計算問題は、第１セットの結合係数及びフィールド係数とは異なる結合係数とフィールド係数の第２セットを有するＮスピンに対する第２のイジングスピンモデル問題である。この場合、第１及び第２の双方のイジングスピンモデル問題が、数Ｎに等しいと考えられる。実施形態によれば、第１のイジングスピンモデル問題の短距離ハミルトニアンは、第２のイジングスピンモデル問題の短距離ハミルトニアンと同じである。

他の実施形態と組み合わせ可能な実施形態によれば、短距離ハミルトニアンは、ｄ体ハミルトニアンであってもよい。この場合、ｄは計算問題から独立していてもよい。更に、相互作用遮断距離は、計算問題から独立していてもよい。他の実施形態と組み合わせ可能な実施形態によれば、短距離ハミルトニアンは、計算問題から独立していてもよい。

実施形態によれば、前記方法は、量子システムを初期量子状態から最終ハミルトニアンの基底状態に向けて進行させるステップを含む。最終ハミルトニアンの基底状態は、最終ハミルトニアンに対しエネルギーを最小化する量子システムの量子状態である。最終ハミルトニアンの基底状態は、最終ハミルトニアンの固有状態であり、特に最小固有値を有する固有状態である。計算問題は問題ハミルトニアンにおいてコード化されるため、且つ、最終ハミルトニアンは問題ハミルトニアン及び短距離ハミルトニアンの総和であるため、最終ハミルトニアンの基底状態は、計算問題の情報を含み、及び／又は計算問題の解をコード化することができる。

最終ハミルトニアンの基底状態は、零度での量子システムの状態であってもよい。特定の理論に縛られることなく量子物理学分野の考え方に従えば、量子システムが絶対零度に達することは不可能であると考えられる。それでもやはり、量子システムを動作温度Ｔ_ｍａｘに冷却する等、初期量子状態から最終ハミルトニアンの基底状態に向けて進行させるステップは、最終ハミルトニアンの基底状態に近づくことを可能にする。動作温度Ｔ_ｍａｘは、量子システムに使用されるキュービットの種類に強く依存する。例えば、超電導キュービットの場合、Ｔ_ｍａｘは５０ｍＫ以下であり、好ましくは１ｍＫ以下であってもよい。最終ハミルトニアンの基底状態に近づけるため、量子システムを初期量子状態から最終量子状態に進行させてもよい。最終量子状態は、動作温度Ｔ_ｍａｘ又はより低温度にある量子システムの状態であってもよい。すなわち、最終量子状態は、動作温度又はより低温度にある最終ハミルトニアンの熱状態であってもよい。したがって、最終量子状態は、最終ハミルトニアンの基底状態を近似する。最終量子状態は、最終ハミルトニアンの基底状態に関する情報を含んでいてもよい。最終量子状態は、計算問題の解に関する情報を含んでいてもよい。

量子システムは、冷却ユニット等により動作温度Ｔ_ｍａｘ又はより低温度に冷却されてもよい。動作温度は、零度でない温度であってもよい。

他の実施形態と組み合わせ可能な実施形態によれば、前記方法は、量子システムを初期ハミルトニアンの基底状態に向けて冷却することにより、量子システムを初期量子状態に初期化するステップを含む。初期ハミルトニアンの基底状態は、初期ハミルトニアンに対するエネルギーを最小化する量子システムの量子状態である。初期ハミルトニアンの基底状態は、初期ハミルトニアンの固有状態であり、特に最小固有値を有する固有状態である。初期ハミルトニアンの基底状態は、零度での量子システムの状態である。初期ハミルトニアンの基底状態に向けて量子システムを冷却することで、初期ハミルトニアンの基底状態に近づくことが可能になる。初期量子状態は、初期ハミルトニアンの基底状態を近似する。

初期ハミルトニアンは、計算問題から独立していてもよい。初期ハミルトニアンは、ｄが１、２、３又は４のｄ体ハミルトニアンであってもよい。初期ハミルトニアンは、単体ハミルトニアンであってもよい。単体ハミルトニアンである初期ハミルトニアンを有することで、超電導キュービットの量子システム等、初期ハミルトニアンを実現するための単純なセットアップが可能になる。

初期ハミルトニアンは、式Ｈ^ｉｎｉｔ＝Σ_ｋａ_ｋσ_ｘ ^（ｋ）を有する単体ハミルトニアンであってもよい。ａ_ｋは複数のキュービットのｋ番目のキュービットの係数であり、σ_ｘ ^（ｋ）はｋ番目のキュービットに作用するパウリ演算子であってもよい。特に、σ_ｘ ^（ｋ）は、第２の空間方向に対応するパウリ演算子であってもよい。第２の空間方向は、第１の空間方向に対して直交していてもよい。パウリ演算子σ_ｘ ^（ｋ）及びパウリ演算子σ_ｚ ^（ｋ）は、非可換、特に反可換の演算子であってもよい。実施形態によれば、各係数ａ_ｋは単一の共通の係数ｈに等しい。初期ハミルトニアンは、式Ｈ^ｉｎｉｔ＝ｈΣ_ｋσ_ｘ ^（ｋ）を有する単体ハミルトニアンであってもよい。

超電導キュービットの場合、超電導キュービットの１次超伝導ループを通過する磁束バイアスは、基本状態の｜０＞及び｜１＞が同じエネルギーを有するように、つまり、これらの基本状態のエネルギーの差異が零になるように、設定されてもよい。更に、２次超伝導ループを通過する磁束バイアスは、基本状態の｜０＞及び｜１＞を結合可能である。したがって、超電導キュービットに対して、式ｈσ_ｘ ^（ｋ）の被加数ハミルトニアンが実現可能である。したがって、複数の超電導キュービットに対して、式Ｈ^ｉｎｉｔ＝ｈΣ_ｋσ_ｘ ^（ｋ）の初期ハミルトニアンが実現可能である。初期ハミルトニアンの基底状態は、係数ｈを量子システムの温度で決定されるエネルギー尺度よりも大きな値に設定することにより、ほぼ確実なものにすることができる。

トラップイオンで実現される量子システムの場合、イオンはレーザを用いた光ポンピングにより初期化できるので、イオンはキュービットの２つの量子基本状態のうちの１つに確実に移動する。これによりエントロピーが減少し、内部状態が冷却される。

低温原子で実現される量子システムの場合、初期量子状態は、大きな離調で基底状態にある原子をリュードベリ状態に励起することにより準備されてもよい。

ＮＶセンターで実現される量子システムの場合、ＮＶセンターは、標準共焦点光学顕微鏡技術（standard optical confocal microscopy technics）を用いて個別に対処してもよい。初期化及び測定は、非共鳴又は共鳴光学励起により実行される。

実施形態によれば、量子システムは、初期量子状態から最終ハミルトニアンの基底状態に向けて進行する。実施形態によれば、量子システムを進行させるステップは、量子アニーリングを実行するステップを含む。量子アニーリングは、プログラム可能な量子アニーリングユニットにより実行される。

量子アニーリングを実行するステップは、量子システムの初期ハミルトニアンから最終ハミルトニアンに移行するステップを含んでもよい。量子アニーリングを実行するステップは、補間ハミルトニアンを経由して初期ハミルトニアンから最終ハミルトニアンに移行するステップを含んでもよい。補間ハミルトニアンは、時間依存性のハミルトニアンであってもよい。補間ハミルトニアンは、初期ハミルトニアン及び最終ハミルトニアンの間を補間するために、１つ以上の補間パラメータを有してもよい。例えば、補間ハミルトニアンは、式Ｈ（ｔ）＝Ａ（ｔ）Ｈ^ｉｎｉｔ＋Ｂ（ｔ）Ｈ^{ｆｉｎａｌ}を有してもよい。Ｈ^ｉｎｉｔは初期ハミルトニアンであり、Ｈ^{ｆｉｎａｌ}は最終ハミルトニアンであり、ｔは時間パラメータであり、Ａ（ｔ）及びＢ（ｔ）は時間パラメータｔに依存する補間係数であってもよい。

量子アニーリングを実行するステップは、初期時間での初期ハミルトニアンから最終時間での最終ハミルトニアンに移行するステップを含んでもよい。初期時間では、補間ハミルトニアンの１つ以上の補間パラメータがそれぞれ１つ以上の初期値に設定されてもよい。初期時間では、補間ハミルトニアンは初期ハミルトニアンに等しくてもよい。例えば、補間ハミルトニアンが式Ｈ（ｔ）＝Ａ（ｔ）Ｈ^ｉｎｉｔ＋Ｂ（ｔ）Ｈ^{ｆｉｎａｌ}を有する実施形態では、補間パラメータＡ（ｔ）は初期時間での初期値１に設定され、補間パラメータＢ（ｔ）は、初期時間での初期値０に設定されてもよい。したがって、補間ハミルトニアンは初期時間でのＨ^ｉｎｉｔに等しい。

量子アニーリングを実行するステップは、初期ハミルトニアンから最終ハミルトニアンに徐々に移行するステップを含む。量子アニーリングを実行するステップは、補間ハミルトニアンの１つ以上の補間パラメータを徐々に変化させるステップを含む。初期時間及び最終時間の間にある中間時間では、補間ハミルトニアンは、初期ハミルトニアン及び／又は最終ハミルトニアンと異なる。補間ハミルトニアンの補間パラメータは、初期時間等の初期値から最終時間等の最終値に徐々に変化してよい。最終時間では、補間ハミルトニアンは最終ハミルトニアンに等しい。例えば、補間ハミルトニアンが式Ｈ（ｔ）＝Ａ（ｔ）Ｈ^ｉｎｉｔ＋Ｂ（ｔ）Ｈ^{ｆｉｎａｌ}を有する実施形態の場合、補間パラメータＡ（ｔ）は、初期時間での初期値１から最終時間での最終値０に徐々に変化してよい。同様に、補間パラメータＢ（ｔ）は、初期時間での初期値０から最終時間での最終値１に徐々に変化してよい。したがって、補間ハミルトニアンＨ（ｔ）は、最終時間での最終ハミルトニアンに等しい。

他の実施形態と組み合わせ可能な実施形態によれば、量子システムは、量子アニーリングを実行している間、５０ｍＫ以下、特に１ｍＫ以下の動作温度に維持される。

他の実施形態と組み合わせ可能な実施形態によれば、初期ハミルトニアンから最終ハミルトニアンへ量子アニーリングを実行するステップは、初期ハミルトニアンを最終ハミルトニアンへ断熱的に移行させるステップを含む。

式ｈΣ_ｋσ_ｘ ^（ｋ）の初期ハミルトニアン、式Σ_ｋＪ_ｋσ_ｚ ^（ｋ）の問題ハミルトニアン、及びプラケットに対応する被加数ハミルトニアンＣ_ｌを用いた式ＣΣ_ｌＣ_ｌの短距離ハミルトニアンを有する複数の超電導キュービットを含む量子システムの場合、量子アニーリングは以下のように実行される。Ｃ＝Ｊ_ｋ＝０である初期化の後、係数ａをａ＝０になるまで徐々に下げる間、Ｃ及びＪ_ｋを徐々に大きくする。

量子アニーリングを伴う実施形態に代えて、又はこれに付け加えて、最終ハミルトニアンの初期量子状態から基底状態に向けて量子システムを移行させるステップは、量子システムを初期量子状態から最終量子状態に冷却するステップを含む。初期量子状態は、初期温度での量子システムの状態であってもよい。最終量子状態は、最終温度での量子システムの状態であってもよい。最終温度は初期温度よりも低い。初期温度での初期量子状態から最終温度での最終量子状態に量子システムを冷却するステップは、初期温度から最終温度へ温度を下げる、例えば徐々に下げることを含んでもよい。他の実施形態と組み合わせ可能な実施形態によれば、最終温度は５０ｍＫ以下、特に１ｍＫ以下であってもよい。他の実施形態と組み合わせ可能な実施形態によれば、初期温度は室温以下、特に２００ケルビン以下であってもよい。

実施形態によれば、前記方法は、複数のキュービットの少なくとも一部分を測定することで最終量子状態の読み出しを取得するステップを含む。いくつかの実施形態によれば、複数のキュービットの一部分が測定され、複数のキュービットの全てが測定されるわけではない。複数のキュービットの一部分は、複数のキュービットの７０％以下であってもよく、具体的には６０％以下、より具体的には５０％以下であってもよい。いくつかの実施形態によれば、複数のキュービットの総数をＮとすると、一部分に含まれるキュービットの数は√Ｎとなる。

複数のキュービットの少なくとも一部分を測定するステップは、その少なくとも一部分における各キュービットを個別に測定するステップを含んでもよい。少なくとも一部分を測定するステップは、キュービットの少なくとも一部分における各キュービットに対するパウリ演算子σ_ｚ等のパウリ演算子を測定するステップを含んでもよい。少なくとも一部分を測定するステップは、複数のキュービットの少なくとも一部分における各キュービットに対し２通りの成果測定（two-outcome measurement）を行うステップを含んでもよい。２通りの成果測定は、０又は１など、２通り生じ得る結果の１つを提供する。キュービットの少なくとも一部分は、測定装置により測定される。

少なくとも一部分を測定するステップは、最終量子状態の読み出しを与える。読み出しは、複数の古典ビットで表される古典的な情報の式を有してもよい。読み出しは、最終量子状態及び／又は最終ハミルトニアンの基底状態に関する情報を明らかにできる。読み出しは、計算問題の試行解、実際の解、又は証拠変数のセット等、解に関する情報を提供できる。読み出しは、計算問題の解であってもよい。

複数のＮ個の超電導キュービットを有する量子システムの場合、複数のキュービットにおけるキュービットの状態｜０＞及び｜１＞を、複数の超電動量子干渉装置、特にＮ個のヒステリシスＤＣ超電動量子干渉装置、及び数が√Ｎであるバイアス線で制御されるＮ個のＲＦ超電動量子干渉装置ラッチを有する測定装置を用いて、高精度に測定可能である。

トラップイオンで実現される量子システムの場合、量子システムの測定を蛍光分光法により実行可能である。この場合、イオンが２つのスピン状態のうちの１つにあると、イオンは短い寿命で遷移する。その結果、駆動された状態のイオンは、他のイオンが暗いままいる間に、多くのフォトンを放出する。放出されたフォトンは、市販のＣＣＤカメラで記録可能である。蛍光分光前の適切な単一キュービットパルス（single-qubit pulse）により、ブロッホ球の任意の方向の測定が可能である。

低温原子で実現される量子システムの場合、キュービットは、基底状態の原子の選択的一掃（selective sweep）及び単一の位置分解能（single site resolution）での蛍光撮像を実行することで測定可能である。

量子ドットで実現される量子システムの場合、キュービットは、高速断熱通過（rapid adiabatic passage）によるパルスシーケンスから読み出し可能である。

実施形態によれば、前記方法は、読み出しから計算問題の解を決定するステップを含んでもよい。上記方法は、読み出しから解を計算するステップを含んでもよい。読み出しは、古典的計算システムに提供されてもよい。古典的計算システムは、その読み出しから計算問題の解を決定又は計算してもよい。

計算問題は、コンピュータ科学の分野で考えられる複雑性クラスＮＰに関連してもよい。「ＮＰ」は、「非決定性多項式時間」の略称である。他の実施形態と組み合わせ可能な実施形態によれば、計算問題は複雑性クラスＮＰに属する。複雑性クラスＮＰは、決定問題を含む。口語的に言えば、複雑性クラスＮＰに属する計算問題の場合、計算問題の解が「イエス」だということを検証可能であることに基づく証拠変数のセットが存在する。ＮＰの計算問題の場合、解が「イエス」であることの検証プロセスは、計算問題のサイズによって多項式的にのみ変化する実行時間を有する検証アルゴリズムによって実行可能である。換言すれば、証拠変数のセットは解に関する情報を含み、その情報は、解が「イエス」であることを検証するための検証アルゴリズムによって多項式実行時間で処理される。複雑性クラスＮＰの正式な定義に関しては、関連するコンピュータ科学の文献を参照して頂きたい。

例えば、巡回セールスマン問題、３彩色問題及びイジングスピンモデル問題は、複雑性クラスＮＰの決定問題の例である。例えば、イジングスピンモデル問題を検討してみる。結合係数及びフィールド係数の所定のセット、及び定数Ｋについてのイジングスピンモデル問題の解が「イエス」である場合、イジングエネルギー関数Ｈ（ｓ_１，ｓ_２，・・・，ｓ_ｎ）に関連するスピン（ｓ_１，ｓ_２，・・・，ｓ_ｎ）がＫよりも小さくなる形態は、証拠変数のセットと見なすことができる。証拠変数（ｓ_１，ｓ_２，・・・，ｓ_ｎ）が与えられると、エネルギーＨ（ｓ_１，ｓ_２，・・・，ｓ_ｎ）は、Ｈ（ｓ_１，ｓ_２，・・・，ｓ_ｎ）の数を計算しＫと比較することにより、確かにＫよりも小さいということが多項式時間で検証される。したがって、イジングスピンモデル問題は複雑性クラスＮＰに含まれる。

検証アルゴリズムが多項式実行時間を有するのに対して、決定問題で「イエス」又は「ノー」である解を決定するタスクは、ＮＰのいくつかの計算問題に対し多項式時間アルゴリズムを有していなくてもよく、指数関数的な実行時間を有しさえすればよい。複雑性クラスＮＰのいくつかの計算問題は、古典的計算システムではコンピュータ的に解決困難であると考えられる。「コンピュータ的に解決困難である」計算問題とは、計算問題の解が「イエス」であるか「ノー」であるかを決定するための多項式実行時間を有する古典的計算システムで実行されるアルゴリズムが存在しない計算問題である。具体的には、巡回セールスマン問題、３彩色問題及びイジングスピンモデル問題は、古典的計算システムでは解決困難であると考えられ、少なくともこれらの問題を多項式実行時間で解けるアルゴリズムは知られていない。

他の実施形態と組み合わせ可能な実施形態によれば、量子システムを用いて解を計算できる計算問題は、ＮＰ完全問題又はＮＰ困難問題である。ＮＰ完全問題は、クラスＮＰに属しており、古典的計算システムではコンピュータ的に解決困難であると考えられる。全てのＮＰ困難問題がＮＰに属するわけではないが、ＮＰ困難問題もまた古典的計算システムではコンピュータ的に解決困難であると考えられる。

ＮＰ完全問題のような複雑性クラスＮＰに属する計算問題についての実施形態によれば、測定の読み出しは計算問題の証拠変数のセット、又は少なくともこのセットの一部分を含んでもよい。

他の実施形態と組み合わせ可能な実施形態によれば、問題コード化設定を決定するステップは、予備の計算問題に計算問題をマッピングするステップを含んでもよい。この場合、予備の計算問題は、長距離スピンモデルの基底状態を決定するステップを含む。予備の計算問題は計算問題に依存する。計算問題を予備の計算問題にマッピングするステップは、計算問題の入力パラメータを予備の計算問題の入力パラメータにマッピングするステップを含んでもよい。計算問題を予備の計算問題にマッピングするステップは、計算問題の解が予備の計算問題の解から決定されるものであってもよい。

実施形態によれば、予備の計算問題は、イジングスピンモデル問題であってもよい。更なる実施形態によれば、計算問題は、例えば巡回セールスマン問題等の複雑性クラスＮＰの問題であってもよい。イジングスピンモデル問題はＮＰ完全問題であるため、巡回セールスマン問題等の複雑性クラスＮＰの全ての問題がイジングスピンモデル問題にマッピングされてもよい。例えば、第１リスト及び第２リストを含む巡回セールスマン問題の場合、第１リスト及び第２リストがイジングスピンモデル問題の結合係数及びフィールド係数のセットにマッピングされてもよい。巡回セールスマン問題の解は、対応する結合係数及びフィールド係数を有するイジングスピンモデル問題の解から計算されてもよい。このようなマッピングは既知である。

実施形態によれば、問題コード化設定を決定するステップは、イジングスピンモデル等の長距離スピンモデルから問題コード化設定を決定するステップを含んでもよい。図９〜図１６を参照して、この決定を実行する具体的な方法について詳細に説明する。

他の実施形態と組み合わせ可能な実施形態によれば、長距離スピンモデルは、ｍが１、２又は３のｍ体相互作用を有する長距離スピンモデルであってもよい。

他の実施形態と組み合わせ可能な実施形態によれば、前記方法は、長距離スピンモデルのスピンの複数の閉ループから短距離ハミルトニアンを決定するステップを含んでもよい。

更なる実施形態によれば、複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法が提供される。

前記方法は、計算問題を量子システムの問題ハミルトニアンにコード化するステップを含む。問題ハミルトニアンは、複数の調整可能なパラメータを含む単体ハミルトニアンである。上記のコード化するステップは、計算問題から、複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む。

前記方法は、初期量子状態の量子システムを初期化するステップを更に含む。前記方法は、量子アニーリングを実行することにより、量子システムを初期量子状態から最終量子状態に移行させるステップを更に含む。量子アニーリングを実行するステップは、量子システムの初期ハミルトニアンから量子システムの最終ハミルトニアンに移行するステップを含む。最終ハミルトニアンは、問題ハミルトニアン及び短距離ハミルトニアンの総和である。この場合、問題ハミルトニアンの複数の調整可能なパラメータは問題コード化設定に存在し、短距離ハミルトニアンはｄ体ハミルトニアンであり、ｄは計算問題から独立している。

前記方法は、複数のキュービットの少なくとも一部分を測定することで最終量子状態の読み出しを取得するステップを更に含む。

前記方法は、読み出しから計算問題の解を決定するステップを更に含む。

更なる実施形態によれば、図１に示す装置４００等のような計算問題の解を計算する装置が提供される。

前記装置は、複数のキュービットを備えた量子システムを有する。

前記装置は、図１に示す冷却ユニット４１０等のような、量子システムをその基底状態に向けて冷却するのに適した冷却ユニットを更に有する。冷却ユニットは、初期量子状態の量子システムを初期化するために、量子システムを初期ハミルトニアンの基底状態に向けて冷却するように構成されてもよい。冷却ユニットは、量子システムが、前記装置に使用されるキュービットの種類に強く依存する動作温度に維持されるように構成されてもよい。例えば、超電導キュービットの場合、動作温度は５０ｍＫ以下で、特に１ｍＫ以下である。

前記装置は、図１に示すプログラム可能な量子アニーリングユニット４３０等のような、量子アニーリングによって量子システムの初期ハミルトニアンを量子システムの最終ハミルトニアンに移行させるのに適したプログラム可能な量子アニーリングユニットを更に有する。最終ハミルトニアンは、問題ハミルトニアン及び短距離ハミルトニアンの総和である。

冷却ユニットは、プログラム可能な量子アニーリングユニットにより、量子アニーリングを行っている間、量子システムを動作温度に維持するように構成されてもよい。

前記装置は、図１に示す測定装置４４０等のような、複数のキュービットの少なくとも一部分を測定するのに適した測定装置を更に有する。

前記装置は、図１に示す古典的計算システム４５０等のような、プログラム可能な量子アニーリングユニット及び測定装置に接続された古典的計算システムを更に有する。古典的計算システムは、計算問題を入力として受信するように構成されてもよい。古典的計算システムは、計算問題を問題ハミルトニアンにコード化するように更に構成されてもよい。この場合、コード化は、問題ハミルトニアンの複数の調整可能なパラメータに対する問題コード化設定を、計算問題から決定することを含んでもよい。古典的計算システムは、問題コード化設定を量子アニーリングユニットに伝達するように更に構成されてもよい。

プログラム可能な量子アニーリングユニットは、古典的計算システムから問題コード化設定を受信するように構成されてもよい。プログラム可能な量子アニーリングユニットは、量子アニーリングによって、初期ハミルトニアンから最終ハミルトニアンに移行するように構成されてもよい。この場合、問題ハミルトニアンの複数の調整可能なパラメータは問題コード化設定に存在する。

古典的計算システムは、測定装置から量子システムの読み出しを受信するように更に構成されてもよい。古典的計算システムは、読み出しから計算問題の解を決定するように更に構成されてもよい。

更なる実施形態によれば、計算問題の解を計算するためのプログラム可能な量子アニーリング装置が提供される。プログラム可能な量子アニーリング装置は、２次元格子に沿って配置された複数の超電導キュービットを含む量子システムを有する。

プログラム可能な量子アニーリング装置は、複数の調整可能な磁束を生成するように構成された複数の磁束バイアスユニットを含む磁束バイアスアセンブリを更に有する。各調整可能な磁束は、複数の超電導キュービットのうちの単一の超電導キュービットに作用する。

プログラム可能な量子アニーリング装置は、プラケットハミルトニアンに従い複数の超電導キュービットを結合するように構成された少なくとも１つの超電導量子干渉装置を含む結合ユニットを更に有する。

プログラム可能な量子アニーリング装置は、磁束バイアスユニット及び結合ユニットに接続された制御装置を更に有する。制御装置は、量子システムの問題ハミルトニアンの複数の調整可能なパラメータに対する問題コード化設定を受信するように構成される。この場合、問題ハミルトニアンは単体ハミルトニアンであり、問題コード化設定は計算問題をコード化する。制御装置はさらに、量子アニーリングにより量子システムの初期ハミルトニアンを量子システムの最終ハミルトニアンに展開するために、磁束バイアスアセンブリ及び結合ユニットを制御するよう構成される。最終ハミルトニアンは、プラケットハミルトニアン及び問題ハミルトニアンの合計である。その場合、問題ハミルトニアンの複数の調整可能なパラメータは、問題コード化設定に存在する。

本明細書において、「プログラム可能な量子アニーリングユニット」は、「プログラム可能な量子アニーリング装置」の同意語として使用される。

図９〜図１６を参照して、本開示事項の更なる側面について説明する。計算問題から問題ハミルトニアンへ、そして対応する最終ハミルトニアンへの具体的なコード化について説明する。すなわち、長距離相互作用を有するイジングスピンモデル問題から、単体の問題ハミルトニアン及びプラケットハミルトニアンの総和である最終量子ハミルトニアンへのコード化である。長距離相互作用を有する（古典的）イジングスピンモデル問題は、ＮＰ完全であり、その量子化は些細なものなので、ここでは古典的な量子イジングスピンとそうでないものとを区別しない。他の古典的計算問題をイジングスピンモデル問題へマッピングすることは公知である。したがって、最終量子ハミルトニアンの基底状態、又は低い動作温度での熱状態は、イジングスピンモデル問題の解、及び逆マッピングを実行することで多くの古典的でＮＰ困難な計算問題の解に関する情報を含むことができる。最終量子ハミルトニアンへの具体的なマッピングは、イジングスピンモデル問題がｄが２以下のｄ体相互作用のみを含む場合、２次元表面、特に２次元格子での量子処理装置（量子プロセッサ）の実現を可能にし、イジングスピンモデル問題がｄが３以下のｄ体相互作用のみを含む場合、３次元空間、特に３次元格子での量子処理装置の実現を可能にする。マッピングは、ｄ体相互作用および任意のｄを有するイジングスピンモデル問題に拡大適用することができる。量子処理装置は、単体の問題ハミルトニアンを通じて完全にプログラム可能であり、拡張性のあるアーキテクチャである。

最初に、ｄが２以下のｄ体相互作用のみを含むイジングスピンモデル問題について考える。具体的なコード化は、最大２体相互作用とこれに対応する結合係数ｃ_ｉｊを有する、ｎ個のスピンに対するイジングスピンモデル問題から開始される。指数ｉ及びｊは１〜ｎの範囲で、ｊはｉよりも小さい。この最初のケースでは、全てのフィールド係数ｃ_ｉは零に等しい。図９は、ｎ＝６のスピンに対するイジングスピンモデル問題を示す。図中、スピンには１〜６の符号が付されている。図９においてスピン対を結ぶ線で示すように、スピン間にはｎ（ｎ−１）／２＝１５対の対相互作用がある。例えば、線１２はスピン１及び２の間の対相互作用を示す。１５対の対相互作用は１５個の結合係数ｃ_ｉｊに対応する。相互作用は長距離相互作用である。

イジングスピンモデルにおける全てのスピン対に対し、対応する量子システムのキュービットが与えられる。例えば、図９に示す１５対の対相互作用を有する６つのスピンの場合、対応する量子システムは１５個のキュービットを有する。イジングスピンモデルのスピンの設定は、対応するキュービットの設定にマッピングされる。この際、キュービットの設定は、スピンの相対的な方向に依存する。同じ方向を向く一対のスピン（平行配列）は、量子基本状態「｜１＞」のキュービットにマッピングされる。反対方向を向く一対のスピン（逆平行配列）は、量子基本状態「｜０＞」のキュービットにマッピングされる。このマッピングを図１０に示す。図１０で、符号０及び１はそれぞれ、量子基本状態の｜０＞及び｜１＞に対応する。

結合係数ｃ_ｉｊは、計算問題（この場合は、イジングスピンモデル問題）をコード化する問題ハミルトニアンの複数の調整可能なパラメータＪ_ｋにマッピングされる。問題ハミルトニアンは、式Σ_ｋＪ_ｋσ_ｚ ^（ｋ）を有する。ｋ＝ｎ^＊ｉ＋ｊであり、ｋはＭ＝ｎ（ｎ−１）／２である１〜Ｍの範囲である。イジングスピンモデル問題は、問題ハミルトニアンの調整可能なパラメータＪ_ｋが、結合係数ｃ_ｉｊに対応するイジングスピンモデルのスピン間の相互作用を表すように、問題ハミルトニアンにマッピングされる。

イジングスピンモデル問題を問題ハミルトニアンにコード化するのに必要なキュービットの数は、ｎ個のスピンに対するイジングスピンモデル問題に比べて二次的に増大する。なぜなら、スピン間の２体相互作用の数がＭ＝ｎ（ｎ−１）／２と等しいからである。いくつかの実施形態によれば、追加の自由度が考慮される。量子システムにおけるキュービットの総数は、Ｍ＋ｎ−２以上であってもよい。この場合、ｎ−２個の追加の補助キュービット及び／又は追加の予備キュービットが後述の理由で追加される。したがって、キュービットの数はスピンの数ｎより大きい。具体的には、キュービットの数はスピンの数ｎにＭ−２の追加の自由度を付加した数でもよい。問題ハミルトニアンは、ローカル相互作用のみで、具体的には外部フィールドとの単体相互作用で、量子処理装置のプログラミングを可能にする。

イジングスピンモデルに比べて、量子システムの増加した自由度は、Ｍ−ｎ個の４体被加数ハミルトニアンＣ_ｌの総和である短距離ハミルトニアンにより相殺される。これを制約ハミルトニアンと称し、キュービットの一部分を固定するための制約を表す。短距離ハミルトニアンは式Σ_ｌＣ_ｌを有する。指数ｌの範囲は１〜（ｎ^２−３ｎ）／２であり、各被加数ハミルトニアンＣ_ｌは下記の式を有する制約ハミルトニアンである。

上記式を参照して、制約ハミルトニアンの２通りの可能な実行について考える。上記式の総和は、補助に基づく実行を表す。この総和は、キュービットが配置される２次元格子のプラケットの４つの要素（北、東、南、西）を実行する。さらに、各Ｓ_ｚ ^ｌは量子システムに含まれる補助キュートリットに作用する演算子である。補助キュートリットには、３つの基本状態からなる基準があり、本実施形態では、｜０＞、｜２＞及び｜４＞と表す。短距離ハミルトニアンの２つ目の実行は、補助キュートリットを必要としない相互作用に基づく実行である。相互作用に基づく実行では、Ｃ_ｌは格子のプラケットを形成するキュービット間の４体相互作用である。さらに、上記式中、Ｃは一定の制約強度等の制約強度を表す。

上述のように、問題ハミルトニアンでのイジングスピンモデルのコード化は、イジングスピンモデルのスピンの設定を量子システムのキュービットの設定にマッピングすることを伴う。この場合、キュービットの設定は、対応するスピンの設定のスピン対の相対的な方向に依存する。一貫性のあるマッピングを提供するために、後述のようなイジングスピンモデルにおける閉ループに関わる側面が検討される。イジングスピンモデルにおけるスピンの各閉ループで、逆平行配列を有するスピン対の数は偶数である。例えば、図９を参照して、破線で示される連結部１４、２４、２３及び１３で形成された閉ループ等について考える。この閉ループは、スピン１、２、３及び４を含む。スピン１、２、３及び４のいずれの設定も零対、２対又は４対の逆平行スピンを含む。スピン１、２、３及び４のいずれの設定も１対又は３対の逆平行スピンを含まない。したがって、スピン１、２、３及び４の全ての設定は、偶数の逆平行スピンを有する。

逆平行のスピン対が量子基本状態｜０＞にあるキュービットにマッピングされるため、イジングスピンモデルのスピンの閉ループに対応する量子システムのキュービットの全てのセットは、偶数個の量子基本状態｜０＞を有する。これにより、量子システムのキュービットの少なくとも一部分に対する制約のセットが提供される。例えば、図９に関する上記の閉ループの場合、対応する４つのキュービットのグループを、イジングモデルのスピン対と量子システムのキュービット間の対応関係を考慮して、参照符号１４、２４、２３及び１３を付して図１１に示す。図１１に示すように、キュービット１４、２４、２３及び１３は、２次元格子１２０のプラケットに対応する。前述のような閉ループの制約を考慮すると、キュービット１４、２４、２３及び１３に対する量子基本状態のいずれの設定も、図１２に示すように、零、２つ又は４つのいずれかの量子基本状態｜０＞を含む。

全ての閉ループに対応する制約が確実に満たされるには、適切な閉ループのサブセットに関連する制約を実行することで十分である。本実施形態によれば、最大４つのスピンからなるグループを含む閉ループの特定の構成要素（building blocks of closed loops）が、全ての制約が確実に満たされるのに十分であるため、イジングスピンモデルから量子システムへ一貫性のあるマッピングが提供される。前記構成要素は、４つの連結部によって連結された４つのスピンからなる閉ループを伴う。１つの連結部は指標距離（index distance）ｓを有し、２つの連結部は指標距離ｓ＋１を有し、もう１つの連結部は指標距離ｓ＋２を有する。ｓは１〜Ｎ−２の範囲であり、スピンｓ_ｉ及びｓ_ｊ間の「指標距離」は、｜ｉ−ｊ｜である。ｓ＝１の閉ループ構成要素のセットは、ｎ−２個の制約を提供する。例えば、図９に示すようにスピン１、２、３及び４の連結部１４、２４、２３及び１３を含む閉ループが、ｓ＝１の閉ループ構成要素である。

更なる側面は、量子システムの境界に関する。いくつかの閉ループ構成要素は、４つの連結部で連結された４つのスピンではなく、３つの連結部で連結された３つのスピンからなるグループを伴う。この点に関し、例えば図９を参照すれば、スピン１、２及び３の連結部１２、２３及び１３を含む閉ループが考えられる。量子システムにおける対応するキュービットのグループは、２次元格子の３角形状のプラケットに沿って配置された３つのキュービット１２、２３及び１３を含む。３つのスピンの閉ループに対応する制約を実行するために、３体制約ハミルトニアンＣ_ｌが、対応する３つのキュービットからなるグループに作用すると考えられる。或いは、図１１に破線の円で示すような、量子基本状態｜１＞で固定されたｎ−２個の予備キュービットの追加の列が量子システムに含まれていてもよい。例えばキュービット１２、２３及び１３に対応する閉ループ等のように、３つのスピンの閉ループに対応する制約を実行するために、制約ハミルトニアンＣ_ｌが対応する３つのキュービットと、予備キュービットの１つ（すなわち、図１１に示す予備キュービット１１０１）に作用すると考えられる。したがって、制約ハミルトニアンＣ_ｌは、前述と同じ式を有する拡大された２次元格子のプラケットに作用する４体ハミルトニアンである。全ての制約ハミルトニアンが２次元格子のプラケットに対応する４体ハミルトニアンであるため、後者の実現には、全ての制約ハミルトニアンを同じ条件で処置できるという利点がある。

制約ハミルトニアンＣ_ｌは、閉ループ構成要素に対応する制約、ひいては全ての閉ループに対応する制約が満たされることを確実にする。したがって、短距離ハミルトニアンは、イジングスピンモデルにおけるスピンの制約から量子システムに課される制約への一貫性のあるマッピングを提供する。

読み出しを与えるために、図１１で示す一部分４２５等のようにキュービットの一部分を測定可能である。量子システムが最終ハミルトニアンの基底状態にある場合、一部分４２５のキュービットは、イジングスピンモデルの基底状態にあるスピンの設定に対応する量子基本状態の設定に存在する。量子システムが基底状態に近い最終ハミルトニアンの熱状態、つまり十分に低い温度にある場合、これが高確率に当てはまる。したがって、一部分４２５を測定することにより、少なくとも高確率でイジングスピンモデル問題の解を決定することが可能になる。量子システムが最終ハミルトニアンの基底状態にうまく近似された最終状態にある場合、一部分４２５を測定することで、試行解を計算できるイジングモデルの基底状態に関する情報が少なくとも提供される。そして、その試行解が実際の解であるかどうか、多項式時間の古典的計算で分析できる。もし試行解が実際の解でなければ、実際の解が見つかるまで計算を繰り返すことができる。

実施形態の更なる利点として、イジングスピンモデルに関する情報は量子システムで冗長的にコード化されるため、計算問題の解が決定される読み出しを提供するためにキュービットの生じ得る様々なグループを測定可能である。

上記の観点から、本実施形態に係る短距離ハミルトニアンの構造は、（i）制約がスピン間の全ての相互作用を網羅し、（ii）制約の数が（ｎ^２−３ｎ）／２で、（iii）短距離ハミルトニアンが、ｄ体相互作用を有する単純な２次元形状で実現可能である。この場合、ｄ＝４であり、相互作用は２次元格子のプラケットに対応する。さらに、イジングスピンモデルに１つのスピンを追加することは量子システムにｎ個のスピンの一列を追加するのと同等であるため、本実施形態は拡張性のある実行を可能にする。

図９〜図１２を参照して説明する実施形態は、ｎ個のスピン間の対相互作用を伴うイジングスピンモデルに関する。この場合のフィールド係数は零である。零でないフィールド係数を有するイジングスピンモデルに対して同様のコード化が考えられる。＋１の値に固定された追加のスピンｓ_ｎ＋１をイジングモデルに含めることができる。そして、零でないフィールド係数は、ｎ個のスピンと追加のスピンｓ_ｎ＋１との間の結合係数として再定式化される。零でないフィールド係数を有するイジングスピンモデルは、フィールド係数が零のイジングスピンモデルにマッピングされる。したがって、上述の方法によって量子システムへのマッピングが適用される。追加スピンｓ_ｎ＋１を追加することは、量子システムに追加のｎ個のキュービットの一列を含めることを伴う。

更に、３つのスピンからなるグループ間の相互作用を伴うイジングスピンモデルに対するコード化が考えられる。この場合、イジングエネルギー関数は下記の式を有していてもよい。
Ｈ（ｓ_１，ｓ_２，・・・，ｓ_ｎ）＝Σ_ｉｊｋｃ_ｉｊｋｓ_ｉｓ_ｊｓ_ｋ
ここで、係数ｃ_ｉｊｋは、スピンｓ_ｉ、ｓ_ｊ及びｓ_ｋ間の３体相互作用を表し、ｉ＞ｊ＞ｋである。このような３体イジングモデルを量子システムにマッピングするステップ、及び対応するイジングスピンモデル問題を量子システムの問題ハミルトニアンへコード化するステップを図１３〜図１６に示す。本実施形態では、量子システムのキュービットは、３体イジングモデルの３つのスピンの組に対応する。３体イジングモデルでは、Ｒ＝ｎ（ｎ−１）（ｎ−２）／６組の３つのスピンが存在する。したがって、キュービットの数はＲ以上である。この場合、前述の２体イジングスピンモデルのマッピングと同様に、補助キュービット及び／又は予備キュービット等の追加のキュービットが含まれてもよい。本実施形態では、図１６に示すように、複数のキュービットが３次元正方格子１６０１に沿って配置される。問題ハミルトニアンは、スピン間で最大２体相互作用を伴うイジングスピンモデルの場合と同様に、式Σ_ｋＪ_ｋσ_ｚ ^（ｋ）を有してもよい。短距離ハミルトニアンは、式Σ_ｌＣ_ｌを有していてもよい。この場合、制約ハミルトニアンＣ_ｌは３次元正方格子のプラケットに対応する。制約ハミルトニアンの数は、２（Ｒ−ｎ）個であってもよい。３体制約ハミルトニアンを伴う、及び／又は、補助キュービット及び／又は予備キュービットの含有を伴う量子システムの境界に関する同様の考えが本実施形態にも当てはまる。

実施形態の更なる利点は、量子システムにおける誤差の防止及び／又は訂正である。例えば、４つのキュービットを有する小型の量子システムでは、Ｊ＝ｍａｘ（｜ｃ_ｉｊ｜）であるＣが１．５Ｊよりも大きい場合、特にＣが約１．５Ｊである場合、静誤差が減少し、零になることもある。一般的に、Ｃが（ｎ−２）ｍａｘ（｜ｃ_ｉｊ｜）よりも大きい場合、静誤差は減少、又はゼロになり得る。上記のｎはイジングスピンモデルのスピンの数を表す。更に、量子システムにおけるキュービットの数Ｎは、例えばＮがｎ（ｎ−１）／２以上の場合のように、イジングスピンモデルにおけるスピンの数よりも大きくなる。したがって、イジングスピンモデルに関する情報は、量子システムにおいて冗長的にコード化される。この冗長さによって、デコヒーレンスによる誤差に対し、低感度の誤差訂正測定の実行が可能になる。トポロジカル量子メモリの誤差訂正スキームと同様に、キュービットの多様な生じ得る組み合わせを測定することで、量子システムのデコヒーレンスを検知して訂正することが可能になる。図１７は、量子システムにおける誤差の防止及び／又は訂正に関する実施形態の利点を示す。図１７では、軸１７０１はスピンの数ｎに関する。軸１７０２は誤差の大きさに関する。更に、曲線１７１０、１７２０及び１７３０はそれぞれ、誤差の変動、情報損失、及びスピンの反転の数を示す。挿入図の軸１７０３はスピンの数に関し、軸１７０３は読み出しの数に関する。

更なる実施形態によれば、量子処理ユニット［ＱＰＵ］が提供される。量子処理ユニットは、計算問題の解を計算する装置の中央処理装置として作動するのに適している。量子処理ユニット［ＱＰＵ］は、量子システムを有する。量子システム［ＱＳ］は、複数のキュービットを有する。量子処理ユニットは、実施形態に係る方法を実行するのに適している。

量子処理ユニットは、初期時間ｔ＝ｔ_０での初期ハミルトニアンＨ^ｉｎｉｔを統合するのに適していてもよい。量子処理ユニットは、最終時間ｔ＝ｔ_ｆｉｎでの最終ハミルトニアンＨ^{ｆｉｎａｌ}を統合するのに更に適していてもよい。ここで、「ハミルトニアンを統合する」とは、本発明の量子計算機能に関連する量子システム［ＱＳ］の量子状態が、量子物理学の枠組みに従った各ハミルトニアンで主に決定されることを意味する。つまり、前記ハミルトニアンによって量子システム［ＱＳ］の量子計算特性が十分に判定され、量子システム［ＱＳ］の実際のハミルトニアンで考えられる更なる数式の項は、本発明の機能的な原理に関して無視できるほどである。

最終ハミルトニアンＨ^{ｆｉｎａｌ}は、問題ハミルトニアンＨ^ｐｒｏｂ及び短距離ハミルトニアンＨ^ＳＲの総和、つまり、Ｈ^{ｆｉｎａｌ}＝Ｈ^ｐｒｏｂ＋Ｈ^ＳＲであってもよい。この場合、問題ハミルトニアンＨ^ｐｒｏｂは、複数の調整可能なパラメータＪ_ｋを有する、つまりはその関数である単体ハミルトニアンであってもよい。調整可能なパラメータＪ_ｋは、計算問題のコード化に適した方法で個別に調整可能であってもよい。

短距離ハミルトニアンＨ^ＳＲは、ｄが２以上のｄ体ハミルトニアンであってもよく、ｄは計算問題から独立していてもよい。この場合、「短距離ハミルトニアン」という語句は、複数のキュービットの相互作用を表すハミルトニアンを意味してもよい。この場合、相互作用遮断距離Ｌ_ｃｕｔよりも長い距離で互いに離れたキュービット間で相互作用は生じない。「ｄ体ハミルトニアン」という語句は、複数のキュービットの相互作用を表すハミルトニアンを意味してもよい。この場合、ｄ＋１個以上のキュービットからなるグループ間で共有の相互作用（joint interactions）は生じない。

量子処理ユニット［ＱＰＵ］は、時間ｔ＝ｔ_ｆｉｎ後の、量子システム［ＱＳ］の複数のキュービットの少なくとも一部分の量子状態を測定するのに適していてもよい。

量子処理ユニット［ＱＰＵ］は、式Ｈ（ｔ）＝Ａ（ｔ）・Ｈ^０＋Ｂ（ｔ）・Ｈ^ｐｒｏｂ＋Ｃ（ｔ）・Ｈ^ＳＲの補間ハミルトニアンＨ（ｔ）を統合するのに適していてもよい。式中、Ａ（ｔ）、Ｂ（ｔ）及びＣ（ｔ）は時間パラメータｔに依存する補間係数であり、Ｈ^０は所望の計算プロセスに対し量子システムを初期化するのに適した初期化ハミルトニアンである。この場合、Ａ（ｔ）、Ｂ（ｔ）及びＣ（ｔ）は、Ｈ（ｔ_０）＝Ｈ^ｉｎｉｔ＝Ｈ^０＋Ｃ（ｔ_０）・Ｈ^ＳＲ及びＨ（ｔ_ｆｉｎ）＝Ｈ^ｆｉｎ＝Ｈ^ｐｒｏｂ＋Ｈ^ＳＲとなるように、Ａ（ｔ_０）＝１、Ａ（ｔ_ｆｉｎ）＝０、Ｂ（ｔ_０）＝０、Ｂ（ｔ_ｆｉｎ）＝１、Ｃ（ｔ_ｆｉｎ）＝１の条件を満たす。Ｃ（ｔ_０）は任意であり、０又は１であってもよい。

複数のキュービットにおけるキュービットは、２次元格子又は３次元格子に沿って配置されてもよい。

複数のキュービットにおけるキュービットは、各格子のプラケットの角に位置していてもよく、「プラケット」という語句は各格子のメッシュを意味してもよい。

問題ハミルトニアンは、式Ｈ^ｐｒｏｂ＝Σ_ｋ［Ｊ_ｋ・σ_ｚ ^（ｋ）］を有してもよい。式中、σ_ｚ ^（ｋ）は第１空間方向ｚ（ｒ_ｋ）に関連し、複数のキュービットのｋ番目のキュービットに作用するパウリ演算子であってもよい。ｚは前記第１空間方向の方向ベクトルを意味してもよく、ｒ_ｋはｋ番目のキュービットの位置ベクトルを示す。

初期化ハミルトニアンＨ^０は、所望の計算プロセスに対し量子システムを初期化するのに適し、且つ好ましくはそれを容易に実施できる式を有してもよい。例えば、初期化ハミルトニアンは、式Ｈ^０＝ｈΣ_ｋ［σ_ｘ ^（ｋ）］を有してもよい。式中、ｈは係数で、σ_ｘ ^（ｋ）は第２空間方向ｘ（ｒ_ｋ）に関連し、複数のキュービットのｋ番目のキュービットに作用するパウリ演算子で、ｘは前記第２空間方向の方向ベクトルを示し、ｒ_ｋはｋ番目のキュービットの位置ベクトルを示す。この場合、ｚ（ｒ_ｋ）及びｘ（ｒ_ｋ）は互いに直交していることが好ましく、パウリ演算子σ_ｚ ^（ｋ）及びσ_ｘ ^（ｋ）は非可換の演算子、特に反可換の演算子であることが好ましい。

短距離ハミルトニアンＨ^ＳＲは、式Ｈ^ＳＲ＝Ｈ^Ｐ：＝Σ_ｌ［Ｃ_ｌ（σ_ｚ ^{（ｌ，１）}，・・・，σ_ｚ ^{（ｌ，ｍ）}，・・・，σ_ｚ ^{（ｌ，Ｍ［ｌ］）}）］を有するプラケットハミルトニアンＨ^Ｐであってもよい。式中、ｌはｌ番目のプラケットの数を示し、（ｌ，１）はｌ番目のプラケットの第１の角に位置するキュービットを示し、（ｌ，ｍ）はｌ番目のプラケットのｍ番目の角に位置するキュービットを示し、（ｌ，Ｍ［ｌ］）はｌ番目のプラケットの最後（つまり、Ｍ［ｌ］番目）の角に位置するキュービットを示す。式中、σ_ｚ ^{（ｌ，ｍ）}は各キュービットのパウリ演算子を示し、各Ｍ［ｌ］体被加数Ｃ_ｌ（σ_ｚ ^{（ｌ，１）}，・・・，σ_ｚ ^{（ｌ，ｍ）}，・・・，σ_ｚ ^{（ｌ，Ｍ［ｌ］）}）はプラケットハミルトニアンＨ^Ｐへのｌ番目のプラケットの寄与度を示す。

短距離ハミルトニアンＨ^ＳＲは、各キュービット−格子のプラケットを形成するキュービットのグループのみが被加数をＨ^ＳＲ＝Ｈ^Ｐに寄与させるような方法で実行されてもよい。そして、例えば、三角形のプラケットのみからなるキュービット−格子の場合、Ｈ^ＳＲ＝Ｈ^Ｐはｄ＝Ｍ＝３であるｄ体ハミルトニアンであり、三角形及び四角形のプラケットからなるキュービット−格子の場合、Ｈ^Ｐは３体及び４体の被加数からなるハミルトニアンであり、したがってそれは４体ハミルトニアンである。

量子処理ユニット［ＱＰＵ］は、量子システムの所望の特徴、特に短距離ハミルトニアンＨ^ＳＲ及びプラケットハミルトニアンのそれぞれの所望の特徴を実行するのに適した追加のキュービット及び／又はキュートリット及び／又はｑ段階量子システム及び／又は更なる構成要素及び装置を含んでもよい。

キュービットが配置される格子は、実質的に平面的な２次元格子、又は実質的に互いに平行で３次元的に積み重ねられた実質的に平面的な２次元副格子からなる３次元格子であってもよい。

キュービットは、それぞれ２次元格子の実質的に三角形の部分に配置されてもよいし、各２次元副格子の基本的に三角形の部分に配置されてもよい。

各キュービットの配置により形成された複数のプラケットは、（多数の）四角形のプラケットのセット、及び（少数の）三角形のプラケットのセットから構成されてもよい。或いは又は更に、各キュービットの配置により形成された複数のプラケットは、例えば、全てのプラケットを四角形にするために量子システムに更なるキュービットを追加することにより、四角形のプラケットのみから構成されてもよい。

プラケットハミルトニアンＨ^Ｐ＝Σ_ｌ［Ｃ_ｌ］の被加数Ｃ_ｌは、下記２つの式のうちの１つであってもよい。

量子システムのキュービットの空間的配置、及び／又は問題ハミルトニアンＨ^ｐｒｏｂの係数Ｊ_ｋを量子システムのキュービットの空間的座標へ割り当てること、及び／又は短距離ハミルトニアンの実行、及び／又は読み出しを提供するために測定されるキュービットの一部分の有利な決定は、公知の２体量子相互作用モデル、又は公知の３体量子相互作用モデル（例えば、２体又は３体の全対全イジングスピンモデル等）を量子システム［ＱＳ］にマッピングするステップに従い実行可能である。特に、前記のマッピングから推論できる制約（例えば、閉ループの制約等）に従い実行可能である。

量子処理ユニットは、量子システム［ＱＳ］をＨ^ｉｎｉｔ：＝Ｈ（ｔ_０）の基底状態に向けて冷却するのに適していてもよい。

量子処理ユニットは、量子システム［ＱＳ］をＨ^ｆｉｎ：＝Ｈ（ｔ_ｆｉｎ）の基底状態に向けて冷却するのに適していてもよい。

量子処理ユニットは、ｔ＝ｔ_０とｔ＝ｔ_ｆｉｎの間の任意の時間において又は全ての時間において、量子システム［ＱＳ］をＨ（ｔ）の基底状態に向けて冷却するのに適していてもよい。

量子処理ユニットは、補間ハミルトニアンＨ（ｔ）をＨ^ｉｎｉｔからＨ^ｆｉｎに徐々に移行させる、特にその移行を断熱的量子アニーリングプロトコルに従い実行するのに適していてもよい。

量子処理ユニットは、２次元格子に沿って配置された複数の超電導キュービットを有してもよい。

量子処理ユニットは、複数の調整可能な磁束を生成するように構成された複数の磁束バイアスユニットを含む磁束バイアスアセンブリを有してもよい。この場合、各調整可能な磁束は複数の超電導キュービットのうちの単一の超電導キュービットに作用する。

量子処理ユニットは、プラケットハミルトニアンＨ^Ｐに従い複数の超電導キュービットを結合するように構成された少なくとも１つの超電導量子干渉装置を含む結合ユニットを有してもよい。

量子処理ユニットは、磁束バイアスユニット及び結合ユニットに接続された制御装置を有してもよい。

制御装置は、量子システムの問題ハミルトニアンＨ^Ｐの複数の調整可能なパラメータに対して問題コード化設定を受信するように構成されてもよい。この場合、問題ハミルトニアンは単体ハミルトニアンであり、問題コード化設定は計算問題をコード化する。

制御装置は、磁束バイアスアセンブリ及び結合ユニットを制御して、量子アニーリングによって量子システムの初期ハミルトニアンＨ^ｉｎｉｔを量子システムの最終ハミルトニアンＨ^ｆｉｎに移行させるように構成されてもよい。この場合、最終ハミルトニアンＨ^ｆｉｎは、プラケットハミルトニアンＨ^Ｐ及び問題ハミルトニアンＨ^ｐｒｏｂの総和であり、問題ハミルトニアンＨ^ｐｒｏｂの複数の調整可能なパラメータＪ_ｋは問題コード化設定に存在する。

実施形態によれば、前記装置は、実施形態に係る方法を実施するのに適している。

更なる実施形態によれば、量子処理装置が提供される。量子処理装置は、計算問題の解を計算するように構成されていてもよい。量子処理装置は、量子処理ユニット（ＱＰＵ）を有する。量子処理ユニットは、キュービットのプラケットを有する。プラケットは、２次元又は３次元格子の基本セルであってもよい。プラケットは、最大４つ又は最大３つのキュービットを含んでもよい。量子処理ユニットは、キュービット及び１つ以上の補助キュービット（例えば、特定の量子状態で準備された１つの補助キュービット等）を含むプラケットを有してもよい。キュービット又は場合によってはキュービット及び補助キュービットは、プラケットの角に配置されてもよい。量子処理装置は、プラケットのいくつか又は全てに配置された補助キュートリット等の補助ｑ段階システムを有してもよい。補助ｑ段階システムは、プラケットのキュービット間又はキュービット及び補助キュービット間の相互作用を仲介してもよい。プラケットは、正方格子又はその一部を成すように配置されてもよい。特に、プラケットは、正方格子の三角形部分を成すように配置されてもよい。

量子処理装置は、単一のキュービット又は補助キュービットと相互作用するユニットを有する第１装置、及びプラケットのキュービット又はプラケットのキュービット及び補助キュービットと相互作用するユニットを有する第２装置を更に有してもよい。第１装置は、プラケットのキュービット又はプラケットのキュービット及び補助キュービットに単体ハミルトニアンを実行するように構成されてもよい。第２装置は、プラケットのキュービット又はプラケットのキュービット及び補助キュービットにプラケットハミルトニアンを実行するように構成されてもよい。

いくつかの実施形態によれば、キュービットは超電導キュービットである。補助キュービットが存在するのであれば、補助キュービットも超電導キュービットであってもよい。プラケットは、２次元正方格子又はその一部を形成するように、特にその三角形部分を形成するように配置される。量子処理ユニットは、複数の調整可能な磁束を生成するように構成された複数の磁束バイアスユニットを備えた磁束バイアスアセンブリを有する。この場合、各調整可能な磁束は、プラケットの単一の超電導キュービットに作用する。磁束バイアスアセンブリは、２次元又は３次元格子のプラケットの超電導キュービット又は超電導キュービット及び超電導補助キュービットに単体ハミルトニアンを実行するように構成されてもよい。量子処理ユニットは、超電導量子干渉装置を有する。この場合、各超電導量子干渉装置はプラケットの超電導キュービット又はプラケットの超電導キュービット及び超電導補助キュービットを結合する。量子干渉装置は、プラケットの超電導キュービット又はプラケットの超電導キュービット及び超電導補助キュービットの制御された相互作用を提供する。量子干渉装置は、２次元又は３次元格子のプラケットにプラケットハミルトニアンを実行するように構成された量子結合ユニットの構成要素であってもよい。

量子処理ユニット（ＱＰＵ）は、計算問題の解を計算する方法を実行するように構成されてもよい。量子処理装置は、冷却ユニット、古典的計算装置、及び制御装置等の構成要素を有してもよい。例えば、制御装置は、磁束バイアスユニット及び結合ユニットに接続されてもよく、（i）量子システムの問題ハミルトニアンの複数の調整可能なパラメータに対する問題コード化設定を受信するように構成されてもよい。この場合、問題ハミルトニアンは単体ハミルトニアンであり、問題コード化設定は計算問題をコード化する。そして、（ii）量子アニーリングによって量子システムの初期ハミルトニアンを量子システムの最終ハミルトニアンに移行させるために磁束バイアスアセンブリ及び結合ユニットを制御するように構成されてもよい。この場合、最終ハミルトニアンはプラケットハミルトニアン及び問題ハミルトニアンの総和であり、問題ハミルトニアンの複数の調整可能なパラメータは問題コード化設定に存在する。

実施形態は、キュービット（つまり、２段階量子ビット）のみでなく、キュートリットやｑが任意であるｑ段階量子システムで実現されてもよい。超電導キュービット、トラップイオン、量子ドット、及びＮＶセンター等を用いて量子システムが実現される別の態様は、２つ以上の状態を考慮することによって、ｑ段階システム、特にキュートリットシステムに拡張可能である。ｑ段階システムの操作及び読み出し、及びこれらの相互作用は、上述の実施形態を拡張させることにより実現できる。

上述の内容は本発明のいくつかの実施形態を対象にしているが、それ以外の更なる実施形態についても、特許請求の範囲によって決まる範囲から逸脱することなく創出されてもよい。

Claims (15)

1. 複数のキュービットを有する量子システムを用いて計算問題の解を計算する方法であって、
   前記計算問題を量子システムの問題ハミルトニアンにコード化するステップであって、前記問題ハミルトニアンが、複数の調整可能なパラメータを有する単体ハミルトニアンであり、前記計算問題から前記複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む、前記コード化するステップと、
   前記量子システムを初期量子状態から前記量子システムの最終ハミルトニアンの基底状態に向けて移行させるステップであって、前記最終ハミルトニアンは前記問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンの前記複数の調整可能なパラメータは前記問題コード化設定に存在し、前記短距離ハミルトニアンはｄが前記計算問題から独立しているｄ体ハミルトニアンである、前記移行させるステップと、
   前記複数のキュービットの少なくとも一部分を測定することで、前記量子システムの読み出しを取得するステップと、
   前記読み出しから前記計算問題の解を決定するステップと、
   を含むことを特徴とする量子処理方法。
2. 前記量子システムを前記初期量子状態から前記最終ハミルトニアンの前記基底状態に向けて移行させるステップが、量子アニーリングによって前記量子システムの初期ハミルトニアンから前記最終ハミルトニアンに移行するステップを含むことを特徴とする請求項１に記載の量子処理方法。
3. 前記量子システムを前記初期ハミルトニアンの基底状態に向けて冷却することにより、前記初期量子状態の前記量子システムを初期化するステップを更に含むことを特徴とする請求項２に記載の量子処理方法。
4. 前記複数の調整可能なパラメータが、前記複数のキュービットに作用する単体フィールドの複数のフィールド強度及び／又は複数のフィールド方向を含むことを特徴とする請求項１から３の何れかに記載の量子処理方法。
5. 前記初期ハミルトニアンから前記最終ハミルトニアンへの量子アニーリングを実行するステップが、前記初期ハミルトニアンを前記最終ハミルトニアンに断熱的に移行させるステップを含むことを特徴とする請求項１から４の何れかに記載の量子処理方法。
6. 前記短距離ハミルトニアンは、ｄ＝４であるｄ体ハミルトニアンであることを特徴とする請求項１から５の何れかに記載の量子処理方法。
7. 前記短距離ハミルトニアンは、前記計算問題から独立していることを特徴とする請求項１から６の何れかに記載の量子処理方法。
8. 前記複数のキュービットが２次元格子又は３次元格子に沿って配置されることを特徴とする請求項１から７の何れかに記載の量子処理方法。
9. 前記複数のキュービットが２次元格子に沿って配置され、前記短距離ハミルトニアンが、前記２次元格子のプラケットに対応する４つのキュービットからなるグループ間での相互作用を伴うことを特徴とする請求項１から８の何れかに記載の量子処理方法。
10. 前記計算問題がＮＰ困難問題であることを特徴とする請求項１から９の何れかに記載の量子処理方法。
11. 前記問題ハミルトニアンが式Σ_ｋＪ_ｋσ_ｚ ^（ｋ）を有し、σ_ｚ ^（ｋ）が前記複数のキュービットのｋ番目のキュービットのパウリ演算子であり、各Ｊ_ｋが係数であり、前記係数Ｊ_ｋが前記問題ハミルトニアンの前記複数の調整可能なパラメータを形成することを特徴とする請求項１から１０の何れかに記載の量子処理方法。
12. 前記問題コード化設定を決定するステップが、
    前記計算問題を予備の計算問題にマッピングするステップであって、前記予備の計算問題が長距離スピンモデルの基底状態（ｍが１、２又は３のｍ体相互作用を有する長距離スピンモデルの基底状態）を決定するステップを含む、前記マッピングするステップと、
    前記長距離スピンモデルから前記問題コード化設定を決定するステップと、
    を含むことを特徴とする請求項１から１１の何れかに記載の量子処理方法。
13. 前記長距離スピンモデルにおけるスピンの複数の閉ループから前記短距離ハミルトニアンを決定するステップを更に含むことを特徴とする請求項１２に記載の量子処理方法。
14. 計算問題の解を計算する装置であって、
    複数のキュービットを有する量子システムと、
    前記量子システムの基底状態に向けて前記量子システムを冷却するのに適した冷却ユニットと、
    量子アニーリングによって前記量子システムの初期ハミルトニアンを前記量子システムの最終ハミルトニアンに移行させるのに適したプログラム可能な量子アニーリングユニットであって、前記最終ハミルトニアンが問題ハミルトニアン及び短距離ハミルトニアンの総和であり、前記問題ハミルトニアンが複数の調整可能なパラメータを有する単体ハミルトニアンである、量子アニーリングユニットと、
    前記複数のキュービットの少なくとも一部分を測定するのに適した測定装置と、
    前記プログラム可能な量子アニーリングユニット及び前記測定装置に接続された古典的計算システムと、を備え、
    前記古典的計算システムは、
    計算問題を入力として受信するステップと、
    前記計算問題を前記問題ハミルトニアンにコード化するステップであって、前記計算問題から前記問題ハミルトニアンの前記複数の調整可能なパラメータに対する問題コード化設定を決定するステップを含む、前記コード化するステップと、
    前記問題コード化設定を前記量子アニーリングユニットに伝達するステップと、を実行するように構成され、
    前記プログラム可能な量子アニーリングユニットは、
    前記古典的計算システムから前記問題コード化設定を受信するステップと、
    量子アニーリングによって前記初期ハミルトニアンを前記最終ハミルトニアンに移行させるステップであって、前記問題ハミルトニアンの前記複数の調整可能なパラメータが前記問題コード化設定に存在する、前記移行させるステップと、を実行するように構成され、
    前記古典的計算システムは、
    前記測定装置から前記量子システムの読み出しを受信するステップと、
    前記読み出しから前記計算問題の解を決定するステップと、を更に実行するように構成されていることを特徴とする量子処理装置。
15. 計算問題の解を計算するためのプログラム可能な量子アニーリング装置であって、
    ２次元格子に沿って配置された複数の超電導キュービットを有する量子システムと、
    複数の調整可能な磁束を生成するように構成された複数の磁束バイアスユニットを有する磁束バイアスアセンブリであって、各調整可能な磁束が前記複数の超電導キュービットのうち単一の超電導キュービットに作用する、磁束バイアスアセンブリと、
    プラケットハミルトニアンに従い、前記複数の超電導キュービットを結合するように構成された少なくとも１つの超電導量子干渉装置を有する結合ユニットと、
    前記磁束バイアスユニット及び前記結合ユニットに接続された制御装置と、備え
    前記制御装置は、
    前記量子システムの問題ハミルトニアンの複数の調整可能なパラメータに対する問題コード化設定を受信するステップであって、前記問題ハミルトニアンが単体ハミルトニアンであり、前記問題コード化設定が計算問題をコード化する、前記受信するステップと、
    前記磁束バイアスアセンブリ及び前記結合ユニットを制御して、量子アニーリングによって前記量子システムの初期ハミルトニアンを前記量子システムの最終ハミルトニアンに移行させるステップであって、前記最終ハミルトニアンが前記プラケットハミルトニアン及び前記問題ハミルトニアンの総和であり、前記問題ハミルトニアンの前記複数の調整可能なパラメータが前記問題コード化設定に存在する、前記制御するステップと、を実行するように構成されていることを特徴とする量子アニーリング装置。

Images