JP2022510394A - 量子回路において２キュービットゲートを低減するための方法および装置 - Google Patents

Abstract

本開示は、量子回路における２キュービットゲートを低減するための方法、装置、コンピュータ可読媒体、および／または手段を説明し、この方法は、量子回路を形成する第一の複数の２キュービット量子ゲートに関連する情報を含むネットリストを受信するステップと、第一の複数の２キュービット量子ゲートと機能的に同等である第二の複数の２キュービット量子ゲートを生成するために、第一の複数の２キュービット量子ゲートに関連する情報に対して制御されたゲート相殺オペレーションを実行するステップと、第二の複数の２キュービット量子ゲートに関する情報を含む新規ネットリストを生成するステップと、第二の複数の２キュービット量子ゲートに基づいて、量子回路の機能をインプリメントするために新規ネットリストを提供するステップと、を含み得る。【選択図】図２４

Description

（関連出願の相互参照）
本特許出願は、２０１９年１１月８日に出願された「量子回路において２キュービットゲートを低減するための方法および装置」と題する米国非仮出願第１６／６７８，８３５号と、２０１８年１２月７日に出願された「量子回路において２キュービットゲートを低減するための方法および装置」と題する米国仮特許出願第６２／７７６，６３４号との優先権を主張し、両出願の内容は、参照によりその全体が本明細書に組み込まれる。

本開示の態様は、一般に、回路の最適化に関連し、より具体的には、連続パラメータを用いた大規模量子回路の自動化された最適化に関連する。

大規模量子回路は、異なる方法で実施することができる。トラップされた原子イオンの使用は、普遍的で完全にプログラム可能な量子コンピューティングマシンを提供してきた量子情報処理（ＱＩＰ）アプローチの１つである。トラップされた原子イオンは、量子情報ネットワーク（ＱＩＮ）の主要なプラットフォームでもある。そのようなシステムまたはネットワークの全体的な通信を改善することができるトラップされた原子イオンに基づくシステムまたはネットワークが望ましい。

したがって、トラップされた原子イオンに基づくＱＩＰシステムおよび超電導回路に基づくＱＩＰシステムを含む様々なタイプのＱＩＰシステムにおける大規模量子回路のより良い最適化（例えば、より少ないゲート数）を可能にする効率的な技術を開発することが望ましい。

以下では、そのような態様の基本的な理解を提供するために、１つ以上の態様の簡略化された概要を提示する。この概要は、考えられる全ての態様の広範な概要ではなく、全ての態様の主要または重要な要素を特定することを意図するものでも、一部または全ての態様の範囲を説明することを意図するものでもない。その目的は、後で提示されるより詳細な説明の前置きとして、１つ以上の態様のいくつかの概念を簡略化された形式で提示することである。

本明細書に記載されるのは、連続パラメータを用いた大規模量子回路の自動最適化のための手法である。例えば、本開示の態様は、古典的コンピュータよりも優れた量子計算で期待されるサイズおよびタイプの量子回路を最適化するための自動化された方法のインプリメンテーションを説明する。ここで説明する手法は、連続ゲートパラメータを処理する方法を示しており、大規模量子回路を最適化できる高速アルゴリズムのコレクションに適用できる。これらの手法は、基礎となる量子アルゴリズムの基本的なレイアウトを維持するために最小限の構造変更を行いながら、以前のアプローチよりも大幅に短い時間でより良い最適化を提供できる。これらの手法によって提供される結果は、既存の量子コンピューティングハードウェアで実行できる計算と、量子コンピューティングハードウェアでのインプリメンテーションがより困難であるが、古典的コンピュータで達成できるものよりも優れていると期待される、より高度な計算との間のギャップを埋めるのに役立つ。

一例では、量子回路を最適化するための方法が説明され、この方法は、量子回路を形成する第一の複数の２キュービット量子ゲートについての情報を含むネットリストを受信するステップと、量子ゲートの第一のリストとの機能等価性を有する量子ゲートの第二のリストを生成するために、量子ゲートの第一のリストについての情報に位相多項式低減オペレーションを実行するステップであって、量子ゲートの第二のリスト中の量子ゲートの数などの計算コストは、量子ゲートの第一のリスト中の量子ゲートの数などの計算コストよりも小さい、ステップと、量子ゲートの第二のリストについての情報を含む新規ネットリストを生成するステップと、量子ゲートの第二のリストを使用して量子回路の機能をインプリメントするための新規ネットリストを提供するステップと、を含む。

いくつかのインプリメンテーションでは、量子回路内の２キュービットゲートを低減するための方法は、量子回路を形成する第一の複数の２キュービット量子ゲートに関連する情報を含むネットリストを受信するステップと、第一の複数の２キュービット量子ゲートと機能的に同等である第二の複数の２キュービット量子ゲートを生成するために、第一の複数の２キュービット量子ゲートに関連する情報に対して制御されたゲート相殺オペレーションを実行するステップであって、第一の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第一の数は、第二の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第二の数よりも大きい、ステップと、第二の複数の２キュービット量子ゲートに関する情報を含む新規ネットリストを生成するステップと、第二の複数の２キュービット量子ゲートに基づいて、量子回路の機能をインプリメントするために新規ネットリストを提供するステップと、を含み得る。

添付の図面は、一部のインプリメンテーションのみを示しており、したがって、範囲を制限するものと見なされるべきではない。

本開示の態様による量子回路の一例を示す図表である。 本開示の態様による位相多項式表現のための量子回路の一例を示す図表である。 本開示の態様によるアダマール（Ｈａｄａｍａｒｄ）ゲート低減のためのルールの例を示す図表である。 本開示の態様による位置転換ルールの例を示す図表である。 本開示の態様による、位相多項式を用いる回転併合による最適化サブルーチンの一例を示す図表である。 本開示の態様による位相多項式表現に適うサブ回路を示す図表である。 本開示の態様によるゲートカウント保持の書換えルールの例を示す図表である。 本開示の態様による、ゲートカウント低減の書換えルールの例を示す図表である。 本開示の態様による図５中の回路のさらなる単純化の一例を示す図表である。 本開示の態様によるトフォリ（Ｔｏｆｆｏｌｉ）ゲートのインプリメンテーションの例を示す図表である。 本開示の態様による、近似量子フーリエ変換（ＱＦＴ：ｑｕａｎｔｕｍ Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ）、クイッパー（Ｑｕｉｐｐｅｒ）ライブラリ加算回路（ａｄｄｅｒ）、およびフーリエベース加算回路（ＱＦＡ：Ｆｏｕｒｉｅｒ－ｂａｓｅｄ ａｄｄｅｒ）に対する合計ゲートカウントの例を示すチャートである。 本開示の態様による、クイッパーライブラリ加算回路に対するＣＮＯＴゲートのいくつかの例を示すチャートである。 本開示の態様による諸加算回路の軽度最適化に対する結果を示す表である。 本開示の態様による、クイッパーライブラリ加算回路に対する重度最適化の結果を示す表である。 本開示の態様による積公式アルゴリズムに対する合計ゲートカウントの例を示すチャートである。 本開示の態様による積公式アルゴリズムの最適化に対する結果を示す表である。 本開示の態様による積公式アルゴリズムの最適化に対する結果を示す表である。 本開示の態様による、各種のアルゴリズムに対する軽度および重度のＴ－ｐａｒ最適化の比較を示す表である。 本開示の態様による、各種のアルゴリズムに対する軽度および重度のＴ－ｐａｒ最適化の比較を示す表である。 本開示の態様による、各種のアルゴリズムに対する軽度および重度のＴ－ｐａｒ最適化の比較を示す表である。 本開示の態様による、連続するパラメータを備える大規模量子回路の自動最適化を実行するためのデバイスの一例を示すブロック図である。 本開示の態様による、連続するパラメータを備える大規模量子回路の自動最適化の方法の一例を示すフローチャートである。 本開示の態様による、最適化された大規模量子回路をインプリメントするためのトラップイオンベースの量子情報処理（ＱＩＰ）システムの一例を示すブロック図である。 本開示の態様による２キュービット量子ゲート低減のためのルールの例を示す図である。 本開示の態様による制御されたゲート相殺のためのルールの例を示す図である。 本開示の態様による制御されたゲート相殺のためのルールの例を示す図である。 本開示の態様による等価量子回路の図を示す。 本開示の態様による、量子回路内の２キュービット量子ゲート低減のための方法の一例を示すフローチャートである。

添付の図面に関連して以下に記載される詳細な説明は、様々な構成の説明として意図されており、本明細書に記載の概念が実施され得る唯一の構成を表すことを意図するものではない。詳細な説明には、様々な概念を完全に理解するための具体的な詳細が含まれている。しかしながら、これらの概念は、これらの特定の詳細がなくても実施され得ることは当業者には明らかであろう。場合によっては、そのような概念を曖昧にすることを避けるために、周知のコンポーネントがブロック図の形式で示される。

本開示は、古典的コンピュータよりも優れた量子計算で期待されるサイズおよびタイプの量子回路を最適化するための自動化された方法のインプリメンテーションを説明する。本開示は、連続ゲートパラメータを処理する方法を示しており、大規模量子回路を最適化できる高速アルゴリズムのコレクションを報告する。検討した一連のベンチマークについて、説明される手法は、ゲート数を大幅に削減できる。特に、本開示の手法は、基礎となる量子アルゴリズムの基本的なレイアウトを維持するために最小限の構造変更を行いながら、以前のアプローチよりも大幅に短い時間でより良い最適化を提供する。これらの手法によって提供される結果は、既存の量子コンピューティングハードウェアで実行できる計算と、量子コンピューティングハードウェアでのインプリメンテーションがより困難であるが、古典的コンピュータで達成できるものよりも優れていると期待される、より高度な計算との間のギャップを埋めるのに役立つ。

（序章）
量子コンピュータまたは量子情報処理（ＱＩＰ）システムは、特定の問題を解決するにあたって従来式のコンピュータをはるかに凌駕する潜在能力を有する。おそらく、これらのもっと知られた応用は整数の因数分解作業であって、既知の最速の従来式アルゴリズムは超多項式であり、ショア（Ｓｈｏｒ）のアルゴリズムは多項式時間でこの問題を解き、広く使われているＲＳＡ暗号システムへの攻撃法を提供する。

ショアのアルゴリズムが発見される前でさえ、量子力学をシミュレートするために量子コンピュータが提案されていた。ハミルトニアン（Ｈａｍｉｌｔｏｎｉａｎ）力学をシミュレートすることによって、量子コンピュータは、物性物理学および高エネルギ物理学、量子化学、および材料科学を探求ことができる。量子シミュレーションの有用なインスタンスは、因数分解問題の従来の難題インスタンスよりも、より小規模な量子コンピュータへのアクセス可能性が高い。

これらのおよび他の潜在的応用は、スケーラブル量子コンピュータの構築に向けての大きな努力をモチベートする助力となってきた。２つの量子コンピューティング技術、超電導回路とトラップイオンとは、現在は程々のサイズではあるが、完全にプログラム可能なユニバーサルなデバイスを可能にするのに十分なほどに成熟してきた。いくつかのグループは、産業界および政府の相当な投資に後援されて、これらのプラットフォームをより大規模なデバイスに積極的に発展させている。こうして、数十、あるいはさらには数百キュービットが関わる量子計算がそう遠くない将来に実行される可能性が高い。

実験的な量子情報処理は、まだ困難な技術課題に留まっており、しばらくの期間は、量子計算のために利用可能なリソースは高価で非常に限定されることになろう。利用可能なハードウェアを最大限に生かすために、できるだけ効率の良い（例えば、最小数のゲートでインプリメントされる）量子アルゴリズムのインプリメンテーションを開発することが必須である。

量子アルゴリズムは、通常、量子回路に関係づけて表され、キュービットに作用する基本量子ロジックゲートのシーケンスとして計算を表現する。所与のアルゴリズムを、基本的オペレーションの利用可能なセットとともにインプリメントする多くの仕方があり、最少のリソースを使用するインプリメンテーションを見出すために有用である。抽象的な意味で、効率的なアルゴリズムを開発し、実践的な効率に目を向けて、そのアルゴリズムを実行することが重要であるが、大規模量子回路は、自動最適化を可能にする技法から利益を得るのに十分な複雑性を有することが多い。

本開示では、マニュアルでのゲートレベルの最適化がもはや実際的には行えないスケールで量子回路のパフォーマンスをできるだけ改良することを目指して、それら量子回路のサイズを低減するためのソフトウェアツール（例えば、量子回路オプティマイザ）としてインプリメントできる様々な技法が説明される。任意の量子回路の全体的最適化はＱＭＡ困難なので、本開示で説明するアプローチは、ゲートカウントを低減するための自動化方式の注意深く選択された経験則のセットを適用することであり、多くの場合相当な節減をもたらす。

本明細書で記述する最適化技法は、いくつかの型の量子回路に適用が可能である。ベンチマーク回路は、量子フーリエ変換、整数加算回路、ガロア体乗数など、因数分解および離散対数を計算するための量子アルゴリズムのコンポーネントを含む。さらに検討されたのは、ハミルトニアンシミュレーションへ積公式アプローチのための回路である。これらのケースの各々において、焦点は、従来式の計算を凌駕するアプリケーションにおいて有用な可能性の高い回路のサイズに当てられた。本明細書に記載の諸技法は、実践者が、所与のアプリケーションにおいて、どのアプリケーションのインプリメンテーションが最も効率的かを理解するのに役立つことができる。

量子回路の最適化に関する既存の研究はあるが、本明細書で検討しているような大規模回路を目標とした自動最適化技法に焦点を当てた研究は限られている。さらに、後記でさらに詳しく説明するように、これまでに報告されたランタイムからの推定は、既存の量子回路オプティマイザがかかる大規模回路に対しうまく機能する可能性が低いことを示唆している。本開示で提案する技法を用いて最適化された回路と、他のアプローチとの直接的な比較は、本提案の技法が、一般に、より短時間でより小さな回路を見出すことを示している。さらに、本開示で提案する技法は、連続するパラメータを備えた量子回路の自動最適化のために使用される。

（全体的考察）
本開示では、大規模量子回路、すなわち、従来式のコンピュータの力の及ばない量子計算において出現する量子回路を効率的に最適化する課題が検討される。本開示は、ソフトウェアソリューション（例えば、大規模量子回路の最適化に対するオフラインソリューション）としてインプリメンテーションおよび実行が可能な軽度および重度最適化アルゴリズム（例えば、オプティマイザの軽度バージョンおよび重度バージョン）に関する２つの最適化アプローチを説明する。これらのアルゴリズムは、基本的最適化の選択されたシーケンスに基づくが、それでもこれらは、ゲートカウントのかなりの低減を成就し、Ｔ－ｐａｒ最適化（後記でさらに詳しく説明する）など、より数学的に洗練されたアプローチを上回る改良を行う。当アプローチの簡潔性は、特にそのオプティマイザの軽度バージョンを使うことによって、非常に高速なランタイムに反映される。

本オプティマイザの重度バージョンは、より大きな回路の最適化が可能であることを実証している。出力をさらに改良するために、可能ならば一部のアダマールゲートを含めて、ＣＮＯＴおよびＲ_Ｚゲートのステージを構成するためのより拡張された（しかして計算的によりきつい）アルゴリズムをインプリメントすることによって、Ｒ_Ｚのカウントを低減するためルーチンを修改することが可能である。別の考えは、本明細書に記載のアルゴリズム中にテンプレートベースのピープホール最適化を組み込むことである。他の考えは、サブ回路再書込みルール（後記でさらに詳しく説明する）のセットを拡張し、他のベンチマーク回路へのアプローチのそのパフォーマンスを検討することである。最後に、各種リソース（例えば、各種のゲート、補助キュービット）の相対コストを検討することにより、これらリソースの有利なトレードオフへの最適化につなげることができよう。

（方法）
本開示で提案する様々な最適化アルゴリズムおよびそれらのインプリメンテーションに関する詳細を以降で説明する。本開示全体を通して、回路に出現するゲートの数を示すために用語ｇが用いられる。以下の背景節では、本開示全体を通して使われる表記の定義が提示される。量子回路の表現と題された節では、本明細書に記載の技法と関連させて用いることの可能な量子回路の３つの別個の表現を説明する。前処理と題された節では、本明細書に記載の様々な最適化アルゴリズムと併せ使用が可能な前処理ステップを説明する。最適化サブルーチンと題された節では、本開示で提案されるアプローチの基本構築ブロックを形成するいくつかのサブルーチンを説明する。さらに、汎用最適化アルゴリズムと題された節では、最適化アルゴリズムの様々なバージョンを形成するために、これらのサブルーチンをどのように組み合せるかを説明する。最後に、特殊用途最適化と題された節では、特定種類の回路を取り扱うために使用が可能な、２つの特殊用途最適化技法を提示する。

（背景）
量子回路は、キュービットの集合に作用する一連の量子ゲートである。量子回路は、便宜上、水平のワイヤがキュービットの時間発展を示し、時間は左から右に伝搬し、ボックス（またはワイヤに繋がれた他の記号）が量子ゲートを表す図によって表現される。例えば、図１の図表１００は簡単な３キュービットの量子回路を表している。図１００中の回路は２つの単一キュービットのｚ回転ゲート（１１０ａ、１１０ｂ）、Ｒ_Ｚ（θ）およびＲ_Ｚ（θ’）、２つの単一キュービットのアダマールゲート（１２０ａ、１２０ｂ）Ｈ、および４つの２キュービットの制御ＮＯＴゲート（１３０ａ、１３０ｂ、１３０ｃ、および１３０ｄ）ＣＮＯＴを含む。

量子回路に対する基本ゲートの簡単なセットは、図表１００中の回路に示されるように、２キュービットＣＮＯＴ、ならびに単一キュービットＮＯＴゲート、単一キュービットのアダマールゲート、および単一キュービットｚ回転ゲートから成ると考えてよい。これらの種類のゲートに対するユニタリ行列は、次の形の表現式を取る。

上式のθ∈０，２π］は回転角度である。位相ゲートＰおよびＴゲートは、それぞれ、Ｒ_Ｚ（π／２）とＲ_Ｚ（π／４）として、Ｒ_Ｚ（θ）から上限検出不可能な大域的位相から得ることができる。回転角度が関係しない場合、ｚ回転は、総称的にＲ_Ｚとして表すことができる。

本明細書で説明する技法は、Ｈ、Ｒ_Ｚ、ＮＯＴ、およびＣＮＯＴゲートのセットによる量子回路を主として提示するが、トフォリゲートを含むことが可能な入力回路も検討することができる。トフォリゲート（例えば、図１０の図表１０００中のゲート１０１０ａ）は、計算ベースの状態で、

をマッピングすることによって表される。また、トフォリゲートは、否定制御を有することも可能である。例えば、最上部制御が否定であるトフォリゲート（例えば、図１０の図表１０００中のゲート１０１０ｂ）は、

として作用し、否定制御の両トフォリゲート（例えば、図１０の図表１０００中のゲート１０１０ｃ）は、

として作用する。

所与の量子回路を実行するコストは、それをインプリメントするために使われる物理システムにより決まる。例えば、物理システムが、超電導回路に基づいているか、またはトラップ原子イオンに基づいているかによって異なったコスト考慮がされることになろう。また、このコストは、物理レベル（未保護状態）インプリメンテーションとロジカルレベル（フォールトトレラント状態）のインプリメンテーションとの間でも大きく変わり得る。物理レベルでは、２キュービットのゲートは、通常、単一キュービットのゲートよりもインプリメンテーションの費用が高い。本明細書で説明する技法は、ＣＮＯＴゲートカウントを検討し、最適化アルゴリズム中のＣＮＯＴゲートの数を最適化することによってこれに対応する。

ロジカルレベルのフォールトトレラントな回路に対し、いわゆるクリフォード（Ｃｌｉｆｆｏｒｄ）演算（アダマール、位相、およびＣＮＯＴゲートによって生成される）は、多くの場合、比較的インプリメンテーションが容易であり、非クリフォード演算は相当なオーバーヘッドを要する可能性がある。こうして、本最適化アルゴリズムにおいては、Ｒ_Ｚゲートの数も検討され、それらのカウントも最適化するよう試みられる。フォールトトレラントなインプリメンテーションにおいて、Ｒ_Ｚゲートは、通常、クリフォードおよびＴゲートから成るディスクリートゲートのセットで近似される。かかる近似を生成するため最適のアルゴリズムが知られている。一般的なＲ_Ｚゲートを近似するために必要なクリフォード＋Ｔゲートの数は、回転の具体的角度よりも、むしろ主として望まれる精度の如何により、したがって、そのＲ_Ｚゲートをクリフォード＋Ｔのフォールトトレラントな回路で近似する前に回路を最適化するのが望ましい。

ＣＮＯＴおよびＲ_Ｚのカウントを最小化することによって、物理およびロジカルレベル両方のインプリメンテーションを目指した最適化を行うことが可能である。これら２つの目標の間でのトレードオフを予期してよく、かかるトレードオフが行われた諸インスタンスがある。但し、本明細書に記載の技法では、Ｒ_ＺおよびＣＮＯＴ両方のカウントだけを目指した最適化が検討される。

（量子回路の表現）
量子回路の以下の３つの表現式は、本明細書で説明する最適化技法または最適化アルゴリズムに関連付けて使用することが可能である。

第一に、回路は、逐次に適用されるゲートのリスト（ネットリスト）として格納することができる。場合によっては、サブルーチン（例えば、回路サブルーチン）に関連させて、回路を特定するのが便利であり、これはブロックとして表すことが可能である。各ブロックは、任意の回数反復することが可能で、回路中に存在するキュービットの任意のサブセットにも適用できる。ブロックを用いる表現は、多くの量子回路が相当な数の繰り返しを示すので、特に簡潔にすることができる。ブロックは、ゲートおよびキュービットのアドレスのリストとして特定される。

ネットリスト表現は、いくつかのベンチマーク回路を特定するために用いられた量子プログラミング言語である、クイッパー（Ｑｕｉｐｐｅｒ）によって生成された」フォーマットを使って入力、出力することができる。このフォーマットは、ブロックを取り扱う機能を含む。同じ機能を有する他の量子プログラミング言語および／またはフォーマットを使うことも可能である。

第二に、有向非巡回グラフ（ＤＡＧ：ｄｉｒｅｃｔｅｄ ａｃｙｃｌｉｃ ｇｒａｐｈ）表現を用いることができる。ＤＡＧの角頂は、回路のゲートであり、エッジはそれらの入力／出力関係を符号化する。ＤＡＧ表現は、ゲートの間の隣接性へのアクセスを容易にする利点がある。

第三に、｛ＮＯＴ、ＣＮＯＴ、Ｔ｝回路の位相多項式表現の一般化も用いることができる。ネットリストおよびＤＡＧ表現と違って、この最後の表現は、ＮＯＴ、ＣＮＯＴ、およびＲ_Ｚゲートだけで構成される回路に適用される。かかる回路は、線形ブール変換および対角相変換の編成として簡潔に表現することができる。例えば、Ｃが、ＣＮＯＴゲートおよび、ゲートＲ_Ｚ（θ_１）、Ｒ_Ｚ（θ_２）、…、Ｒ_Ｚ（θ_ｉ）だけから成る回路であるとする。このとき、ｎキュービットベースの状態｜ｘ_１，ｘ_２，…，ｘ_ｎ〉でのＣのアクションは、次の形の式：

を有する。上式のｈ：｛０，１｝^ｎ→｛０，１｝^ｎは線形可逆関数であり、

は、アフィンブール関数ｆ_ｉ：｛０，１｝^ｎ→｛０，１｝と係数低減モジュロ２πとの線形組み合せである。項ｐ（ｘ_１，ｘ_２，…，ｘ_ｎ）は、回路Ｃに関連する位相多項式と言われる。例えば、図２の図表２００中に示された、４つの単一キュービットｚ回転ゲート（２１０ａ、２１０ｂ、２１０ｃ、および２１０ｄ）と、３つの２キュービットＣＮＯＴゲート（２３０ａ、２３０ｂ、および２３０ｃ）とを有する回路は、以下のマッピング表現：

で表すことができ、ここで、

である。いくつかのインスタンスでは、位相多項式表現は、｛ＣＮＯＴ、Ｔ｝回路に対してだけ考慮され、こうして、式（３）中の全てのθ_ｉは、π／４の整数倍数であり、関数ｆ_ｉは線形である。

上記の３つの回路表現の任意の２つの間で、時間的に線形に回路中のゲートの数で変換することが可能である。ネットリストが与えられると、対応するゲート毎のＤＡＧを作成することが可能である。逆に、標準的トポロジカル分類によってＤＡＧをネットリストに変換することも可能である。ネットリストと、｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝回路の位相多項式表現との間の変換は、例えば、例として、深さ最適の量子回路の高速合成に使われる中間突合せ型アルゴリズムに関連して記述される直接的な一般化を用いることによって可能となる。

（前処理）
主たる最適化手順（例えば、最適化アルゴリズム）を実行する前に、当該回路を、さらなる最適化をやり易くするために前処理することができる。本オプティマイザは、ＮＯＴ、ＣＮＯＴ、およびトフォリゲートに加え、ＨおよびＲ_Ｚ（θ）ゲートを扱うことができるので、この前処理は、ＮＯＴ、ＣＮＯＴ、トフォリ、Ｈ、および／またはＲ_Ｚ（θ）ゲートから成る入力回路に適用することが可能である。いくつかの例として、後記でより詳しく説明するクイッパー加算回路およびＴ－ｐａｒ回路のベンチマークがある。例えば、ＮＯＴゲートは、それらを、トフォリゲート制御と、トフォリゲートおよびＣＮＯＴゲートの標的とを通り抜けて交換することによって、できるだけ右方に押しやることができる。トフォリゲート制御を通り抜けてＮＯＴゲートを押しやる場合、その制御は、否定となる（または、それが当初に否定であった場合は、その否定が除去される）。この手順が隣接するＮＯＴゲートの対に及ぶ場合、ＮＯＴゲートを回路から除去することができる。かかる相殺が見出せない場合、その制御否定の変更は、逆戻りさせることが可能で、当該ＮＯＴゲートは、その当初の位置に戻すことができる。さらに、ＮＯＴゲートとＨゲートとの間のこの位置転換の関係は、Ｈが後続するＮＯＴは、ＨにＺゲートが後続することであり、Ｚゲートは、Ｒ_Ｚ（θ＝π）である。また、Ｒ_Ｚ（θ）が後続するＮＯＴは、ＮＯＴが後続するＲ_Ｚ（－θ）であるようなものである。したがって、ＮＯＴ、ＣＮＯＴ、トフォリ、Ｈ、およびＲ_Ｚゲートから成るゲートの標準的セットに対し、ＣＮＯＴおよびＲ_Ｚゲートを通り抜けてのＮＯＴの伝搬を用いることによって位相の全てのアフィン関数を線形関数に低減することに関連して後記で説明する詳細を加えて、前処理ステップの一部として、ＮＯＴゲートは、所与の回路の開始部または終端部に移動することが可能である。

このＮＯＴゲートの伝搬は、本提案のオプティマイザまたは最適化アルゴリズムの２つの側面を利用する。第一に、否定制御を有している可能性のあるトフォリゲートを受入れ、後記でさらに詳しく説明するように、Ｔ／Ｔ^†極性の選択の自由度を利用することによって、クリフォード＋Ｔ回路中への該ゲートの分解を最適化することが可能である。第二に、ＮＯＴゲートの相殺は位相多項式表現を（位相多項式表現（式３）中の関数ｆ_ｉの一部を、元のアフィンの代わりに線形にすることによって）簡素化するので、かかる相殺は、後記で説明するサブルーチン４および５が（これらのサブルーチンが位相多項式表現中に一致する項を見出すことに依拠しているので）最適化を見出すことの可能性をより高くする。

この前処理ステップの複雑性は、このステップが回路を通る単一回の通過を行うだけなので、Ｏ（ｇ）である。

（最適化サブルーチン）
本開示の最適化アルゴリズムは、以下でさらに詳細に説明する様々なサブルーチンに依拠している。各々のアイテムに対し、回路中のゲートの数ｇの関数として最悪の場合の時間複雑性が報告されている（簡潔化のため、キュービットの数および他のパラメータへの依拠は無視されている）。但し、実践上でのソフトウェアのランタイムは、後記でさらに説明するように、これらサブルーチンを注意深く順序付け制限することによって、最適化することが可能である。

（サブルーチン１：アダマールゲートの低減）
アダマールゲートは、位相多項式最適化（後記サブルーチン４および５）に関与せず、また、ゲート交換を妨げる傾向がある。しかして、図３の図表３００中に描かれた回路アイデンティティまたはルールが、アダマールゲートのカウントを低減するために用いられる。このアイデンティティまたはルールは、図表３００の下部にルール３１０ａおよび３１０ｂと、図表３００の上部にルール３２０ａ、３２０ｂ、および３２０ｃを含む。ルール３１０ａおよび３１０ｂは、たとえ中間のＣＮＯＴが無作為に多くのＣＮＯＴゲートを備える回路に置き換えられても、それらの全てが当初のＣＮＯＴゲートの標的を共有しているならば、適用が可能である。

これらのルールの各々の適用は、上限４つまでＨのカウントを低減する。所与のアダマールゲートに対し、ＤＡＧ表現を使って、それがこれらの回路アイデンティティの１つに当てはまるかどうかを一定の時間ごとにチェックすることが可能である。しかして、回路中の全てのアダマールゲートを通って単一回通過することによって、このサブルーチンを複雑度Ｏ（ｇ）でインプリメントすることが可能である。

（サブルーチン２：単一キュービットゲートの相殺）
一般に、量子回路のＤＡＧ表現を用いて、ゲートとその逆ゲートが隣接しているかどうかが直接に判断される。そうである場合、ゲートカウントを低減するために両方のゲートを除去することが可能である。さらに一般的には、ゲートＵと可換であるサブ回路Ａによって分離されている２つの単一キュービットゲートＵとＵ^†とを相殺することが可能である。一般に、ゲートＵが回路Ａと可換であるかどうかの判断は、計算量が膨大となる可能性がある。その代わりに、十分な（だが必須ではない）転換に対する条件を提供する特定のルールのセットを適用することが可能である。このアプローチは迅速であり、量子回路を単純化するため利用可能な多くの位置転換を見出すように見える。

具体的に、回路中の各ゲートＵに対し、オプティマイザまたは最適化アルゴリズムは、Ｕ^†の何らかのインスタンスとの可能な相殺を探索する。これを行うために、このアプローチは、図４の図表４００中パターンの１つによって実証されるように、Ｕが、引き続くゲートのセットを通り抜けて可換かどうかを繰り返しチェックすることになる。図表４００は、Ｒ_Ｚゲートを右方に位置転換するための、上側のルール４１０ａ、４１０ｂ、および４１０ｃと、ＣＮＯＴゲートを右方に位置転換するためのルール４２０ａ、４２０ｂ、および４２０ｃなどの位置転換ルールを含む。何らかの段階で、許容された通信パターンによって、Ｕを右方に移動することができない場合、Ｕと合致するＵ^†との相殺はできず、当初の構成が復元される。それ以外の場合は、ＵとＵ^†の何らかのインスタンスとの相殺は成功する。前述した位置転換スキームは、右方への特定方向だけに適用される必要はなく、左方の特定方向にも適用されてよい。

ｇ個のゲートＵの各々に対し、本サブルーチンは、それがＯ（ｇ）を通して後続の位置と可換であるかどうかをチェックする。しかして、全体的なゲート相殺ルールの複雑性はＯ（ｇ^２）となる。この複雑性は、固定数の後続のゲートを通しての位置転換だけを検討することによって、ｇに対し線形にできようが、実践上はその必要がないことが判明している。

本サブルーチンの些少なバリエーションを使って、相反ゲートを相殺するのでなく、回転ゲートを併合することが可能である。具体的に、２つの回転Ｒ_Ｚ（θ_１）およびＲ_Ｚ（θ_２））を単一の回転Ｒ_Ｚ（θ_１＋θ_２）に組み合せ、１つのＲＺゲートを削除することができる。

（サブルーチン３：２キュービットゲートの相殺）
このサブルーチンは、Ｕが２キュービットゲートであることを除き、サブルーチン２と類似であり、検討対象の回路では、これらのゲートは一般ＣＮＯＴである。同様に、本サブルーチンの複雑性はＯ（ｇ^２）であるが、サブ回路Ａに対し最大限のサイズを取らせることによって、Ｏ（ｇ）に低減することが可能である。

（サブルーチン４：位相多項式を用いる回転併合）
このケースでは、ＮＯＴ、ＣＮＯＴ、およびＲ_Ｚゲートから成るサブ回路を検討する。その位相多項式の２つの別個の項が、何らかのｉ≠ｊに対し、次式：

を満たす場合、対応する回転Ｒ_Ｚ（θ_ｉ）とＲ_Ｚ（θ_ｊ）とを併合することが可能である。例えば、図２の図表２００中に示された回路において、第一および第四の回転の両方は、回路の位相多項式表現によって明らかなように、値ｙを携持するキュービットに適用される。しかして、図２の図表２００中の回路は、図５の図表５００中に示される変換を経て、その中でこれら２つの回転は組み合される。図表５００中の右側図に変換された回路は、Ｒ_Ｚゲート５１０ａ、５１０ｂ、および５１０ｃを含み、Ｒ_Ｚゲート５１０ｃは、組み合せされた２つの回転と、ＣＮＯＴゲート５３０ａ、５３０ｂ、および５３０ｃとを有する。

言い換えれば、回路のこの位相多項式表現は、２つの回転、この場合Ｒ_Ｚ（θ_１）とＲ_Ｚ（θ_４）が、たとえ、これらがこの回路の異なる部分に出現していても、入力の同じアフィン関数に適用される場合を示している。次に、これらの回転は単一の回転に組み合せることができ、回路を改善する。この特定の例において、前述の位置転換技法を用いて、単純化を代替的に達成することができたが、但し、これは一般的には適用できない。関連するアフィン関数が出現する回路中のいかなるポイントにも組み合せ回転を配置するフレキシビリティがある。具体的には、組み合せ回転は、最初の（最左方の）かかる回転の場所に配置される。

サブルーチン４を適用するためには、サブ回路が｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝ゲートだけから成ることを識別する必要がある。このサブ回路は、所定のＣＮＯＴゲートから開始して、一回に１キュービットで構築される。このゲートの最初のキュービットに対し、ＤＡＧ表現を用いて、このキュービットに作用する全ての先行するおよび後続するＮＯＴ、ＣＮＯＴ、およびＲ_Ｚゲートを通してスキャンが行われ、それらがサブ回路に加えられる。回路の開始部または終了部で、アダマールゲートに遭遇した場合、その端部点はマークされ、（各キュービットが１つの開始端部点および１つの終了端部点を有するように）その方向への探索は停止される。このキュービットと、まだ遭遇していない何らかのキュービットとの間の各ＣＮＯＴゲートに対し、そのゲートが新しく遭遇したキュービットに作用するアンカーポイントがマークされる。このプロセスは、最初のＣＮＯＴゲートによって作用された第二のキュービットに対して遂行され、あらゆるアンカーポイントから開始され、新規のキュービットとの遭遇がなくなるまで、プロセスが繰り返される。

もたらされたサブ回路は、ＮＯＴ、ＣＮＯＴ、およびＲ_Ｚゲートだけから成るが、これは多項式位相表現を有さないことがあり、具体的には、サブ回路から出て、再び入る、ワイヤ上の中間アダマールゲートがこれを防ぐことができる。位相多項式形式を適用するために、以下のプルーニング手続を使って、これが生じないようにすることが必要となる。所定の最初のＣＮＯＴゲートから開始して、ネットリスト中のその前後両方のゲートが、端部点に遭遇するまで連続的に検討される。なお、このプロセスで到達したあらゆるＮＯＴおよびＲ_Ｚゲートは、位相多項式が適用されるのを防ぐことがなく、含めることが可能なので、ＣＮＯＴゲートだけを検討する必要がある。遭遇したＣＮＯＴゲートの制御および標的キュービットの両方が端部境界内にあれば、この手続を継続してよい。制御キュービットが端部境界の外側にあるが、標的キュービットは内側にある場合、検査対象のＣＮＯＴゲートが境界の外側になり、当該サブ回路から、該ゲートとその標的キュービットに作用している一切の後続ゲートとが除外されるように、その標的キュービットの端部点が移動される。また一方、制御が境界の内側にあり、標的が外側にある場合、例外処理が行われ、（ＣＮＯＴゲートが当該サブ回路に含まれてはいないが）端部点は移動されない。この例外処理は、図６の図表６００中に示されるように、位相多項式表現に適ったより大きな｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝のサブ回路をもたらす。

図表６００中の回路に特定された例において、上段の（ｑ_１）および中段の（ｑ_２）キュービットに作用している第一のＣＮＯＴゲート（例えば、ＣＮＯＴゲート６３０ａ）から探索を開始することが可能である。ｑ_１から左に進むと、Ｈゲート（例えば、６２０ａ）が見出され、ここに端部点がマークされる。ｑ_１から右に進むと、２つのＣＮＯＴゲート（例えば、ＣＮＯＴゲート６３０ｂ、６３０ｃ）が見出され、１つのＲ_Ｚゲート（例えば、Ｒ_Ｚゲート６１０ａ）が見出され、次いで、Ｈゲート（例えば、Ｈゲート６２０ｂ）が見出され、ここに端部点がマークされる。遭遇したどのＣＮＯＴゲートも、ｑ_１またはｑ_２を残りのキュービットｑ_３に連結していない。同じ手順が、当初のＣＮＯＴゲート（例えば、ＣＮＯＴゲート６３０ａ）からｑ_２上で繰り返される。左側にＲ_Ｚゲートが（例えば、Ｒ_Ｚゲート６１０ｂ）、次いでＨゲート（例えば、Ｈゲート６２０ｃ）が見出され、ここに端部点がマークされる。右に進むと、ｑ_２およびｑ_３に作用するＣＮＯＴゲート（例えば、ＣＮＯＴゲート６３０ｄ）が見出される。このＣＮＯＴゲートは、追加の接続部を示しており、したがって、アンカーポイントは、このＣＮＯＴゲートのｑ_３（標的）側上のｑ_３にマークされる。ｑ_２ワイヤ上のさらに右方に、いずれも追加の接続部を示していない３つのＣＮＯＴゲート（例えば、ＣＮＯＴゲート６３０ｂ、６３０ｅ、および６３０ｃ）と、Ｒ_Ｚゲート（例えば、Ｒ_Ｚゲート６１０ｃ）と、最後にＨゲート（例えば、Ｈゲート６２０ｄ）があり、ここで、端部点がマークされる。次のｑ_３は、前述のアンカーポイントから開始して調査される。左に向かうと、他のキュービットへのさらなる接続のないＨゲート（例えば、Ｈゲート６２０ｅ）が見出され、端部点がマークされる。右に向かうと、Ｈゲート（例えば、Ｈゲート６２０ｆ）が直近に見出され、端部点がマークされる。

サブ回路が構築されると、次に、ネットリスト表現を調べ、それをプルーニングする（ｐｒｕｎｅ ｉｔ）。この経路において、ｑ_２およびｑ_３上で作用している第四のＣＮＯＴゲートに遭遇し（例えば、ＣＮＯＴゲート６３０ｅ）、この制御はサブ回路の境界内にあるが、標的は境界外である。この場合、このプルーニング手続において説明した例外扱いが適用される。これは、当該｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝領域中の最後のＣＮＯＴゲートは含まれるが、（図表６００中に点線の境界で示されるように）第四のＣＮＯＴゲートは除外されることを確実にする。しかして、この回路中に表されている最後のＲ_Ｚゲートは、ｑ_２ライン上の回路の初端、最左端のＨの右に再配置することが可能で、後記でさらに詳しく説明するように、Ｒ_Ｚの併合に基づく位相多項式を可能にする。

有効な｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝サブ回路が識別されたならば、該サブ回路の位相多項式を生成することができる。各Ｒ_Ｚゲートに対し、その位相が適用される関連アフィン関数、および各ゲートに適用される回路の位置を決定することが可能である。次いで、記録されたアフィン関数のリストが分類される。最後に、全てのＲ_Ｚゲートの繰り返しを見出して併合し、併合されたＲ_Ｚを、所望の線形関数を計算するサブ回路中の最初の位置に配置することが可能である。

この手順は、Ｏ（ｇ）のサブ回路を考慮しており、それらの各々を処理するコストが、複雑性Ｏ（ｇ ｌｏｇ ｇ）による分類作業によって占められ、サブルーチン４に対し、Ｏ（ｇ２ｌｏｇ ｇ）の全体的複雑性を課している。しかしながら、実践上、これらのサブ回路は、通常、対象となるサブ回路が多いほど小さくなり、したがって、実際の複雑性はもっと低い。さらに、｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝サブ回路を識別する際に、そのプロセスは、通常、以前に識別されたどの｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝サブ回路にもまだ含まれてないＣＮＯＴゲートから開始され、したがって、サブ回路の数は、実践上はｇよりもずっと少なくできる。望ましい場合、サブ回路の最大サイズを制限することによって、全体的な複雑性をＯ（ｇ）に低下させることが可能である。

最後のステップとして、位相の全てのアフィン関数を線形関数に低減することが可能である。これは、以下：

のように、ＣＮＯＴおよびＲ_Ｚゲートを通り抜けるＮＯＴの伝搬を使って達成することができる。

この手順の適用は、各アフィン関数

が、対応する線形関数

に変形されることを確実にし、これにより、さらなる位相衝突が引き起こされる尤度を改善する。

（サブルーチン５：浮動Ｒ_Ｚゲート）
サブルーチン４では、Ｒ_Ｚゲートに関連するアフィン関数を追跡した。さらに一般的には、Ｒ_Ｚゲートの存在の如何に関わらず、サブ回路中に生じた全てのアフィン関数、およびそれらそれぞれの位置が記録された。しかして、Ｒ_Ｚゲートが回路中に既に出現している位置でなく、Ｒ_Ｚゲートを配置することができる全ての可能な位置を識別することが可能である。この「浮動」Ｒ_Ｚゲート配置の構図において、３つの最適化サブ－サブルーチンを用いることができる。その３つとは、２キュービットゲート相殺、ゲートカウント保存書換えルール、およびゲートカウント低減書換えルールである。

これらサブ－サブルーチンの一番目は、これではＲ_Ｚゲートが浮動であり、焦点が特定の識別されるサブ回路に置かれていることを除き、基本的にはサブルーチン４と同じである。このアプローチは、Ｒ_Ｚゲートを、途上の全ての可能なＲ_Ｚゲートの位置を把握することによって、相殺を促進するように諸Ｒ_Ｚゲートを配置することを可能にする。特に、特定の位置にＲ_Ｚゲートを配置しないと、２つのＣＮＯＴゲートを相殺できる場合、その位置はＲ_Ｚゲートの可能な位置のリストから除去され、その代わりに、ＣＮＯＴ相殺が実行される。

次のゲートカウントを保存する書換えルール（例えば、図７の図表７００を参照）が、さらなる最適化を見出す試みとして適用される。例えば、図表７００は、ゲートカウント保存ルール７１０ａおよび７１０ｂを含む。これらの置き換えは、ゲートを削除しないが、他の箇所での最適化が可能になるように当該回路を修改する。これらの書換えルールは、外部のライブラリファイルによって提供され、サブ回路は、ＤＡＧ表現を用いて、それらに該ルールが適用できるように識別される。これらの置き換えは、それが、前述の浮動Ｒ_Ｚゲートによる２キュービット相殺サブルーチンのもう一回のラウンドを通しての２キュービットゲートカウントの低減につながる場合にだけ適用される。なお、これらの書換えルールは、回路中の特定の位置での特定の浮動Ｒ_Ｚゲートにだけ適用が可能である。このサブルーチンは、浮動Ｒ_Ｚゲートを用いて、ゲートカウントの低減につながるＲ_Ｚゲート位置のこれらの組み合せを選択する。

最後のサブ－サブルーチンは、ゲートカウントを低減する書換えルールを適用する（サブルーチン５において用いられるゲートカウント低減書換えルールについて、例えば、図８の図表８００を参照のこと）。図表８００は、ゲートカウント低減書換えルール８１０ａ、８１０ｂ、８１０ｃ、８１０ｄ、および８１０ｅを含む。これらのルールも、外部のライブラリファイルを介して提供されてよい。これらのルールはそれら自体がゲートカウントを低減するので、この書換えは、一般に、適したパターンが見出されたときはいつでも、実施される。この３ステップのサブルーチン（すなわち、サブルーチン５中の３つの最適化サブ－サブルーチン）の複雑性は、サブ回路の数がＯ（ｇ）であり、各サブ回路内で、２キュービット相殺（サブルーチン３）がＯ（ｇ^２）の複雑性を有するので、上限がＯ（ｇ^３）である。これら書換えルール（例えば、図７中のゲートカウント保存書換えルール、図８中のゲートカウント低減書換えルール）は、複雑性Ｏ（ｇ）で適用することが可能である。なぜなら、サブルーチン１におけるように、回路中のゲート群を通る単一回のパスで十分だからである。ここでも、実践上、サブ回路の数とサブ回路のサイズとは、通常、逆関係にあり、このことは、観測される複雑性を約ｇ分の一に低下させる。この複雑性は、サブ回路の最大サイズを制限することにより、Ｏ（ｇ^２）に低下させることが可能である。実際上、この複雑性は、２キュービットゲートの相殺におけるサブ回路Ａの最大サイズを制限することによって、Ｏ（ｇ ｌｏｇ ｇ）にさらに低下させることができる（分類作業もなお複雑性Ｏ（ｇ ｌｏｇ ｇ）を有することになろう）。

この最適化がどのように作用するかの一例を説明するために、図５の図表５００中の右手の回路を検討する。Ｒ_Ｚ（θ_２）（例えば、ゲート５１０ａ）が、上段キュービットの回路の最終部で実行し、最初の２つのＣＮＯＴの相殺を可能にし、図９の図表９００中の右手側に示された最適化された回路に導くことができるのを観察するが、これは、３つのＲ_Ｚゲート９１０ａ、９１０ｂ、および９１０ｃと、単一のＣＮＯＴゲート９３０ａとを含むことによって、さらになお単純化されている。

（汎用最適化アルゴリズム）
本開示中に記載された最適化アルゴリズムまたは技法は、前述した種々のサブルーチンを注意深く選択された順序で単に適用する。オプティマイザまたは最適化アルゴリズムの２つのバージョンがあり、それらは、軽度バージョン（または単に軽度）および重度バージョン（または単に重度）と言われる。一般に、重度バージョンは、より多くのサブルーチンを適用し、ランタイムがより長いことを犠牲にして、より良好な最適化結果を生成する。前処理ステップは、オプティマイザの軽度および重度バージョン両方において、用いることができる。

オプティマイザの軽度バージョンは、次のシーケンスまたは順序で、最適化サブルーチンを適用する。
１、３、２、３、１、２、４、３、２。

続いて、このシーケンスはさらなる最適化の達成がなくなるまで繰り返される。このシーケンスは、最初にアダマールゲートを削減しながら｛ＣＮＯＴ、Ｒ_Ｚ｝ゲートを露出するという原則に少なくとも部分的に基づき（サブルーチン１）、相殺ルーチンでより大きな低減を可能にし（サブルーチン３、２、３）、特に、単一キュービットゲートの低減を促進するため２つのＣＮＯＴゲートを解放し、またその逆を行う。最初の４つの最適化サブルーチンの後、置き換えルール（サブルーチン１）を適用することによって、さらなる低減を可能にすることができる。次いで、さらなる単一キュービット型の相殺および併合が探索される（サブルーチン２）。これは、さらなるＲ_Ｚカウントの最適化を探すための｛ＮＯＴ、ＣＮＯＴ、Ｒ_Ｚ｝サブ回路領域のより迅速な識別を可能にし（サブルーチン４）、その後、本プロセスは、ゲートの残存する相殺をチェックする（サブルーチン３、２）。

オプティマイザの重度バージョンは、次のシーケンスまたは順序を適用する。
１、３、２、３、１、２、５。

同様に、このシーケンスもさらなる最適化の達成がなくなるまで繰り返される。オプティマイザの重度バージョンのこの最適化シーケンスの最初の６つのステップは、オプティマイザの軽度バージョンのものと同じである。違いは、オプティマイザの重度バージョンでは、用いられるアプローチが浮動Ｒ_Ｚゲートを利用し（サブルーチン５）、さらなるゲート相殺を露出するためのゲートカウント保存書換えルールおよび何らかの残存する非効率性を除去するためのゲートカウント低減書換えルールの使用を含め、ＣＮＯＴゲートのさらなる低減を可能ならしめる、Ｒ_Ｚゲートに対する位置の発見を可能にすることである。

上記で説明した諸順序またはシーケンスは、例示のために提示されたものであり、前述よりも少ないまたは多いサブルーチンが関与する他の順序またはシーケンスも同様に使用することが可能である。

なお、回路表現の変換によって被る計算オーバーヘッドは些少で済むことが多い。全ての変換は、回路のサイズに比例する時間で行うことができる（詳細については、例えば、量子回路の表現と題された節を参照）。諸表現は、必要に応じた程度に整合性を保てばよい。サブルーチン１～サブルーチン３において、迅速に低減箇所を見出すために、ＤＡＧ表現を使って、個別のゲートにアクセスすることができる。これによって、最適化プロセスを続ける前に、ゲートカウントの低減を記録するために、ＤＡＧ表現だけを更新することが可能になる。サブルーチン４およびサブルーチン５において、それら両方の表現は、低減が見出されたときはいつでも、ＤＡＧおよびネットリスト表現両方の整合性を維持しながら、同時に臨機応変に更新することができる。これは、両方のルーチンが、位相多項式表現を使う低減に適ったサブ回路を識別するので有用である。この識別プロセスは、最新のＤＡＧ表現を必要とし、位相多項式表現の生成には、最新のネットリスト表現が必要である。なお、位相多項式表現は、識別されたサブ回路中の最適化を支援するためだけに用いられ、位相多項式表現をネットリストまたはＤＡＧのいずれかに戻して変換する必要はない。しかして、この位相多項式表現は、対応するサブ回路最適化プロセスが完了したときに、安全に消去されて（ｐｕｒｇｅｄ）よい。

（特殊用途最適化）
前述した汎用的オプティマイザ（例えば、汎用最適化アルゴリズム）に加え、特定の構造を備えた回路を改良するために、２つの特化型最適化を用いることができる。

ＬＣＲオプティマイザ：積公式シミュレーションアルゴリズムなどの一部の量子アルゴリズムは、一定のブロックの複数回の繰り返しを含む。かかる回路を最適化するためには、最初に、単一のブロックにオプティマイザを実行して、その最適化バージョン、Ｏを得るのがベストであり得る。複数のブロックに亘っての単純化を見出すために、回路Ｏ^２を最適化することが可能で、その結果はＬＲと呼ぶことができ、このＬは、Ｏ^２の最適化中のＯの最大プレフィックスである。次のステップは、Ｏ^３の最適化である。最適化がブロック間の境界の近辺で行われるという条件の下で、Ｏ^３の最適化バージョンからプレフィックスＬおよびサフィックスＲを除去することが可能で、残った回路をＣと呼ぶ。かかるＬ、ＣおよびＲが発見可能であるとすれば（実践上ではあり得る）、Ｏ^ｔをＬＣ^ｔ－２Ｒに単純化することが可能である。

トフォリ分解：多くの量子アルゴリズムは、トフォリゲートを使って自然に説明される。本開示に記載されたオプティマイザまたは最適化アルゴリズムは、正および負の両方の制御でトフォリゲートを取り扱うことが可能である。目標は、ゲートのセット｛ＮＯＴ、ＣＮＯＴ、Ｈ、Ｒ_Ｚ｝に関する回路を表現することなので、トフォリゲートは、これらの基本的なゲートに関連させて分解される。本明細書に記載の技法は、最適化の品質を向上させるために、これを行う各種のやり方を利用する。

具体的に、トフォリゲートは、図１０の図表１０００中に示された特徴を用いる１キュービットゲートおよび２キュービットゲートに関連させて拡大することが可能で、それらの回路分解において、（トフォリゲートは自己反転なので）ＴとＴ^†とを交換することによって、所望のトフォリゲートを得ることも可能であることを念頭に置く。前述したように、図表１０００中の特徴は、トフォリゲート１０１０ａ、上部に否定制御を備えたトフォリゲート１０１０ｂ、および両方が否定制御であるトフォリゲート１０１０ｃを含む。オプティマイザは、当初、各トフォリ分解において、Ｔ／Ｔ^†の極性（すなわち、どちらのゲートがダガーを含みどちらが含まないかの選択）を未定のまま残す。オプティマイザは、未定のＴおよびＴ^†ゲートを、それらの相対的極性を保持しながら、所与の量子回路内で単にそれらの位置を移動することによって記号的に処理する。この最適化は、これら未定のＴおよびＴ^†ゲートの移動によってもはやゲートカウントを低減できなくなった時に、完了すると見なされる。最後に、最適化された回路中のＴのカウントを最小化することを目標に、各トフォリゲートの極性を（それらの間の固定された関係を条件として）選択することが可能である。この最小化は、最長マッチを行う方法（ｇｒｅｅｄｙ ｗａｙ）で行われ、ほぼ最適化された回路中の関連するＴ／Ｔ^†ゲートの出現の順序で、Ｔのカウントができるだけ多く低減されるように、各トフォリゲートに対する極性を選択する。

全体的には、この極性の選択プロセスはＯ（ｇ）の時間を要する。極性を選択した後、サブルーチン３およびサブルーチン２を実行することが可能である。というのは、特定の極性の選択が、該選択がなければ、所望の位置転換を阻止していた未定のゲートの存在に起因してできなかった、ＣＮＯＴゲートおよび単一キュービットゲートのさらなる相殺につながり得るからである。

（最適化の結果）
本明細書に記載のオプティマイザまたは最適化アルゴリズムの諸態様のインプリメンテーション例は、フォートランプログラミング言語を用いて作成され、３セットのベンチマーク回路を用いてテストされた。本開示中に示された全ての結果は、２．９ＧＨｚインテルコアｉ５プロセッサ、および８ＧＢの１８６７ＭＨｚ ＤＤＤＲ４メモリを備えたマシンを使い、ＯＳ Ｘ ＥＩ Ｃａｐｉｔａｎを実行して得られた。当然のことながら、これらの結果は、例示目的で、相対的なパフォーマンスを示すために提示され、また、これらの結果は、異なるハードウェア構成を使って行われた場合に、変わり得る。

テストの一部として、ショアの整数因数分解アルゴリズムのコンポーネントを含む様々な量子回路、具体的には量子フーリエ変換（ＱＦＴ：ｑｕａｎｔｕｍ Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ）および整数加算回路を検討した。また、ハミルトニアンシミュレーションへの積公式（ＰＦ：ｐｒｏｄｕｃｔ ｆｏｒｍｕｌａ）アプローチのための量子回路も検討した。両方のケースにおいて、テストの焦点は、従来式の計算を凌駕する、諸アプリケーションにおいて有用な可能性のあるサイズの量子回路に置かれ、各種の型の加算回路および積公式による実験が実行された。最後に、（ガロア体乗算回路のファミリを含め）様々な演算回路、および多重制御トフォリゲートのインプリメンテーションから成るベンチマーク量子回路のセットも検討した。比較のため、最適化の前および後のこれらの回路を記載したデータファイルが利用可能である。

オプティマイザまたは最適化アルゴリズムの正確性は、最適化の前と後との様々なテスト回路の機能的等価性（すなわち、対応するユニタリ行列の同等性）を検証することによってチェックされた。かかるテストは、一般に少数のキュービットを用いる回路に対し実行可能である。これらのテストは、図１３の図表１３００中に示された表１および図１４の図表１４００中に示された表２中の全ての８キュービットのベンチマークに対して行われた。表１は、加算回路の軽度最適化を、上段にＱＦＴ（またはＱＦＡ）に基づくインプレース加算回路に対する結果、および下段にクイッパーライブラリ加算回路に対する結果とともに含む。また、これらのテストは、それぞれ図１６Ａおよび１６Ｂの図表１６００ａおよび１６００ｂ中に示された表３Ａおよび表３Ｂ中の全ての１０キュービットのベンチマークに対しても行われた。表３Ａは、ＣＮＯＴゲートの低減による積公式アルゴリズムの最適化を示し、表３Ｂは、Ｒ_Ｚゲートカウントの低減による積公式アルゴリズムの最適化を示す。ソフトウェアのランタイム範囲は、０．００４秒（第一オーダー、ｎ＝１０）～０．１３７秒（第六オーダー、ｎ＝１００）である。クリフォードゲート低減は、（ｎにかかわらず、第一オーダー公式に対し）アダマールについては６２．５％、位相ゲートについては７５％から、（これもｎにかかわらず、第六オーダー公式に対し）アダマールについては７５％、位相ゲートについては８５％までの範囲である。記号「（×１０００）」は、第一オーダー公式に対するゲートカウントが、（丸め誤差なしの）千の単位であることを示す。記号「（Ｌ）」は、前述した最適化の標準的な軽度バージョンを示す。

さらに、図１７Ａ～１７Ｃの図表１７００ａ～１７００ｃ中に示された表４Ａ～４Ｃは、以下のベンチマーク、Ｍｏｄ５_４、ＶＢＥ―Ａｄｄｅｒ_３、ＣＳＬＡ－ＭＵＸ_３、ＲＣ－Ａｄｄｅｒ_６、Ｍｏｄ－Ｒｅｄ_２１、Ｍｏｄ－Ｍｕｌｔ_５５、Ｔｏｆｆ－Ｂａｒｅｎｃｏ_３．．５、Ｔｏｆｆ－ＮＣ_３．．５、ＧＦ（２^４）－Ｍｕｌｔ、およびＧＦ（２^５）－Ｍｕｌｔと比較した、本オプティマイザの結果を表す。例えば、表４Ａ～４Ｃは、Ｔ－ｐａｒとの比較（後記でさらに詳しく説明する）を含み、これらベンチマークのアルゴリズムの名称は、Ｔｏｆｆ－ＢａｒｅｎｃｏおよびＴｏｆｆ－ＮＣが乗算制御のトフォリゲートのインプリメンテーションを示すのに使われていることを除き、周知のものである。記号「（Ｌ）」は、本最適化の標準的な軽度バージョンを示し、一方、記号「（Ｈ）」は、本最適化の標準的な重度バージョンを示す。記号－｜｜－は、最適化の軽度バージョンに対し、最適化の重度バージョンに改善がないことを示す。

（ＱＦＴおよび加算回路）
ＱＦＴは、量子計算における基本的サブルーチンであり、指数関数的高速化のため多くの量子アルゴリズムに出現する。正確なｎキュービットＱＦＴに対する標準的回路は、Ｒ_Ｚゲートを用い、その一部はｎが指数関数的に小さい角度を有する。非常に小さな回転角度を有するゲートを省くことによって、高精度の近似ＱＦＴが行えることがよく知られている。本開示のオプティマイザまたは最適化アルゴリズムの目的のためには、最大π／２^１３までの角度の回転は省いてもよく、これは、対象以下のサイズの量子回路に対しては十分な精度の近似ＱＦＴを確実にする。これらの小さな回転は、それらの省略が結果中に報告される改善に寄与しないので、最適化の前に除去される。

図１１の図表１１００中に示されるチャートは、近似量子フーリエ変換（ＱＦＴ、挿入図）、クイッパーライブラリ加算回路、およびフーリエベースの加算回路（ＱＦＡ）に対するゲートカウントを示す。白／黒の記号は、最適化の前／後のゲートカウントを表し、四角／丸／三角の記号は、それぞれ、クイッパーライブラリ加算回路／ＱＦＡ／ＱＦＴに対するゲートカウントを表す。５１２以上のキュービットを備えるＱＦＴに対して、３６％を超える節減レートが観察される。この最適化は、もっぱら、フォールトトレラントなインプリメンテーションで最も高価なリソースであるＲ_Ｚゲートの数の低減から来ている。

図１１の図表１１００中のチャートに示されるように、２つの型の整数加算回路、すなわち、クイッパーライブラリ中にインプリメントされているようなインプレースモジュロ（ｉｎ－ｐｌａｃｅ ｍｏｄｕｌｏ）２^ｑ加算回路とＱＦＴ（またはＱＦＡ）に基づくインプレース加算回路とを検討した。ＱＦＡ回路は、前述したように、π／２^１３より小さい角度による回転が除去された近似ＱＦＴを使用する。加算回路は、整数因数分解に対するショアの量子アルゴリズム（Ｓｈｏｒ’ｓ ｑｕａｎｔｕｍ ａｌｇｏｒｉｔｈｍ）の基本的コンポーネントである。示された結果は、クイッパー加算回路、およびＬが４～１１の範囲にある２Ｌキュービットに作用する回路のＱＦＡに対する、最適化の前と後とのゲートカウントを報告している。例えば、１，０２４ビットの数を因数分解するためのショアの整数因数分解アルゴリズムでは、Ｌ＝１０の加算回路が使われる。関連するＲＳＡ－１０２４への課題は、まだ解決されていない。

加算回路の軽度（バージョン）の最適化の結果が図１１および図１３に示されている。クイッパーライブラリ加算回路に対しては、フルライトオプティマイザ（ｆｕｌｌ Ｌｉｇｈｔ ｏｐｔｉｍｉｚｅｒ）が使用された。ＱＦＡ最適化に対しては、標準的なライトオプティマイザの最後の３つのサブルーチン４、３、２を省いてサブルーチンのシーケンス１、３、２、３、１、２を有する変更されたライトオプティマイザ（ｍｏｄｉｆｉｅｄ Ｌｉｇｈｔ ｏｐｔｉｍｉｚｅｒ）が代わりに使用された。これは、小さなインスタンス（ｎ≦２５６）で、それらのサブルーチンからさらなるゲート節減が見出されなかった際に行われた。

簡素化されたクイッパーライブラリ加算回路は、大幅な差でＱＦＡを凌駕し、これが実践上選好される可能性があることを暗示している。このクイッパーライブラリ加算回路に対し、最大５．２分の１までのＴゲートのカウントの低減が観察された。この低減は、回路構造のいかなる事前知識も用いることなく、もっぱら自動化手段によって（すなわち、最適化アルゴリズムによって）得られた。ショアの整数因数分解アルゴリズムは冪剰余のコストに左右され、次いで、これは、主として整数の加算に依拠しているので、この最適化は、全体的因数分解アルゴリズムを実行するコストを５分の１未満に低減する。

また、重度オプティマイザも、ＱＦＴおよび加算回路に適用された。ＱＦＴおよびＱＦＡ回路に対して、重度の設定はゲートカウントを改善しなかった。クイッパー加算回路に対する重度最適化の結果が、図１４の図表１４００中の表２に示されている。ＣＮＯＴカウントの低減は、２．７分の１であることが分かったが、これと比較すると、軽度最適化での低減は、わずかに１．７分の１であった。図１２の図表１２００中のチャートは、最適化の前（例えば、プレ最適化）、軽度最適化の後、および重度最適化の後のクイッパーライブラリ加算回路の合計ＣＮＯＴカウントを表し、これら２つの型の最適化によるＣＮＯＴカウントの低減を示している。白／ドットパターン／黒の四角の記号は、それぞれ、プレ最適化／ポスト軽度最適化／ポスト重度最適化でのゲートカウントを表す。
量子シミュレーション

ハミルトニアン力学（Ｈａｍｉｌｔｏｎｉａｎ ｄｙｎａｍｉｃｓ）をシミュレートするための第一の明示的多項式時間量子アルゴリズムが導入され、そのアプローチが後により高次の積公式（ＰＦ）に一般化され、改善された漸近法計算量をもたらした。本開示は、オーダー１、２、４および６（１より高いオーダーに対して、標準的スズキ積公式の構成のオーダーは偶数である）のＰＦアルゴリズムに対する最適化の前と後とのゲートカウントを報告している。具体的に、これらのアルゴリズムは、周期的境界条件を有する磁界中の一次元のハイゼンベルグモデルに対しインプリメントされ、システムをそのサイズに比例する時間の間展開し、アルゴリズムのパラメータは、積公式近似の誤差に対する既知の限界を用いて、ハミルトニアンシミュレーション誤差が最大でも１０－３となることを保証するように選択された。

積公式アルゴリズムの軽度最適化の結果が（例えば、図１６Ａおよび１６Ｂの）表３Ａおよび表３Ｂに報告され、図１５の図表１５００中のチャートに図示されている。この図の白色／黒色の記号は、それぞれ、最適化の前／後のゲートカウントを表し、四角／丸の記号は、第二／第四オーダーに対するゲートカウントを表す。これらの積公式アルゴリズムに対し、重度最適化は、軽度最適化を上回るさらなる改善を示さなかった。第二の、第四、および第六オーダーアルゴリズムでは、ＣＮＯＴカウントで約３３．３％の低減、Ｒ_Ｚカウントで約２８．５％の低減が認められ、おおむね、物理レベルおよびロジカルレベルのインプリメンテーションに関連する低減に対応している。第一オーダー公式アルゴリズムは、ＣＮＯＴまたはＲ_Ｚゲート最適化を示さなかった。全ての積公式アルゴリズムにおいて、位相ゲートおよびアダマールゲートの数は、おおむね３分の１～６分の１に、大きく低減された。

（他のアプローチとの比較）
量子回路最適化は、周知の分野である。しかしながら、回路の最適化についてのこれまでの研究は、本開示で検討しているような、従来式のコンピュータを凌駕するような型の大規模量子回路を検討してこなかった。例えば、Ａｍｙ、Ｍａｓｌｏｖ、およびＭｏｓｃａによる論文（「Ｐｏｌｙｎｏｍｉａｌ－ｔｉｍｅ Ｔ－ｄｅｐｔｈ ｏｐｔｉｍｉｚａｔｉｏｎ ｏｆ Ｃｌｉｆｆｏｒｄ ＋Ｔ ｃｉｒｃｕｉｔｓ ｖｉａ ｍａｔｒｏｉｄ ｐａｒｔｉｔｉｏｎｉｎｇ」、ＩＥＥＥ Ｔｒａｎｓ．Ｃｏｍｐｕｔ．Ａｉｄｅｄ Ｄｅｓ．Ｉｎｔｅｇｒ．Ｃｉｒｃｕｉｔｓ Ｓｙｓｔ．，ｖｏｌ．３３，１４７６－１４８９頁（２０１４））では、ｇのゲート回路の最適化の複雑性はＯ（ｇ^３）であり、大規模回路の最適化を非現実的だとしている。キュービットの数がｎ＝１０～３５、およびゲートカウントが６０～３６８に対し、実行時間が０．０７～１．８８３秒の範囲との例があるが、本明細書に記載の最適化アルゴリズムは、表１（図１３）に示すように、上限ｎ＝２５６までで、約２３，０００個のゲートのクイッパー加算回路を最適化する際に、同じ位の時間で実行した。また、最適ゲートライブラリを使ったピープホール（ｐｅｅｐ ｈｏｌｅ）最適化に基づく、量子回路最適化の他の例もあるが、これは高コストで、場合によっては、２０キュービット、１，０００個のゲートの回路に対し１００秒より長くかかる。

本開示で提案する最適化アルゴリズムを用いて生成された結果を、従前に報告されている結果と比べるために、ＴおよびＣＮＯＴカウントの重み付き組み合せを検討した。Ｔゲートは、状態蒸留（ｓｔａｔｅ ｄｉｓｔｉｌｌａｔｉｏｎ）を使ってフォールトトレラントにインプリメントするには相当に高価とはなり得るが、ＣＮＯＴゲートのコストを無視するのは、著しい過小見積もりにつながる可能性がある。解析では、フォールトトレラントなＴゲートをインプリメントするのに、局部的にフォールトトレラントなＣＮＯＴゲートよりも４６～３５０倍高価となり得ることを示唆しており、そのコスト比率が１：５０であることに関し、少なくとも１つの推奨をしている。実際のオーバーヘッドは、フォールトトレランスのスキーム、誤差モデル、計算のサイズ、アーキテクチャ上の制限、Ｔゲートのインプリメンテーションを最適化できる程度、およびＴ状態産生が、そのコストが（部分的に）割引できるようにオフラインで行えるかどうかを含め、多くの細目により決まる。おおまかな比較のための１つのアプローチは、＃Ｔ＋０．０１・ｌｏｇ ｎ・＃ＣＮＯＴで定義される総計コスト指標を目途に作業をすることであり、この式の＃Ｔは使われるＴゲートの数であり、０．０１はＴゲートに対するＣＮＯＴゲートの相対効率割合であり、ｎは計算中のキュービットの数であり、＃ＣＮＯＴは使われるＣＮＯＴゲートの数である。この式で、ｌｏｇ ｎの因数は、現実的アーキテクチャ中のキュービットの間のゲートを実行するための一般的なコストを過小評価している（実際のコストは、三次元では^３√ｎに、または二次元では√ｎにより近くなり得る）。本明細書に記載のアプローチは、当初の回路の構造を保持しているので、この指標は、ロングレンジゲートを導入している可能性のある（後記で述べるＴ－ｐａｒアプローチなどの）他のアプローチに対し控え目の比較を提示しているはずである。したがって、この総計コストに対する有利性を示すことで、本開示に記載の本提案の最適化技法のメリットを、おおまかに実証することができる。

本開示に記載の技法を用いた結果を、マトロイド分割（ｍａｔｒｏｉｄ ｐａｒｔｉｔｉｏｎｉｎｇ）に基づく技法を用いてＴのカウントおよびＴの深さを低減することを狙った、Ａｍｙ、Ｍａｓｌｏｖ、およびＭｏｓｃａにより報告された結果と、直接に比較することが可能である。そのアプローチは、本開示ではＴ－ｐａｒと称しているものである。本開示に記載の技法の使用による結果が、上記の論文に現れたベンチマークのセットを最適化するのに使われ、表４Ａ～４Ｃ（図１７Ａ～１７Ｃ）に示されるように、その結果をＴ－ｐａｒの最適化と比較した。

これらのベンチマーク回路は、３つのカテゴリに分類される。第一のセット（図１７Ａ中の表４Ａ）は、算術演算の選抜から成る。これらの回路に対し、ベンチマークに比べて、より良いまたはこれに合致するＴのカウントが得られ、同時にかなり良いＣＮＯＴのカウントも得られた。なお、回路ＣＳＬＡ－ＭＵＸ３は、Ｔ－ｐａｒがそれを不正確に最適化すると考えられているので、比較から除かれた。本明細書に記載の最適化技法の利点を例証するため、前述の総計コスト指標を用いて、本提案の技法はＲＣ－Ａｄｄｅｒ６回路のコストを７１．９１から４９．７０に低減した。コスト面での改善は、しかして、ほとんどが低減されたＴゲートのカウントに起因しており約３１％であった。

ベンチマークの第二のセットは、乗算制御のトフォリゲートから成る（図１７Ｂ中の表４Ｂ）。本提案のオプティマイザは、Ｔ－ｐａｒオプティマイザによって得られたＴカウントと同等であり、ＣＮＯＴカウントを著しく低減したが、本提案のオプティマイザも従前のアプローチも最もよく知られたインプリメンテーションを見出すことはできなかった。この最もよく知られたインプリメンテーション中では非常に異なった回路構造が用いられていることを考えれば、これは驚くべきことではない。

ベンチマークの第三のセットは、ガロア体乗算回路を包含する（図１７Ｃ中の表４Ｃ）。一例では、重度オプティマイザは、そのランタイムが軽度オプティマイザのランタイムの２００倍を超えた時点で終了された。かかる時間切れは、ガロア体乗算回路の４つの最大インスタンスに重度オプティマイザを適用したときに生じた。テストされたケースでは、重度オプティマイザに軽度オプティマイザを上回る利点はなく、したがって、該４つの最大インスタンスに重度オプティマイザは適用されなかった（対応するエントリは表４Ｃ中の左の空白部である）。Ｔカウントは、これもＴ－ｐａｒオプティマイザと同等であったが、ＣＮＯＴカウントは、かなり低下し、明らかに好ましい回路をもたらした。例えば、最適化済みＧＦ（２^６４）乗算回路は、１８０，８９２個のＣＮＯＴゲートの使用を示していたが、本開示の技法から得られた最適化インプリメンテーションは、わずか２４，７６５個のＣＮＯＴゲートを使用し、しかして、総計コストは、Ｔカウントの変化がなかったにも関わらず、３０，１６８．５９から１８，３２６．４２に、すなわち、約３９％低減された。この低減は、ほとんど、または全面的に、ＣＮＯＴゲートから来たものである。この比較は、Ｔカウントと理論的に予測された現実的な総計コスト見積もりとの間の不一致が、実践上でも示されることを実証している。軽度オプティマイザの効率性は、未解決で残っている楕円曲線離散対数問題のインスタンスに対応して、ＧＦ（２^１３１）およびＧＦ（２^１６３）乗算量子回路の最適化を可能にした。報告されたＴ－ｐａｒ最適化のランタイムを考えると、このサイズのインスタンスは、Ｔ－ｐａｒオプティマイザには手に負えないように見える。

Ｔカウント最適化に対する新規のツールは、（ＨｅｙｆｒｏｎおよびＣａｍｐｂｅｌｌによって）提案されている。この新規ツールで使われるアプローチは、本開示で考慮されている全面的なユニタリ回路と対照的に、測定と従来式のフィードバックとに依拠している。さらに、この新規ツールは、ＣＮＯＴのカウントを設けておらず、ＴおよびＣＮＯＴゲート両方に対するカウントを行うツールとの直接的比較は不可能である。再度、このツールは、Ｔカウントの最適化だけを目指しており、本開示の技法は、この単純なコスト指標からは外れている。例えば、この新規のＴカウント最適化ツールに基づき最適化されたＱＦＴ４回路は、４４のキュービットを使う４キュービットＱＦＴ変換をインプリメントし、ＣＮＯＴゲートのオーバーヘッドが大きいに違いないことを示唆している。さらなる重要な違いは、スケーラビリティである。本明細書に記載の技法が、大規模回路を最適化するのに適し、適用される一方、該新規のＴカウント最適化ツールは、非常に小さい回路だけ取扱い、例えば、このツールで編成された最大のＧＦ乗算回路は７ビットケースであるが、本明細書に記載のものは、１３１および１６３ビットによるＣＦ乗算に取り組むことが可能であり、未解決のＣｅｒｔｉｃｏｍチャレンジに対応している。別の差異は、本明細書で提案する技法は、入力回路で既に利用可能なキュービット間の相互作用だけを用いることである。これは、同じアーキテクチャ中の最適化済み回路を入力回路として実行することを可能にし、限定されたアーキテクチャよりも有用な量子コンピュータであり得る。対照的に、この新規Ｔカウント最適化ツールは、新しい相互作用を導入している。最後に、本明細書で提案する技法は、任意のＲ_Ｚゲートを介して回路を取り扱えるが、該新規のＴカウント最適化ツールは、クリフォード＋Ｔ回路に限定されている。

（全体的パフォーマンス）
本開示に記載の最適化アルゴリズムにより生成された数値的な最適化結果の諸例が、表１、表２、表３Ａおよび３Ｂ、ならびに表４Ａ～４Ｃに提示されている。これらの表は、従来式のコンピュータの到達域を超える実際上の量子計算に関連するベンチマークを包含している。表１および表２中には、ショアの整数因数分解アルゴリズムの、従来式では手に負えないインスタンス中で使われている１，０２４および２，０４８キュービットＱＦＴおよび整数加算回路がある。表３Ａおよび３Ｂ中には、量子力学の直接的な従来式シミュレーションでは、現在実行不可能な、ｎ≧５０の全インスタンスが含まれる。表４Ａ～４Ｃ中には、サイズ１３１および１６３のバイナリフィールドに対するガロア体乗算回路があり、これらは、未解決のＣｅｒｔｉｃｏｍのＥＣＣチャレンジ問題への量子アタックに関連している。これは、本提案のオプティマイザまたは最適化アルゴリズムが、実務に適する十分に大きな量子回路を取り扱う能力を有することを例証している。

本提案のオプティマイザまたは最適化アルゴリズムは、回路の最適化に対する従前の研究よりも、もっと広範に適用することが可能である。これは、トフォリゲート（これは否定制御を有し得る）などの複合ゲートを容易に受け入れる。また、これは、ハミルトニアンシミュレーションおよび因数分解を含め、必然的にＲ_Ｚゲートを使用するアルゴリズムに対する有用な機能である連続するパラメータを備えたゲートを取り扱う。トラップイオンおよび超電導回路の両方を含め、多くの量子情報処理技術は本来的にかかるゲートをサポートしており、しかして、本明細書に記載の提案のアプローチは、物理レベルの回路を最適化するため有用であり得る。

フォールトトレラントな量子計算は、一般に、クリフォード＋Ｔなどのディスクリートのゲートのセットに依拠しており、最適なＲ_Ｚゲートのクリフォード＋Ｔインプリメンテーションが既に知られている。そうであっても、フォールトトレラントな設定において、連続するパラメータによって回路を最適化する能力は、同様に貴重である。その理由は、ディスクリートなフォールトトレラントなセットに編成される前の、連続的にパラメータ化された自然なゲートのセットに対する最適化は、より小さな最終回路をもたらす可能性が高いからである。

最後に、従前のアプローチと違って、本開示で提案するオプティマイザは、当初の量子回路の構造を保持するように構成される。具体的に、最適化された量子回路によって使われる２キュービット相互作用のセットは、当初の量子回路中で使われていたものサブセットである。これは、前処理ステップも、最適化アルゴリズムも、いかなる新規の２キュービットゲートも導入しないために行える。（例えば、使われる相互作用のセットを劇的に増加させるＴ－ｐａｒとは全く対照的に）、使われる相互作用の型を制御下に置くことによって、本提案の最適化アルゴリズムからもたらされる最適化インプリメンテーションは、限定された連結性を備えるアーキテクチャにより良好に適応する。具体的に、限定された連結性を備えるハードウェア上の当初の量子回路のレイアウトを所与として、本プロパティによって、最適化された量子回路に対しても同じレイアウトを用いることが可能になる。さらに、他の最適化ツール（例えば、前述した新規Ｔカウント最適化ツール）と違って、本開示で提案するオプティマイザは、使われるＣＮＯＴゲートの数を増加しない。これは、適切な実務上の考慮であり得る。なぜなら、ロングレンジＣＮＯＴゲートがＴゲートよりも高価なことさえあり得、Ｔの最適化だけに焦点を合わせることで、そのコストがＣＮＯＴゲートによって左右される回路をもたらす可能性があるからである。

ここで図１８の図表１８００を参照すると、一インプリメンテーションによる、例示のコンピュータデバイス１８１０が示されており、コンピュータデバイス１８１０は、様々な最適化サブルーチンとオプティマイザの軽度バージョンおよび重度バージョンを含め、前述した様々な最適化アルゴリズムを実行するために使用が可能である。コンピュータデバイス１８１０は、例えば、単一コンピューティングデバイス、複数コンピューティングデバイス、または分散型コンピューティングシステムを表していてよい。コンピュータデバイス１８１０は、本明細書で説明したように、連続するパラメータを備える大規模量子回路の自動最適化を含め、量子回路の自動最適化を実行または遂行するように構成されてよい。さらに、コンピュータデバイス１８１０は、量子回路についての情報（例えば、ネットリスト）を受信し、量子回路（または、量子回路の機能）をインプリメントするために必要な量子ゲートの数が、該量子回路（または、量子回路の機能）をインプリメントするために当初に（例えば、最適化の前に）必要としていた量子ゲートの数よりも少なくなるように最適化した後、新しい情報（例えば、新規ネットリスト）を生成するように構成することができる。

一例において、コンピュータデバイス１８１０は、１つ以上の最適化機能または本明細書に記載のオペレーションに関連する処理機能を実行するプロセッサ１８４８を含むことが可能である。プロセッサ１８４８は、プロセッサまたはマルチコアプロセッサの単一のまたは複数のセットを含むことができる。さらに、プロセッサ１８４８は、集積処理システムおよび／または分散型処理システムとしてインプリメントされてもよい。一インプリメンテーションにおいて、プロセッサ１８４８は、例えば、中央処理ユニット（ＣＰＵ：ｃｅｎｔｒａｌ ｐｒｏｃｅｓｓｉｎｇ ｕｎｉｔ）、グラフィックス処理ユニット（ＧＰＵ：ｇｒａｐｈｉｃｓ ｐｒｏｃｅｓｓｉｎｇ ｕｎｉｔ）、テンソル処理ユニット（ＴＰＵ：ｔｅｎｓｏｒ ｐｒｏｃｅｓｓｉｎｇ ｕｎｉｔ）またはこれらの種類の１つ以上のユニットの組み合せを含んでよく、それらを本明細書に記載の最適化機能またはオペレーションの１つ以上を行うように構成することができる。

一例において、コンピュータデバイス１８１０は、本明細書に記載の諸機能を実行するために、プロセッサ１８４８によって実行可能な命令を格納するメモリ１８５０を含むことができる。一インプリメンテーションにおいて、メモリ１８５０は、例えば、本明細書に記載の最適化機能またはオペレーションの１つ以上を実行するためのコードまたは命令を格納するコンピュータ可読ストレージ媒体に対応していてよい。

さらに、コンピュータデバイス１８１０は、本明細書に記載したようなハードウェア、ソフトウェア、およびサービスを用いる１つ以上のパーティとの通信の確立および維持を提供する通信コンポーネント１８５２含むことができる。通信コンポーネント１８５２は、コンピュータデバイス１８１０上の諸コンポーネントの間、および、コンピュータデバイス１８１０と、通信ネットワークに亘って配置された諸デバイスおよび／またはコンピュータデバイス１８１０にシリアルにまたはローカルに連結された諸デバイスのような外部デバイスとの間で通信を搬送することが可能である。一例において、通信コンポーネント１８５２は、コンピュータデバイス１８１０と、量子情報処理（ＱＩＰ）システムとの間で、最適化されたネットリストまたは類似情報がコンピュータデバイス１８１０によって生成された後それがＱＩＰシステムに供給されるように、情報の通信を提供することができる。別の例において、通信コンポーネント１８５２は、１つ以上のバスを含むことができ、送信機と受信機とにそれぞれ関連付けられ、外部デバイスをインターフェース接続するよう作動可能な、送信チェーンコンポーネントおよび受信チェーンコンポーネントをさらに含むことが可能である。

さらに、コンピュータデバイス１８１０は、データストア１８５４を含むことができ、これは、本明細書に記載のインプリメンテーションに関連して用いられる情報、データベース、およびプログラムの大容量ストレージを提供するハードウェアおよび／またはソフトウェアの任意の適切な組み合せであってよい。例えば、データストア１８５４は、オペレーティングシステム１８４０および／または最適化アプリケーション１８３０のためのデータリポジトリであってよい。一インプリメンテーションにおいて、データストア１８５４は、メモリ１８５０を含むことができる。

データストア１８５４および／またはメモリ１８５０は、最適化前回路、最適化後回路に関連する情報、最適化の過程で生成された中間情報、サブルーチンおよびオプティマイザの様々なバージョンを含む最適化アルゴリズム、本明細書に記載の最適化オペレーションに関連する任意のルールを格納するために使われてよい。

また、コンピュータデバイス１８１０は、コンピュータデバイス１８１０のユーザからの入力を受信するよう動作可能で、ユーザへのプレゼンテーションのための出力を生成するようにさらに動作可能なユーザインターフェースコンポーネント１８５６を含むことができる。ユーザインターフェースコンポーネント１８５６は、以下に限らないが、キーボード、テンキー、マウス、タッチセンシティブディスプレイ、デジタイザ、ナビゲーションキー、ファンクションキー、マイクロフォン、声認識コンポーネント、ユーザから入力を受信する能力のある任意の他のメカニズム、またはこれらの任意の組み合せを含め、１つ以上の入力デバイスを含むことができる。さらに、ユーザインターフェースコンポーネント１８５６は、以下に限らないが、ディスプレイ、スピーカ、触覚フィードバックメカニズム、プリンタ、ユーザに出力を提示する能力のある任意の他のメカニズム、またはこれらの任意の組み合せを含め、１つ以上の出力デバイスを含むことができる。

一インプリメンテーションにおいて、ユーザインターフェースコンポーネント１８５６は、オペレーティングシステム１８４０および／または最適化アプリケーション１８３０のオペレーションに対応して、メッセージを送信および／または受信することが可能である。さらに、プロセッサ１８４０は、オペレーティングシステム１８４０および／または最適化アプリケーション１８３０を実行することができ、メモリ１８５０またはデータストア１８５４は、それらを格納することができる。さらに、最適化アプリケーション１８３０は、本明細書に記載の最適化機能またはオペレーションの１つ以上を遂行するため、格納されたコードまたは命令に基づいて、実行することが可能である。例えば、最適化アプリケーション１８３０は、オプティマイザの軽度バージョンまたはオプティマイザの重度バージョンの間で選択することができ、これにより、サブルーチンの適切なシーケンス（およびサブルーチンの各々内の任意のオペレーション）を実行することができる。

図１９は、本開示の態様による量子回路の自動最適化を行うためのプロセスまたは方法１９００を示すフローチャートである。方法１９００の諸態様は、図１８に示されたコンピュータデバイス１８１０中のハードウェアおよび／またはソフトウェアによって実行することができる。

ブロック１９０５で、方法１９００は、量子回路を形成する量子ゲートの第一のリストについての情報を含むネットリストを受信するステップを含む。量子ゲートのリストは、量子回路が同じ量子ゲートの重複を有し得るので、諸要素の繰り返しを許容してよい。

ブロック１９１０で、方法１９００は、量子ゲートの第一のリストとの機能等価性を有する量子ゲートの第二のリストを生成するために、量子ゲートの第一のリストについての情報に位相多項式低減オペレーションを実行するステップを含み、量子ゲートの第二のリスト中の量子ゲートの数などの計算コストは、量子ゲートの第一のリスト中の量子ゲートの数等の計算コストよりも小さい。

ブロック１９１５で、方法１９００は、量子ゲートの第二のリストについての情報を含む新規ネットリストを生成するステップを含む。

ブロック１９２０で、方法１９００は、量子ゲートの第二のリストを用いて、量子回路の機能をインプリメントするために、新規ネットリストを提供するステップを含む。

方法１９００の別の態様において、位相多項式低減オペレーションを行う前に、前処理オペレーションを実行することができる。この前処理オペレーションは、ＮＯＴゲート、ＣＮＯＴゲート、トフォリゲート、アダマールゲート、およびＲ_Ｚゲートに適用することが可能である。

方法１９００の別の態様において、位相多項式低減オペレーションを行う前に、アダマールゲートの低減オペレーションを実行することができる。

方法１９００の別の態様において、位相多項式低減オペレーションを行う前に、単一キュービットゲートの相殺オペレーションを実行することができる。

方法１９００の別の態様において、位相多項式低減オペレーションを行う前に、２キュービットゲートの相殺オペレーションを実行することができる。

方法１９００の別の態様において、位相多項式低減オペレーションを行うステップは、書換えルールのセットを実行するステップを含んでよい。書換えルールのセットは、ゲートカウント保持の書換えルールまたはゲートカウント低減の書換えルールのうちの一方またはその両方を含んでよい。

方法１９００の別の態様において、位相多項式低減オペレーションとともに、１つ以上のゲート相殺オペレーションまたはゲート低減オペレーションを、反復的に行うことができる。

方法１９００の別の態様において、位相多項式低減オペレーションを含む最適化オペレーションの固定シーケンスを反復して行うことができ、この位相多項式低減オペレーションは、この固定シーケンス中の最初の最適化オペレーションではなく、この固定シーケンスの中で一回だけ行われる。

図２０は、本開示の態様による、ＱＩＰシステム２００５の一例を示す図表２０００を示す。ＱＩＰシステム２００５は、また、量子コンピューティングシステム、量子コンピューティングネットワーク、コンピューティングデバイスと呼ばれることもある。一態様において、ＱＩＰシステム２００５は、量子コンピューティングオペレーションまたはアルゴリズムをインプリメントまたは実行するために用いることができ、該システムに対する、インプリメントされている量子ゲートの適合性は、例えば、量子ビットとして用いられるトラップイオンに印加される調整された安定なレーザ出力を有する能力に依存する。ＱＩＰシステム２００５は、図１８中のコンピュータデバイス１８１０、および／または図１８中のコンピュータデバイス１８１０中の最適化アプリケーション１８３０（例えば、本提案のオプティマイザまたは最適化アルゴリズム）を実行することによって生成された最適化結果を受信したまたはインプリメントした量子コンピュータの量子コンピュータインプリメンテーションに対応してよい。

このＱＩＰシステム２００５は、量子コンピュータのトラップ原子イオンバージョンを表し、光コントローラ２０２０によってイオン化された原子種（例えば、トラップイオン）を捕捉するイオントラップ２０７０を有するチャンバ２０５０に、原子種を供給する供給源２０６０を含むことができる。光コントローラ２０２０中の光源２３０は、１つ以上のレーザ源を含むことができ、レーザ源は、原子種のイオン化と、原子種の制御（例えば位相制御）と、光コントローラ２０２０中のイメージングシステム２０４０中で動作しているイメージ処理アルゴリズムによってモニタし追跡が可能な原子イオンの蛍光発光と、のために用いることができる。

イメージングシステム２０４０は、原子イオンがイオントラップ２０７０に供給されている間に（カウントするために）、または原子イオンがイオントラップ２０７０に供給された後に（原子イオンの状態をモニタするために）、それらの原子イオンをモニタするための高解像度イメージャ（例えばＣＣＤカメラ）を含むことができる。一態様において、イメージングシステム２０４０は、光コントローラ２０２０と分離してインプリメンテーションが可能であるが、但し、イメージ処理アルゴリズムを使って原子イオンを検出し、識別し、および標識するための蛍光発光の使用は、光コントローラ２０２０と調整する必要があり得る。

また、ＱＩＰシステム２００５は、前述の最適化技法を使用する量子アルゴリズム（例えば、ＱＦＴ、量子シミュレーション）を行うため、ＱＩＰシステム２００５の他の部分（図示せず）とともに動作することが可能なアルゴリズムコンポーネント２０１０も含んでよい。アルゴリズムコンポーネント２０１０は、量子回路またはその等価物のインプリメンテーションを可能にするために、ＱＩＰシステム２００５の様々なコンポーネントに（例えば、光コントローラ２０２０に）命令を与えることができる。すなわち、アルゴリズムコンポーネント２０１０は、各種のコンピューティングプリミチブを、例えば、イオントラップ２０７０中のトラップイオンをキュービットとして使って、物理的表現にマッピングすることを可能にできる。

ＱＩＰシステム２００５および／またはアルゴリズムコンポーネント２０１０を使用して、化学的または分子的用途に関連する量子シミュレーションをインプリメントすることもできる。例えば、一般に、分子などの量子オブジェクトのハミルトニアンは、以下：

のように表すことができる。

上式で、ｐ、ｑ、ｒ、およびｓは電子軌道の整数および／またはラベルを表し、

は生成演算子であり、ａ_ｑ、ａ_ｒ、およびａ_ｓは、消滅演算子である。ハミルトニアンは、例えば、Ｊｏｒｄａｎ－Ｗｉｇｎｅｒ変換：

を介して、量子コンピューティングにより適した形式に変換できる。

小さいｔ（例えば、１よりも著しく小さい）の場合、変換は、一次式が

であるトロッター分割プロセスによって近似することができる。

ｔが大きい場合、式は、

となる。

こうして、括弧内の

は、十分に大きなｒに対して、

でトロッター化され得る。個々のハミルトニアン

は、Ｊｏｒｄａｎ－Ｗｉｇｎｅｒ変換形式（上記を参照）などの変換形式であり得る。

いくつかのインプリメンテーションでは、分子のハミルトニアンをインプリメントすることは、２キュービット量子ゲートなどの各種のゲートの使用を含み得る。２キュービットＣＮＯＴゲートなどの２キュービット量子ゲートは、インプリメンテーションの物理レベルで、アダマールゲートなどの１キュービット量子ゲートよりもリソースを大量に消費する可能性がある。したがって、２キュービット量子ゲートの数の低減は、プレフォールトトレラントな量子コンピューティングにとって重要である可能性がある。

図２１において、量子化学シミュレーションで頻繁に発生する、効果的なハミルトニアン演算子である

の例を示している。ｅ^－ｉΣθをインプリメントする量子回路（上式で、Σは、加数２１１０が図２１のシーケンス２１００に現れるパウリ行列の積の合計である）、トロッター式を使用すると、１つ以上の２キュービット量子ゲートを有する。この例では、奇数のσ_Ｘおよび／またはσ_Ｙを検討するが、例えば、ハミルトン力学シミュレーションを考慮することから生じる可能性のある、積内のそのような行列の偶数も、以下のように検討することができる。加数２１１０を並べ替えることにより、Σを、図２１のシーケンス２１５０に示される順序で現れる加数２１１０で表現することが可能になり得る。シーケンス２１００に示される加数順序でｅ^－ｉΣθをインプリメントする第一の量子回路と、シーケンス２１５０に示される加数順序でｅ^－ｉΣθをインプリメントする第二の量子回路は、トロッターエラーまで同一の機能を実行することができるが、シーケンス２１５０に示される加数を有する第二の量子回路は、特定の低減サブルーチンを使用することによって、より少ない２量子ビット量子ゲートを利用することができる。その例を以下でより詳細に説明する。

（サブルーチン６；制御されたゲート相殺（２））
いくつかのインプリメンテーションでは、シーケンス２１５０の加数２１１０ｈは、図２２Ａに示される量子回路２２００を導くことができる。

いくつかのインプリメンテーションでは、シーケンス２１５０内の加数２１１０ｇは、図２２Ａの量子回路２２３０を導くことができる。図２２Ａ～Ｂは、本開示の一態様による２キュービット量子ゲートの相殺ルーチンを示している。相殺ルーチンは機能を保持する場合がある。第一の量子回路２２００は、６つのＣＮＯＴゲートおよび単一キュービットｚ回転ゲートを含み得、第二の量子回路２２３０は、６つのＣＮＯＴゲートおよび単一キュービットｚ回転ゲートを含み得る。元の項が特定のキュービットに

を含む場合は常に、６つのＣＮＯＴおよび１つのＲ_Ｚ回路は、キュービットライン上のＨによって共役になり得る。元の項が特定のキュービットに

を含む場合は常に、６つのＣＮＯＴおよび１つのＲ_Ｚ回路は、キュービットライン上で、左からＰ^†Ｈによって、また右からＨＰによって共役になり得る。第一の量子回路２２００は、第二の量子回路２２３０の前に実行され得る。

引き続き、図２２Ａを参照すると、非限定的な例では、２つの隣接するＨゲート２２６０は、Ｈゲート相殺オペレーションを使用して相殺され得る。同様に、２つの隣接するＨゲート２２６２は、Ｈゲート相殺オペレーションを使用して相殺され得る。次に、ＣＮＯＴゲート２２６４は、ＨＰＨ位置転換オペレーションを使用して、第一のＨゲート、Ｐゲート、および第二のＨゲートの組み合わせ２２６６で（一時的に）位置転換され得る。ＨＰＨ位置交換は、ＨＰＨゲートの組み合わせ２２６６とＣＮＯＴゲート２２６４との間で実行の順序を切り替えることができる。ＣＮＯＴゲート２２６４の位置交換とＨＰＨゲートの組み合わせ２２６６後、ＣＮＯＴゲート２２６４の制御キュービットは、ＣＮＯＴゲート２２６８の制御キュービットに隣接し得る。ＣＮＯＴゲート２２６４の標的キュービットは、ＣＮＯＴゲート２２６８の標的キュービットに隣接し得る。ＣＮＯＴゲート２２６４、２２６８は両方とも、ＣＮＯＴゲート相殺オペレーションを使用して相殺することができる。

図２２Ｂに目を向けると、特定の例では、第一の量子回路２２００および第二の量子回路２２３０における２キュービット量子ゲートは、相殺することができる。ＨＰＨゲートの組み合わせ２２６６は、ＨＰＨ位置交換オペレーションを使用して、ＣＮＯＴゲート２２７０と位置交換することができる。ＣＮＯＴゲート２２７０の位置交換およびＨＰＨゲートの組み合わせ２２６６後、ＣＮＯＴゲート２２７０の制御キュービットは、ＣＮＯＴゲート２２７２の制御キュービットに隣接し得る。ＣＮＯＴゲート２２７０の標的キュービットは、ＣＮＯＴゲート２７７２の標的キュービットに隣接し得る。ＣＮＯＴゲート２２７０、２２７２は両方とも、ＣＮＯＴゲート相殺オペレーションを使用して相殺することができる。

引き続き、図２２Ｂを参照すると、いくつかのインプリメンテーションでは、第一のＨゲート、Ｐ^†ゲート、および第二のＨゲート（すなわち、ＨＰ^†Ｈ）の組み合わせは、ＨＰ^†Ｈ変換（ＨＰ^†Ｈ→ＰＨＰ）を使用して、第一のＰゲート２２７４、Ｈゲート、第二のＰゲート２２７６（すなわち、ＰＨＰ）の組み合わせに変換することができる。第一の位相ゲート２２７４は、ＣＮＯＴゲート２２７８と位置交換することができ、第二のＰゲート２２７６は、Ｐゲート位置交換オペレーションを使用して、別個にＣＮＯＴゲート２２８０と位置交換することができる。２つのＰゲート位置交換オペレーションの後、ＣＮＯＴゲート２２７８、Ｈゲート、およびＣＮＯＴゲート２２８０の組み合わせ２２８２は、ＣＮＯＴゲート置換を使用して、３つのＨゲート、Ｐ^†ゲート、ＣＮＯＴゲート、および２つのＰゲートを含む組み合わせ２２８４によって置き換えることができる。特定の例では、第一のＨゲート、Ｐゲート、および第二のＨゲート（すなわち、ＨＰＨ）は、ＨＰＨ変換（ＨＰＨ→Ｐ^†ＨＰ^†）において、第一のＰ^†ゲート、Ｈゲート、および第二のＰ^†ゲート（Ｐ^†ＨＰ^†）に変換することができる。いくつかの例では、Ｐ^†ゲート位置交換を使用して、Ｐ^†ゲートをＣＮＯＴゲート（図示せず）と位置交換させることができる。

図２２Ａ～Ｂに示されるこの非限定的な例では、６つのＣＮＯＴゲート（つまり、２キュービット量子ゲート）が１つのＣＮＯＴゲート（６：１の低減比）に削減された。例えば、異なる順序が使用されている場合、あるいはハードウェアの制約のために特定のオペレーションが置き換えられない場合は、５：１、４：１、３：１、または２：１などの他の低減比が可能である。制御されたゲート相殺サブルーチンは、２キュービットの量子ゲート（例えば、ＣＮＯＴゲート）の数を減らしながら、量子回路の機能を維持することができる。２キュービット量子ゲートの数を低減することで、量子回路をインプリメントするのに使用されるリソースを節約できる可能性がある。

いくつかのインプリメンテーションでは、ハードウェアにＣＮＯＴゲートなどの制御されたゲートを直接インプリメントすることは、困難または非効率的である可能性がある。図２３に示すように、例えば、１つ以上のＭφｌｍｅｒ－Ｓφｒｅｎｓｅｎ（ＭＳ）ゲートおよび／または１つ以上の単一回転ゲートを使用して、式２３００に基づいてＣＮＯＴゲートをインプリメントすることができ、ここで、ｓ＝±１はゲートのイオンが適用される相互作用パラメータの符号であり、ν＝±１は任意に選択できる。同様に、ＭＳゲートは、式２３５０に基づいて、４つのＨゲート、２つのＣＮＯＴゲート、および単一キュービットｚ回転ゲートによってインプリメントされ得る。

いくつかの例では、２キュービットゲートの数が同じままであるか、増加し得る場合でさえ、第二の等価回路が、第一の元の回路と比較して、ＣＮＯＴゲート（またはπ／２などの大きな角度のＭＳゲート）の数を減らしながら、小さな角度のＭＳゲートをより頻繁に適用できる場合は、第二の回路をインプリメントすることが望ましい場合がある。例えば、第二の回路は、１以上の熱雑音を含む現実的な雑音源、不正確な制御などが存在する中で、小さな角度のＭＳゲートが大きな角度のＭＳゲートよりもこのようなエラーの原因となるエラーが少ない場合、より優れた性能を示すことがある。

上記のオペレーションをインプリメントする他の例も、また、２キュービット量子ゲートの数を減らすことができる。上述したサブルーチン６およびオペレーションは、例えば、コンピュータデバイス１８１０および／またはＱＩＰシステム２００５によってインプリメントされ得る。

図２４は、本開示の態様による、量子回路内の２キュービット量子ゲートを低減するためのプロセスまたは方法２４００を示すフローチャートである。方法２４００の態様は、図１８に示されるコンピュータデバイス１８１０内のハードウェアおよび／またはソフトウェアによって、および／または図２０のＱＩＰシステム２００５によって、実行され得る。

ブロック２４０５において、方法２４００は、量子回路を形成する第一の複数の２キュービット量子ゲートに関連する情報を含むネットリストを受信するステップを含む。この情報は、制御されたゲート、Ｒ_Ｚゲート、アダマールゲート、位相ゲート、または逆位相ゲートのうちの１つ以上に関連する情報を含み得る。例えば、コンピュータデバイス１８１０は、量子回路構成に関連する情報を含むネットリストを受信することができる。

ブロック２４１０において、方法２４００は、第一の複数の２キュービット量子ゲートと機能的に同等である第二の複数の２キュービット量子ゲートを生成するために、第一の複数の２キュービット量子ゲートに関連する情報に対して制御されたゲート相殺オペレーションを実行するステップを含み、第一の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第一の数は、第二の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第二の数よりも大きい。制御されたゲート相殺オペレーションの例は、図２２Ａ～Ｂに関連するオペレーションを含み得る。最適化アプリケーション１８３０は、２キュービット量子ゲートの数を減らすために、図２２Ａ～Ｂに関連するオペレーションを実行することができる。

ブロック２４１５において、方法２４００は、第二の複数の２キュービット量子ゲートに関する情報を含む新規ネットリストを生成するステップを含む。

ブロック２４２０において、方法２４００は、第二の複数の２キュービット量子ゲートに基づいて、量子回路の機能をインプリメントするための新規ネットリストを提供するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＨゲート相殺オペレーションを実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＨＰＨ位置交換オペレーションを実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＣＮＯＴゲート相殺オペレーションを実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＨＰ^†Ｈ変換を実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＰゲート位置交換オペレーションを実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＨＰＨ変換を実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＰ^†ゲート位置交換オペレーションを実行するステップを含む。

代替のインプリメンテーションでは、方法２４００は、１つ以上のＣＮＯＴゲート置換を実行するステップを含む。

本開示を、図示のインプリメンテーションによって提示してきたが、当業者は、これらの実施形態に対するバリエーションがあり得、それらのバリエーションは、本開示の範囲内にあることになろうことを容易に認識していよう。したがって、当業者は、添付の特許請求の範囲から逸脱することなく、多くの修正をなすことができよう。

Claims (20)

1. 量子回路の２キュービットゲートを低減するための方法であって、
   前記量子回路を形成する第一の複数の２キュービット量子ゲートに関連する情報を含むネットリストを受信するステップと、
   前記第一の複数の２キュービット量子ゲートと機能的に同等である第二の複数の２キュービット量子ゲートを生成するために、前記第一の複数の２キュービット量子ゲートに関連する前記情報に対して制御されたゲート相殺オペレーションを実行するステップであって、前記第一の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第一の数は、前記第二の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第二の数よりも大きい、ステップと、
   前記第二の複数の２キュービット量子ゲートに関する情報を含む新規ネットリストを生成するステップと、
   前記第二の複数の２キュービット量子ゲートに基づいて、前記量子回路の機能をインプリメントするための前記新規ネットリストを提供するステップと、
   含む、方法。
2. 前記情報が、前記制御されたゲート、Ｒ_Ｚゲート、Ｈゲート、Ｐゲート、またはＰ^†ゲートのうちの１つ以上に関連する情報を含む、請求項１に記載の方法。
3. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＨゲート相殺オペレーションを実行するステップをさらに含む、請求項１に記載の方法。
4. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＨＰＨ位置交換オペレーションを実行することをさらに含む、請求項１に記載の方法。
5. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＣＮＴゲート相殺オペレーションを実行するステップをさらに含む、請求項４に記載の方法。
6. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＨＰ^†Ｈ変換を実行するステップをさらに含む、請求項１に記載の方法。
7. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＰゲート位置交換オペレーションを実行するステップをさらに含む、請求項６に記載の方法。
8. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＨＰＨ変換を実行するステップをさらに含む、請求項１に記載の方法。
9. 前記制御されたゲート相殺オペレーションを実行するステップは、１つ以上のＰ^†ゲート位置交換オペレーションを実行するステップをさらに含む、請求項８に記載の方法。
10. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＣＮＯＴゲート置換を実行するステップをさらに含む、請求項１に記載の方法。
11. １つ以上のプロセッサによって実行されると、前記１つ以上のプロセッサに
    量子回路を形成する第一の複数の２キュービット量子ゲートに関連する情報を含むネットリストを受信することと、
    前記第一の複数の２キュービット量子ゲートと機能的に同等である第二の複数の２キュービット量子ゲートを生成するために、前記第一の複数の２キュービット量子ゲートに関連する前記情報に対して制御されたゲート相殺オペレーションを実行することであって、前記第一の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第一の数は、前記第二の複数の２キュービット量子ゲートにおける２キュービット量子ゲートの第二の数よりも大きい、ことと、
    前記第二の複数の２キュービット量子ゲートに関する情報を含む新規ネットリストを生成することと、
    前記第二の複数の２キュービット量子ゲートに基づいて、前記量子回路の機能をインプリメントするための前記新規ネットリストを提供することと、
    を誘発する命令を有する、コンピュータ可読媒体。
12. 前記情報が、前記制御されたゲート、Ｒ_Ｚゲート、Ｈゲート、Ｐゲート、またはＰ†ゲートのうちの１つ以上に関連する情報を含む、請求項１１に記載のコンピュータ可読媒体。
13. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＨゲート相殺オペレーションを実行することをさらに含む、請求項１１に記載のコンピュータ可読媒体。
14. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＨＰＨ位置交換オペレーションを実行することをさらに含む、請求項１１に記載のコンピュータ可読媒体。
15. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＣＮＯＴゲート相殺オペレーションを実行することをさらに含む、請求項１４に記載のコンピュータ可読媒体。
16. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＨＰ†Ｈ変換を実行することをさらに含む、請求項１１に記載のコンピュータ可読媒体。
17. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＰゲート位置交換オペレーションを実行することをさらに含む、請求項１６に記載のコンピュータ可読媒体。
18. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＨＰＨ変換を実行することをさらに含む、請求項１１に記載のコンピュータ可読媒体。
19. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上の逆Ｐゲート位置交換オペレーションを実行することをさらに含む、請求項１８に記載のコンピュータ可読媒体。
20. 前記制御されたゲート相殺オペレーションを実行することは、１つ以上のＣＮＯＴゲート置換を実行することをさらに含む、請求項１１に記載のコンピュータ可読媒体。
    

Images