特許 7324873 トラップイオン型量子コンピュータのための振幅、周波数及び位相変調同時もつれゲート

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-08-02

(45)【発行日】2023-08-10

(54)【発明の名称】トラップイオン型量子コンピュータのための振幅、周波数及び位相変調同時もつれゲート

(51)【国際特許分類】

   G06N 10/40 20220101AFI20230803BHJP

【ＦＩ】

G06N10/40

【請求項の数】 20

(21)【出願番号】P 2021568970

(86)(22)【出願日】2020-05-21

(65)【公表番号】

(43)【公表日】2022-07-25

(86)【国際出願番号】 US2020034010

(87)【国際公開番号】W WO2020237056

(87)【国際公開日】2020-11-26

【審査請求日】2022-01-13

(31)【優先権主張番号】62/851,280

(32)【優先日】2019-05-22

(33)【優先権主張国・地域又は機関】US

(31)【優先権主張番号】62/851,394

(32)【優先日】2019-05-22

(33)【優先権主張国・地域又は機関】US

(31)【優先権主張番号】16/854,043

(32)【優先日】2020-04-21

(33)【優先権主張国・地域又は機関】US

(73)【特許権者】

【識別番号】520132894

【氏名又は名称】イオンキュー インコーポレイテッド

(74)【代理人】

【識別番号】100120891

【弁理士】

【氏名又は名称】林 一好

(74)【代理人】

【識別番号】100165157

【弁理士】

【氏名又は名称】芝 哲央

(74)【代理人】

【識別番号】100205659

【弁理士】

【氏名又は名称】齋藤 拓也

(74)【代理人】

【識別番号】100126000

【弁理士】

【氏名又は名称】岩池 満

(74)【代理人】

【識別番号】100185269

【弁理士】

【氏名又は名称】小菅 一弘

(72)【発明者】

【氏名】ナム ユンソン

(72)【発明者】

【氏名】ブルメル レインホールド

(72)【発明者】

【氏名】グルチェシャク ニコデム

【審査官】渡辺 順哉

(56)【参考文献】

【文献】特表２０２１－５０８３８２（ＪＰ，Ａ）

【文献】特表２０２１－５２７２５３（ＪＰ，Ａ）

【文献】CHOI, T. ほか，Optimal quantum control of multi-mode couplings between trapped ion qubits for scalable entanglement，arXiv[online]，2014年01月08日，pp.1-4，[retrieved on 2023.01.30], Retrieved from the Internet: <URL: https://arxiv.org/pdf/1401.1575.pdf>

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｎ １０／００

(57)【特許請求の範囲】

【請求項1】

量子コンピュータを使用して計算を実行する方法であって、
量子コンピュータにおいてトラップイオンの鎖内の複数の関与イオンの中でイオンの複数ペアを同時にもつれさせるために使用される複数のレーザパルスのそれぞれのパルス関数のフーリエ係数を古典的コンピュータによって計算するステップであって、前記トラップイオンのそれぞれはキュービットを定義する２つの周波数分離状態を有する、ステップと、
計算された前記フーリエ係数を有する前記複数のレーザパルスのそれぞれの前記パルス関数を前記古典的コンピュータによって復調して、振幅及び離調周波数値を計算する、ステップと、
それぞれが計算された前記振幅及び計算された前記離調周波数値を有する前記複数のレーザパルスを前記複数の関与イオンにシステムコントローラによって適用するステップであって、前記複数の関与イオンのそれぞれは、前記複数のレーザパルスのうちの少なくとも１つを受け取る、ステップと、
を含む、方法。

【請求項2】

前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、もつれるイオンの前記複数ペアのそれぞれの間に、ゼロ以外のもつれ相互作用が存在し、もつれることのないイオンのペアのそれぞれは分離されるという条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を前記古典的コンピュータによって選択するステップ
をさらに含む、請求項１に記載の方法。

【請求項3】

前記ゼロ以外のもつれ相互作用がゼロ～π／２の間の値を有する、請求項２に記載の方法。

【請求項4】

前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、トラップイオンの前記鎖内の前記トラップイオンがそれらの元の位置及び運動量値に戻る条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を前記古典的コンピュータによって選択するステップ
をさらに含む、請求項１に記載の方法。

【請求項5】

前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、前記複数のレーザパルスによって前記複数の関与イオンに提供される出力が最小である条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を前記古典的コンピュータによって選択するステップ
をさらに含む、請求項４に記載の方法。

【請求項6】

前記フーリエ係数を計算するステップは、前記複数のレーザパルスのゲート持続時間値を複数の時間セグメントに分割することを含む、請求項１に記載の方法。

【請求項7】

デジタルコンピュータのプロセッサによって、前記デジタルコンピュータのメモリに記憶されているソフトウェアプログラムを実行するステップであって、実行された前記ソフトウェアプログラムは、少なくとも１つの計算を実行する必要があり、少なくとも１つの計算を実行することは、
前記デジタルコンピュータの前記プロセッサによって、トラップイオンの前記鎖に実装される量子アルゴリズムを選択するステップと、
選択した前記量子アルゴリズムを一連のユニバーサル論理ゲートにコンパイルするステップと、
前記一連のユニバーサル論理ゲートをパルスに変換して、トラップイオンの前記鎖内の複数の関与イオンに適用するステップと、
トラップイオンの前記鎖内の前記イオンのキュービット状態の集団を測定するステップと、
前記キュービット状態の測定された前記集団に基づいて、前記デジタルコンピュータの前記プロセッサによって、トラップイオンの前記鎖内のイオンの前記キュービット状態に対応する量子情報を処理するステップと、
を含む、ステップと、
少なくとも一つの計算の処理結果に基づいて、選択した前記量子アルゴリズムの解を生成するステップと、
をさらに含む、請求項１に記載の方法。

【請求項8】

コンピュータプログラム命令を含むコンピュータ可読媒体であって、前記コンピュータプログラム命令は、プロセッサによって実行されると、前記プロセッサに、
トラップイオンの鎖内の複数の関与イオンの中でイオンの複数ペアを同時にもつれさせるために使用される複数のレーザパルスのそれぞれのパルス関数のフーリエ係数を計算するステップであって、前記トラップイオンのそれぞれはキュービットを定義する２つの周波数分離状態を有し、計算された前記フーリエ係数を有する前記複数のレーザパルスのそれぞれの前記パルス関数を復調して、振幅及び離調周波数値を計算する、ステップと、
それぞれが計算された前記振幅及び計算された前記離調周波数値を有する前記複数のレーザパルスを前記複数の関与イオンに適用するステップであって、前記複数の関与イオンのそれぞれは、前記複数のレーザパルスのうちの少なくとも１つを受け取る、ステップと、
を行わせる、コンピュータ可読媒体。

【請求項9】

前記コンピュータプログラム命令は、前記プロセッサに、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、もつれるイオンの前記複数ペアのそれぞれの間に、ゼロ以外のもつれ相互作用が存在し、もつれることのないイオンのペアのそれぞれは分離されるという条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに行わせる、請求項８に記載のコンピュータ可読媒体。

【請求項10】

前記ゼロ以外のもつれ相互作用がゼロ～π／２の間の値を有する、請求項９に記載のコンピュータ可読媒体。

【請求項11】

前記コンピュータプログラム命令は、前記プロセッサに、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、トラップイオンの前記鎖内の前記トラップイオンがそれらの元の位置及び運動量値に戻る条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに行わせる、請求項９に記載のコンピュータ可読媒体。

【請求項12】

前記コンピュータプログラム命令は、前記プロセッサに、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、前記複数のレーザパルスによって前記複数の関与イオンに提供される出力が最小である条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに行わせる、請求項１１に記載のコンピュータ可読媒体。

【請求項13】

前記フーリエ係数を計算するステップは、前記複数のレーザパルスのゲート持続時間値を複数の時間セグメントに分割することを含む、請求項９に記載のコンピュータ可読媒体。

【請求項14】

トラップイオンの鎖であって、前記トラップイオンのそれぞれは、キュービット状態と励起状態を定義する２つの超微細状態を有する、トラップイオンの鎖と、
内部に記憶されたいくつかの命令を有するメモリを備えるコントローラと、
を備える、量子コンピューティングシステムであって、前記命令は、前記量子コンピューティングシステムに、
トラップイオンの鎖内の複数の関与イオンの中でイオンの複数ペアを同時にもつれさせるために使用される複数のレーザパルスのそれぞれのパルス関数のフーリエ係数を計算するステップと、
計算された前記フーリエ係数を有する前記複数のレーザパルスのそれぞれの前記パルス関数を復調して、振幅及び離調周波数値を計算する、ステップと、
それぞれが計算された前記振幅及び計算された前記離調周波数値を有する前記複数のレーザパルスを前記複数の関与イオンに適用するステップであって、前記複数の関与イオンのそれぞれは、前記複数のレーザパルスのうちの少なくとも１つを受け取る、ステップと、
を含む操作を実行させる、量子コンピューティングシステム。

【請求項15】

前記トラップイオンのそれぞれは、^２Ｓ_１／２の超微細状態を有する^１７１Ｙｂ^＋である、
請求項１４に記載の量子コンピューティングシステム。

【請求項16】

前記操作は、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、もつれるイオンの前記複数ペアのそれぞれの間に、ゼロ以外のもつれ相互作用が存在し、もつれることのないイオンのペアのそれぞれは分離されるという条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに含む、請求項１４に記載の量子コンピューティングシステム。

【請求項17】

前記ゼロ以外のもつれ相互作用がゼロ～π／２の間の値を有する、請求項１６に記載の量子コンピューティングシステム。

【請求項18】

前記操作は、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、トラップイオンの前記鎖内の前記トラップイオンがそれらの元の位置及び運動量値に戻る条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに含む、請求項１６に記載の量子コンピューティングシステム。

【請求項19】

前記操作は、
前記複数のレーザパルスのそれぞれの前記パルス関数を復調する前に、前記複数のレーザパルスによって前記複数の関与イオンに提供される出力が最小である条件に基づいて、前記複数のレーザパルスのそれぞれの前記パルス関数の前記フーリエ係数を選択するステップ
をさらに含む、請求項１８に記載の量子コンピューティングシステム。

【請求項20】

前記フーリエ係数を計算するステップは、前記複数のレーザパルスのゲート持続時間値を複数の時間セグメントに分割することを含む、請求項１４に記載の量子コンピューティングシステム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、一般に、イオントラップ型量子コンピュータにおいて、もつれゲート操作を実行する方法に関し、より具体的には複数のペアワイズ（ｐａｉｒ－ｗｉｓｅ）もつれゲート操作を同時に実行するパルスを構築する方法に関する。

【背景技術】

【0002】

量子コンピューティングでは、古典的デジタルコンピュータにおける「０」と「１」を表すビットに類似した量子ビット又はキュービットの状態をほぼ完全な精度で決定するために、量子ビット又はキュービットを計算プロセス中に準備し、操作し、測定（読み出し）する必要がある。キュービットの制御が不完全であると、計算プロセスにおいて誤差が蓄積することがあり、信頼性の高い計算を実行できる量子コンピュータのサイズが制限される。

【0003】

大規模な量子コンピュータを構築するために提案されている物理システムの中に、電磁界によってトラップされて真空中に浮遊するイオンの鎖（例えば、電荷を帯びた原子）がある。イオンは、数ＧＨｚ範囲内の周波数によって分離され、キュービットの計算状態（「キュービット状態」と呼ばれる）として使用することができる内部超微細状態を有する。これらの超微細状態は、レーザから提供される放射線を使用して（場合によっては本明細書ではレーザビームとの相互作用と呼ばれることもある）制御し読み取ることができる。イオンは、また、このようなレーザ相互作用を使用して、運動基底状態の近くまで冷却することができる。個々のキュービットは、２つの超微細状態のいずれかに高精度で光学的に励起し（キュービットの準備）、レーザビームにより２つの超微細状態間で操作することができ（単一キュービットのゲート操作）、共鳴レーザビームの適用時に蛍光によってそれらの内部超微細状態が検出される。１ペアのイオンは、イオン間のクーロン力の相互作用により発生する、トラップイオンの鎖の集合運動モードにイオンを結合するレーザパルスを使用して、キュービット状態に依存する力によって制御可能にもつれることができる（２キュービットのゲート操作）。一般に、もつれは、イオン（又は粒子）のペア又はグループが生成、相互作用、又は空間的近接性を共有するときに発生する。そのため、イオンが大きな距離で隔てられていても、各イオンの量子状態を他のイオンの量子状態とは独立して記述できない。量子コンピュータのサイズが大きくなると、１ペアのイオン間の２キュービットゲート演算の実装が複雑になるため、実装に伴う誤差や、レーザ出力などの実装に必要なリソースが増加する。

【0004】

古典的コンピュータでは難解な問題を解決するアルゴリズムを実装できる可能性のある量子コンピュータのサイズを大きくするためには、キュービットを正確に制御して、最小限のリソースで所望の計算プロセスを実行するための手順が必要である。

【発明の概要】

【0005】

量子コンピュータを使用して計算を実行する方法は、第一の方向に整列された複数のトラップイオンのそれぞれに個別に適用されるために使用される複数のレーザパルスを生成するステップであって、トラップイオンのそれぞれはキュービットを定義する２つの周波数分離状態を有する、ステップと、生成された複数のレーザパルスを複数のトラップイオンに適用して、複数のトラップイオンに対して同時にペアワイズもつれゲート操作を実行するステップと、を含む。複数のレーザパルスを生成するステップは、複数のレーザパルスによって引き起こされる複数のトラップイオンにおけるペアワイズもつれ相互作用の値に基づいて、複数のレーザパルスのそれぞれの振幅値及び離調周波数値を調整することを含む。

【0006】

コンピュータプログラム命令を含む非一時的なコンピュータ可読媒体は、命令がプロセッサによって実行されると、プロセッサに、第一の方向に整列された複数のトラップされたイオンのそれぞれに個別に適用されるために使用される複数のレーザパルスを生成するステップであって、トラップイオンのそれぞれはキュービットを定義する２つの周波数分離状態を有する、ステップと、生成された複数のレーザパルスを複数のトラップイオンに適用して、複数のトラップイオンに対して同時にペアワイズもつれゲート操作を実行するステップと、を行わせる。複数のレーザパルスを生成するステップは、複数のレーザパルスによって引き起こされる複数のトラップイオンにおけるペアワイズもつれ相互作用の値に基づいて、複数のレーザパルスのそれぞれの振幅値及び離調周波数値を調整することを含む。

【0007】

量子コンピューティングシステムは、第一の方向に整列された複数のトラップイオンであって、トラップイオンのそれぞれはキュービットを定義する２つの周波数分離状態を有する、複数のトラップイオンと、内部に記憶されたいくつかの命令を有する不揮発性メモリを備えるコントローラと、を含み、命令は、プロセッサによって実行されると、量子コンピューティングシステムに、複数のトラップイオンのそれぞれに個別に適用されるために使用される複数のレーザパルスを生成するステップと、生成された複数のレーザパルスを複数のトラップイオンに適用して、複数のトラップイオンに対して同時にペアワイズもつれゲート操作を実行するステップと、を実行させる。複数のレーザパルスを生成するステップは、複数のレーザパルスによって引き起こされる複数のトラップイオンにおけるペアワイズもつれ相互作用の値に基づいて、複数のレーザパルスのそれぞれの振幅値及び離調周波数値を調整することを含む。

【図面の簡単な説明】

【0008】

本開示の上記特徴を詳細に理解することができるように、上で簡単に要約された本開示のより具体的な記載は、いくつかが添付の図面に示されている実施形態を参照することによって説明することができる。しかしながら、添付の図面は、本開示の典型的な実施形態のみを説明しており、その範囲を限定すると見なされるべきではないことに留意されたい。なぜなら、本開示は、他の同等に有効な実施形態を認めることができるからである。

【0009】

【図1】一実施形態に従うイオントラップ型量子コンピュータの部分図である。

【図2】一実施形態に従って、イオンを鎖に閉じ込めるためのイオントラップの概略図を示す。

【図3】図３Ａ、図３Ｂ、及び図３Ｃは、５つのトラップイオンの鎖のいくつかの概略的な集合横運動モード構造を示す。

【図4】一実施形態に従って、トラップイオンの鎖内の各イオンの概略エネルギー図を示す。

【図5】ブロッホ球の表面上の点として表されるイオンのキュービット状態を示す。

【図6】図６Ａ及び図６Ｂは、一実施形態に従って、各イオンの運動側波帯スペクトル及び運動モードの概略図を示す。

【図7】一実施形態に従って、制御パラメータのセット（ＥＡＳＥプロトコル）を決定するための方法を示すフローチャートを示す。

【0010】

理解を容易にするために、可能な場合には、図に共通する同一の要素を示すために同一の参照番号を使用する。図及び以下の説明では、Ｘ軸、Ｙ軸、及びＺ軸を含む直交座標系を使用する。図面の矢印で表される方向は、便宜上、正の方向であると想定される。いくつかの実施形態で開示された要素は、具体的な明記なく、他の実装で有益に利用されてよいと考えられる。

【発明を実施するための形態】

【0011】

本明細書に記載の実施形態は、一般に、量子計算中にイオンの複数ペアに対してもつれ操作を同時に実行するためのパルスを設計、最適化、及び送達する方法及びシステムに関し、より具体的には、効率的な方法で構築でき、さらに、もつれゲート操作を実行するのに必要なレーザ出力を減らすことのできるパルスに関する。

【0012】

トラップイオンを使用して量子コンピューティングを実行できるシステム全体には、古典的コンピュータ、システムコントローラ、及び量子レジスタが含まれる。古典的コンピュータは、グラフィックス処理ユニット（ＧＰＵ）などのユーザインターフェイスを使用して実行する量子アルゴリズムの選択、選択した量子アルゴリズムの一連のユニバーサル論理ゲートへのコンパイル、一連のユニバーサル論理ゲートを量子レジスタに印加するためのレーザパルスへの変換、中央処理ユニット（ＣＰＵ）を使用してレーザパルスを最適化するパラメータの事前の計算などのサポート及びシステム制御タスクを実行する。量子アルゴリズムを分解して実行するタスクを実行するためのソフトウェアプログラムは、古典的コンピュータ内の不揮発性メモリに記憶されている。量子レジスタには、様々なハードウェアと結合されたトラップイオンが含まれ、これらのハードウェアには、トラップイオンの内部超微細状態（キュービット状態）を操作するレーザや、トラップイオンの内部超微細状態（キュービット状態）を読み取る音響光学変調器が含まれる。システムコントローラは、古典的コンピュータから、量子レジスタで選択されたアルゴリズムの実行の開始時に出力最適パルスの事前計算されたパラメータを受け取り、量子レジスタで選択されたアルゴリズムを実行するために使用されるすべての側面の制御に関連する様々なハードウェアを制御し、量子レジスタの読み取り値を返し、こうして、アルゴリズムの実行の最後に、量子コンピューティングの結果を古典的コンピュータに出力する。

【0013】

本明細書に記載の方法及びシステムは、論理ゲートを量子レジスタに適用されるレーザパルスに変換するためのプロセスと、量子レジスタに適用され、量子コンピュータの性能を向上させるために使用されるレーザパルスを最適化するパラメータを事前計算するためのプロセスを含む。

【0014】

任意の量子アルゴリズムを分解することができるユニバーサル論理ゲートのいくつかの既知のセットのうち、一般的に｛Ｒ，ＸＸ｝と表記されるユニバーサル論理ゲートのセットは、本明細書に記載されているトラップイオンの量子コンピューティングシステムに固有のものである。ここで、Ｒゲートは、トラップイオンの個々の量子状態の操作に対応し、ＸＸゲート（「もつれゲート」、または「ペアワイズもつれゲート」とも呼ぶ）は、２つのトラップイオンのもつれ操作に対応する。当業者にとって明らかであるように、Ｒゲートは、ほぼ完全な忠実度で実装できるが、ＸＸゲートの形成は、複雑なので、ＸＸゲートの忠実度を向上させ、量子コンピュータ内の計算の誤差を回避又は削減するためには、いくつかの要因を挙げれば、トラップイオンの所定のタイプと、トラップイオンの鎖内のイオンの数と、トラップイオンがトラップされるハードウェア及び環境との最適化が必要である。以下では、向上した忠実度を有するＸＸゲートの形成に基づいて計算を実行するために使用されるパルスを生成し，最適化する方法を説明する。

【0015】

量子コンピュータのサイズが大きくなるにつれて、量子計算を実行するために使用されるもつれゲート操作がますます複雑さになり、これらのもつれゲート操作を実行するために使用されるパルスもますます複雑になる。そのような複雑なパルスを実装するために必要なレーザ出力は、その後増加し、したがって、利用可能なレーザ出力は、実装できる量子コンピュータのサイズを制限する可能性がある。本開示で説明される方法及びシステムは、パルスの構築を単純化し、パルスを実装するのに必要なレーザ出力をさらに低減して、量子コンピュータをより大きなサイズにスケールアップして、より複雑な計算操作を実行できるようにする。これは、所与の出力バジェットに対して、もつれゲートの実行が高速であることを意味する。提供されるレーザ出力に比例する誤差は、必要なレーザ出力が小さいほど減少する。

【0016】

（一般的なハードウェア構成）
図１は、一実施形態に係るイオントラップ型量子コンピュータシステム、又はシステム１００の部分図である。システム１００は、古典的コンピュータ１０１と、システムコントローラ１１８と、Ｚ軸に沿って延びる、トラップイオン（例えば、５つを示す）の鎖１０２である量子レジスタとを含む。古典的コンピュータ１０１は、中央処理ユニット（ＣＰＵ）、メモリ、及びサポート回路（及びＩ／Ｏ）を含む。メモリは、ＣＰＵに接続されており、読み取り専用メモリ（ＲＯＭ）、ランダムアクセスメモリ（ＲＡＭ）、フロッピーディスク、ハードディスク、及びその他の形式のデジタルストレージなどで、ローカル及びリモートで、すぐに利用できるメモリの１つ及び複数であり得る。ソフトウェア命令、アルゴリズム、及びデータは、ＣＰＵに命令するために符号化され、メモリ内に記憶される。サポート回路（図示せず）も、従来の方法でプロセッサをサポートするためにＣＰＵに接続されている。サポート回路は、従来のキャッシュ、電源、クロック回路、入力／出力回路、サブシステムなどを含み得る。

【0017】

例えば、開口数（ＮＡ）が０．３７の対物レンズなどのイメージング対物レンズ１０４は、イオンからＹ軸に沿って蛍光を収集し、個々のイオンを測定するために、各イオンをマルチチャネル光電子増倍管（ＰＭＴ）１０６にマッピングする。Ｘ軸に沿って提供される、レーザ１０８からの非共伝搬ラマンレーザビームは、イオンに対して操作を実行する。回折ビームスプリッタ１１０は、マルチチャネル音響光学変調器（ＡＯＭ）１１４を使用して個別に切り替えられる静的ラマンビーム１１２のアレイを作成し、かつ個々のイオンに選択的に作用するように構成される。グローバルラマンレーザビーム１１６は、すべてのイオンを一度に照射する。システムコントローラ（「ＲＦコントローラ」とも呼ばれる）１１８は、ＡＯＭ１１４を制御する。システムコントローラ１１８は、中央処理ユニット（ＣＰＵ）１２０、読み取り専用メモリ（ＲＯＭ）１２２、ランダムアクセスメモリ（ＲＡＭ）１２４、記憶ユニット１２６などを含む。ＣＰＵ１２０は、ＲＦコントローラ１１８のプロセッサである。ＲＯＭ１２２は、様々なプログラムを記憶し、ＲＡＭ１２４は、様々なプログラム及びデータの作業メモリである。記憶ユニット１２６は、ハードディスクドライブ（ＨＤＤ）又はフラッシュメモリなどの不揮発性メモリを含み、電源が切られても様々なプログラムを記憶する。ＣＰＵ１２０、ＲＯＭ１２２、ＲＡＭ１２４、及び記憶ユニット１２６は、バス１２８を介して相互接続されている。ＲＦコントローラ１１８は、ＲＯＭ１２２又は記憶ユニット１２６に記憶され、ＲＡＭ１２４を作業領域として使用する制御プログラムを実行する。制御プログラムは、データの受信、分析、及び本明細書で説明されたイオントラップ型量子コンピュータシステム１００を作成するために使用される方法及びハードウェアの全ての態様の制御に関連する様々な機能を実行するためにプロセッサによって実行することができるプログラムコード（例えば、命令）を含む１つ以上のソフトウェアアプリケーションを含む。

【0018】

図２は、一実施形態に係る、鎖１０２内にイオンを閉じ込めるイオントラップ２００（ポールトラップとも呼ばれる）の概略図を示す。閉じ込め電位は、静的（ＤＣ）電圧と無線周波数（ＲＦ）電圧の両方によって印加される。静的（ＤＣ）電圧ＶＳがエンドキャップ電極２１０及び２１２に印加されて、Ｚ軸（「軸方向」又は「長手方向」とも呼ばれる）に沿ってイオンを閉じ込める。鎖１０２内のイオンは、イオン間のクーロン相互作用のために、軸方向にほぼ均等に分布している。いくつかの実施形態では、イオントラップ２００は、Ｚ軸に沿って延びる４つの双曲線形状の電極２０２、２０４、２０６、及び２０８を含む。

【0019】

操作中、（振幅Ｖ_ＲＦ／２を有する）正弦波電圧Ｖ_１は、対向する一対の電極２０２、２０４に印加され、正弦波電圧Ｖ_１から１８０°の位相シフト（及び振幅Ｖ_ＲＦ／２を有する正弦波電圧Ｖ_２は、駆動周波数ω_ＲＦで対向する他対の電極２０６、２０８に印加されて、四重極電位を生成する。いくつかの実施形態では、正弦波電圧は、対向する一対の電極２０２、２０４のみに印加され、対向する他対の２０６、２０８は、接地される。四重極電位は、トラップされた各イオンに対してＺ軸に垂直なＸ－Ｙ平面（「半径方向」又は「横方向」とも呼ばれる）に有効な閉じ込め力を生成し、その閉じ込め力は、ＲＦ電界が消失する鞍点（すなわち、軸方向（Ｚ方向）の位置）からの距離に比例する。各イオンの半径方向（すなわち、Ｘ－Ｙ平面の方向）の運動は、半径方向の鞍点に向かう復元力を伴う調和振動（経年運動と呼ばれる）として近似され、それぞれ以下でより詳細に説明されるようなばね定数ｋｘとｋｙによってモデル化できる。いくつかの実施形態では、半径方向のばね定数は、四重極電位が半径方向に対称である場合に等しいものとしてモデル化される。しかしながら、望ましくない場合には、半径方向のイオンの運動は、物理的なトラップ構成のある程度の非対称性、電極の表面の不均一性による小さなＤＣパッチ電位などのために歪む場合があり、これら及び他の外部の歪みの原因により、イオンは、鞍点から中心を外れる場合がある。

【0020】

（トラップイオン構成及びキュービット情報）
図３Ａ、図３Ｂ、及び図３Ｃは、例えば、５つのトラップイオンの鎖１０２のいくつかの概略的な集合横運動モード構造（単に「運動モード構造」とも呼ばれる）を示す。本明細書では、エンドキャップ電極２１０及び２１２に印加された静的電圧Ｖｓによる閉じ込め電位は、半径方向の閉じ込め電位と比較して弱い。トラップイオンの鎖１０２の横方向の集合運動モードは、イオントラップ２００によって生成された閉じ込め電位とトラップイオン間のクーロン相互作用との組み合わせによって決定される。トラップイオンは、集合横方向運動（「集合横運動モード」、「集合運動モード」、又は単に「運動モード」と呼ばれる）を起こし、各モードには、それに関連する異なるエネルギー（及び同等に、周波数）がある。以下では、エネルギーがｐ番目に低い運動モードを│ｎ_ｐｈ＞_ｐと呼び、ここで、ｎ_ｐｈは、運動モードの運動量子の数（エネルギー励起の単位で、フォノンと呼ばれる）を表し、所定の横方向の運動モードの数Ｐは、鎖１０２内のトラップイオンの数Ｎに等しい。図３Ａ～図３Ｃは、鎖１０２内に配置された５つのトラップイオンによって経験され得る異なるタイプの集合横運動モードの例を概略的に説明する。図３Ａは、最も高いエネルギーを有する一般的な運動モード│ｎ_ｐｈ＞_Ｐの概略図であり、ここで、Ｐは、運動モードの数である。一般的な運動モード│ｎ＞_Ｐでは、すべてのイオンは、横方向に同位相で振動する。図３Ｂは、２番目に高いエネルギーを有する傾斜運動モード│ｎ_ｐｈ＞_Ｐ－１の概略図である。傾斜運動モードでは、両端のイオンは、横方向に位相がずれて（すなわち、反対方向に）移動する。図３Ｃは、傾斜運動モード│ｎ_ｐｈ＞_Ｐ－１よりもエネルギーが低く、イオンがより複雑なモードパターンで移動する高次運動モード│ｎ_ｐｈ＞_Ｐ－３の概略図である。

【0021】

なお、上記特定の構成は、本開示のイオンを閉じ込めるトラップのいくつかの可能な例のうちの１つに過ぎず、本開示によるトラップの可能な構成、仕様などを限定するものではない。例えば、電極の形状は、上記双曲線電極に限定されない。他の例では、調和振動として半径方向にイオンの運動を引き起こす実効電界を生成するトラップは、複数の電極層が積層され、対角線上にある２つの電極にＲＦ電圧が印加される多層トラップでも、あるいは全ての電極がチップ上の単一平面に配置されている表面トラップであってよい。さらに、トラップは、複数のセグメントに分割することができ、その隣接するペアが１つ以上のイオンを往復させてリンクすることも、あるいは光子相互接続によって結合することもできる。トラップは、また、微細加工されたイオントラップチップ上に互いに近接して配置された個々のトラップ領域のアレイであってよい。いくつかの実施形態では、四重極電位は、上記ＲＦ成分に加えて、空間的に変化するＤＣ成分を有する。

【0022】

図４は、一実施形態に係る、トラップイオンの鎖１０２内の各イオンの概略エネルギー図４００を示す。一例では、各イオンは、ω_０１／２π＝１２．６４２８２１ＧＨｚの周波数差（「キャリア周波数」と呼ばれる）に対応するエネルギー分割を有する^２Ｓ_１／２超微細状態（すなわち、２つの電子状態）を有する正のイッテルビウムイオン^１７１Ｙｂ^＋であってよい。キュービットは、│０＞と│１＞で表される２つの超微細状態で形成され、超微細基底状態（すなわち、^２Ｓ_１／２超微細状態のうちの低エネルギー状態）が│０＞を表すために選択される。以下、「超微細状態」、「内部超微細状態」及び「キュービット」という用語は、│０＞と│１＞を表すために交換可能に使用されることがある。各イオンは、ドップラー冷却又は分解サイドバンド冷却などの既知のレーザ冷却方法で、フォノン励起なし（すなわち、ｎ_ｐｈ＝０）で任意の運動モードｐの運動基底状態│０＞ｐの近くまで冷却し（すなわち、イオンの運動エネルギーが低下することができる）、次にキュービット状態が光ポンピングによって超微細基底状態│０＞で準備することができる。ここで、│０＞は、トラップイオンの個々のキュービット状態を表し、下付き文字ｐが付いた│０＞_ｐは、トラップイオンの鎖１０２の運動モードｐの運動基底状態を表す。

【0023】

各トラップイオンの個々のキュービット状態は、例えば、励起された^２Ｐ_１／２レベル（｜ｅ＞で表される）を介して３５５ナノメートル（ｎｍ）のモードロックレーザ（ｍｏｄｅ－ｌｏｃｋｅｄ ｌａｓｅｒ）によって操作することができる。図４に示すように、レーザからのレーザビームは、ラマン構成で一対の非共伝搬レーザビーム（周波数ω_１を有する第一のレーザビームと周波数ω_２を有する第二のレーザビーム）に分割され、図４で説明するように、｜０＞と｜ｅ＞の間の遷移周波数ω_０ｅに関して、一光子遷移離調周波数Δ＝ω_１－ω_０ｅによって離調されてよい。二光子遷移離調周波数δは、トラップイオンに第一及び第二のレーザビームによって提供されるエネルギー量の調整を含み、それらを組み合わせて使用すると、トラップイオンが超微細状態｜０＞と｜１＞との間で移動する。一光子遷移離調周波数Δが二光子遷移離調周波数（単に「離調周波数」とも呼ばれる）δ＝ω_１－ω_２－_ω０１（以下、±μで表され、μは正の値である）よりもはるかに大きい場合、それぞれ状態｜０＞と｜ｅ＞の間、状態｜１＞と｜ｅ＞の間でラビフロップが発生する単一光子ラビ周波数Ω_０ｅ（ｔ）とΩ_１ｅ（ｔ）（時間に依存し、第一と第二のレーザビームの振幅と位相によって決定される）、励起状態｜ｅ＞からの自然放出率、２つの超微細状態│０＞と│１＞の間のラビフロップ（「キャリア遷移」と呼ばれる）は、二光子ラビ周波数Ω（ｔ）で誘導される。二光子ラビ周波数Ω（ｔ）は、Ω_０ｅΩ_１ｅ／２Δに比例する強度（すなわち、振幅の絶対値）を有し、ここで、Ω_０ｅとΩ_１ｅは、それぞれ第一と第二のレーザビームによる単一光子ラビ周波数である。以下、キュービットの内部超微細状態（キュービット状態）を操作するためのラマン構成におけるこの非共伝搬レーザビームのセットは、「複合パルス」又は単に「パルス」と呼ばれてよく、結果として生じる二光子ラビ周波数Ω（ｔ）の時間依存パターンは、パルスの「振幅」又は単に「パルス」と呼ばれてよく、それらは、以下で図示され、さらに説明される。離調周波数δ＝ω_１－ω_２－ω_０１は、複合パルスの離調周波数又はパルスの離調周波数と呼ばれてよい。第一及び第二のレーザビームの振幅によって決定される二光子ラビ周波数Ω（ｔ）の振幅は、複合パルスの「振幅」と呼ばれてよい。

【0024】

なお、本明細書に提供される説明で使用される特定の原子種は、イオン化されたときに安定し、かつ明確に定義された２レベルエネルギー構造と、光学的にアクセス可能な励起状態とを有する原子種の一例にすぎないため、本開示のイオントラップ型量子コンピュータの可能な構成、仕様などを限定することを意図するものではない。例えば、他のイオン種は、アルカリ土類金属イオン（Ｂｅ^＋、Ｃａ^＋、Ｓｒ^＋、Ｍｇ^＋、及びＢａ^＋）又は遷移金属イオン（Ｚｎ^＋、Ｈｇ^＋、Ｃｄ^＋）を含む。

【0025】

図５は、方位角φ及び極性角θを有するブロッホ球５００の表面上の点として表されるイオンのキュービット状態を視覚化するのを助けるために提供される。上記のように、複合パルスを適用すると、キュービット状態│０＞（ブロッホ球の北極として表される）と│１＞（ブロッホ球の南極として表される）との間でラビフロップが発生する。複合パルスの持続時間と振幅を調整すると、キュービット状態を│０＞から│１＞に（すなわち、ブロッホ球の北極から南極へ）反転させるか、あるいはキュービット状態│１＞から│０＞に（すなわち、ブロッホ球の南極から北極に）反転させる。複合パルスのこの適用は、「πパルス」と呼ばれる。さらに、複合パルスの持続時間と振幅を調整することにより、キュービット状態│０＞を、２つのキュービット状態│０＞と│１＞が加算され、同位相で均等に重み付けされた重ね合わせ状態│０＞＋│１＞（重ね合わせ状態の正規化係数は、一般性を失うことなく、以下省略される）に変換することができ、そして、キュービット状態│１＞を、２つのキュービット状態│０＞と│１＞が加算され、均等に重み付けされているが、位相がずれる重ね合わせ状態│０＞－│１＞に変換することができる。複合パルスのこの適用は、「π／２パルス」と呼ばれる。より一般的には、加算されて均等に重み付けされた２つのキュービット状態│０＞と│１＞の重ね合わせは、ブロッホ球の赤道上にある点によって表される。例えば、重ね合わせ状態│０＞±│１＞は、方位角φがそれぞれゼロとπである赤道上の点に対応する。方位角φの赤道上の点に対応する重ね合わせ状態は、│０＞＋ｅ^ｉφ│１＞（例えば、φ＝±π／２の場合は│０＞±ｉ│１＞である）として表される。赤道上の２点間の変換（すなわち、ブロッホ球のＺ軸の周りの回転）は、複合パルスの位相をシフトすることで実装できる。

【0026】

イオントラップ型量子コンピュータでは、運動モードは、２つのキュービット間のもつれを仲介するデータバスとして機能することができ、このもつれは、ＸＸゲート操作を実行するために使用される。つまり、２つのキュービットのそれぞれが運動モードともつれて、そして、以下に説明するように、もつれは、運動側波帯励起を使用することによって、２つのキュービット間のもつれに転送される。図６Ａ及び図６Ｂは、一実施形態に係る、周波数ω_ｐを有する運動モード│ｎ_ｐｈ＞ｐでの鎖１０２内のイオンの運動側波帯スペクトルの図を概略的に示す。図６Ｂに示すように、複合パルスの離調周波数がゼロの場合（すなわち、第一と第二のレーザビーム間の周波数差がキャリア周波数δ＝ω_１－ω_２－ω_０１＝０に調整される場合）、キュービット状態│０＞と│１＞の間で単純なラビフロップ（キャリア遷移）が発生する。複合パルスの離調周波数が正の場合（すなわち、第一と第二のレーザビーム間の周波数差が、キャリア周波数よりも高く調整されている場合、δ＝ω_１－ω_２－ω_０１＝μ＞０、青側波帯と呼ばれる）、組み合わされたキュービット運動状態│０＞│ｎ_ｐｈ＞_ｐと│１＞│ｎ_ｐｈ＋１＞_ｐの間でラビフロップが発生する（すなわち、キュービット状態│０＞が│１＞に反転する場合、│ｎ_ｐｈ＞_ｐで表されるｎ_ｐｈ－フォノン励起を伴うｐ番目の運動モードから│ｎ_ｐｈ＋１＞_ｐで表されるｎ_ｐｈ＋１－フォノン励起を伴うｐ番目の運動モードへの遷移が発生する）。複合パルスの離調周波数が負の場合（すなわち、第一と第二のレーザビーム間の周波数差が、運動モード│ｎ_ｐｈ＞_ｐの周波数ω_ｐによってキャリア周波数よりも低く調整されている場合、δ＝ω_１－ω_２－ω_０１＝－μ＜０、赤側波帯と呼ばれる）、組み合わされたキュービット運動状態│０＞│ｎ_ｐｈ＞_ｐと│１＞│ｎ_ｐｈ－１＞_ｐの間のラビフロップが発生する（すなわち、キュービット状態│０＞から│１＞に反転する場合、運動モード│ｎ_ｐｈ＞_ｐから、フォノン励起が１つ少ない運動モード│ｎ_ｐｈ－１＞_ｐへの遷移が発生する）。キュービットに適用された青側波帯のπ／２パルスは、組み合わされたキュービット運動状態│０＞│ｎ_ｐｈ＞_ｐを、│０＞│ｎ_ｐｈ＞_ｐと│１＞│ｎ_ｐｈ＋１＞ｐの重ね合わせに変換する。キュービットに適用された赤側波帯のπ／２パルスは、組み合わされたキュービット運動状態│０＞│ｎ_ｐｈ＞_ｐを、│０＞│ｎ_ｐｈ＞_ｐと│１＞│ｎ_ｐｈ－１＞_ｐの重ね合わせに変換する。二光子ラビ周波数Ω（ｔ）が離調周波数δ＝ω_１－ω_２－ω_０１＝±μと比較して小さい場合、青側波帯遷移又は赤側波帯遷移を選択的に駆動することができる。したがって、キュービットは、π／２パルスなどの適切なタイプのパルスを適用することにより、所望の運動モードでもつれることができ、その後、別のキュービットともつれることができ、２つのキュービットのペアワイズもつれ相互作用をもたらす。イオントラップ型量子コンピュータでＸＸゲート操作を実行するには、キュービットのペアワイズもつれ相互作用が必要である。

【0027】

上記のように、組み合わされたキュービット運動状態の変換を制御及び／又は指示することにより、２つのキュービット（ｉ番目及びｊ番目のキュービット）に対してＸＸゲート操作を実行することができる。一般に、（最大もつれを有する）ＸＸゲート操作は、２キュービット状態｜０＞_ｉ｜０＞_ｊ、｜０＞_ｉ｜１＞_ｊ、｜１＞_ｉ｜０＞_ｊ及び｜１＞_ｉ｜１＞_ｊをそれぞれ次のように変換する。

【数1】

例えば、２つのキュービット（ｉ番目とｊ番目のキュービット）が両方とも最初に超微細基底状態｜０＞（｜０＞_ｉ｜０＞_ｊで表される）にあり、その後、青側波帯のπ／２パルスがｉ番目のキュービットに適用される場合、ｉ番目のキュービットと運動モード｜０＞_ｉ｜ｎ_ｐｈ＞_ｐの組み合わせ状態は、｜０＞_ｉ｜ｎ_ｐｈ＞_ｐと｜１＞_ｉ｜ｎ_ｐｈ＋１＞_ｐの重ね合わせに変換されるため、２つのキュービットと運動モードの組み合わせ状態は、｜０＞_ｉ｜０＞_ｊ｜ｎ_ｐｈ＞_ｐと｜１＞_ｉ｜０＞_ｊ｜ｎ_ｐｈ＋１＞_ｐの重ね合わせに変換される。赤側波帯のπ／２パルスがｊ番目のキュービットに適用される場合、ｊ番目のキュービットと運動モード｜０＞_ｊ｜ｎ_ｐｈ＞_ｐの組み合わせ状態は、｜０＞_ｊ｜ｎ_ｐｈ＞_ｐと｜１＞_ｊ｜ｎ_ｐｈ－１＞_ｐの重ね合わせに変換されるため、組み合わせ状態｜０＞_ｊ｜ｎ＋１＞_ｐは、｜０＞_ｊ｜ｎ_ｐｈ＋１＞_ｐと｜１＞_ｊ｜ｎ_ｐｈ＞ｐの重ね合わせに変換される。

【0028】

したがって、ｉ番目のキュービットにπ／２青側波帯のパルスを適用し、ｊ番目のキュービットにπ／２赤側波帯のパルスを適用すると、２つのキュービットと運動モード｜０＞_ｉ｜０＞_ｊ｜ｎ_ｐｈ＞_ｐの組み合わせ状態を｜０＞_ｉ｜０＞_ｊ｜ｎ_ｐｈ＞_ｐと｜１＞_ｉ｜１＞_ｊ｜ｎ_ｐｈ＞_ｐの重ね合わせに変換することができ、２つのキュービットは、今やもつれ状態にある。当業者にとって明らかであるように、フォノン励起の初期数ｎとは異なる数（すなわち、｜１＞_ｉ｜０＞_ｊ｜ｎ_ｐｈ＋１＞_ｐと｜０＞_ｉ｜１＞_ｊ｜ｎ_ｐｈ－１＞_ｐ）のフォノン励起を有する運動モードともつれる２キュービット状態は、十分に複雑なパルスシーケンスによって除去できるため、ＸＸゲート操作後の２つのキュービットと運動モードの組み合わせ状態は、ｐ番目の運動モードでのフォノン励起の初期数ｎｐｈがＸＸゲート操作の終了時に変化しないので、もつれが解消された（ｄｉｓｅｎｔａｎｇｌｅｄ）と考えてもよい。したがって、ＸＸゲート操作の前後のキュービット状態は、一般に、運動モードを含まずに、以下で説明する。

【0029】

より一般的には、側波帯のパルスを持続時間τ（「ゲート持続時間」と呼ばれる）にわたって適用することによって変換され、振幅Ω^（ｍ）（ｔ）及び離調周波数μ^（ｍ）（ｔ）と、振幅Ω^（ｎ）（ｔ）及び離調周波数μ^（ｎ）（ｔ）とをそれぞれ有するｍ番目とｎ番目のキュービットの組み合わせ状態は、ペアワイズもつれ相互作用χ^{（ｍ，ｎ）}（τ）の観点から次のように記述することができる。

【数2】

ここで、

【数3】

である。

【数4】

は、ｍ番目のイオンと周波数ω_ｐを有するｐ番目の運動モードの間の結合強度を定量化するラムディッケパラメータであり、ψ^（ｍ）（ｔ）はパルスの累積位相

【数5】

であり、

【数6】

は、一般性を失うことなく、簡単にするために以下ゼロ（０）と見なすことができる初期位相であり、Ｐは運動モードの数（鎖１０２内のイオンの数Ｎに等しい）である。

【0030】

（同時ゲートもつれ操作のためのパルスの構築）
上記の２つのキュービットのペアワイズもつれ相互作用を使用して、ＸＸゲート操作（「ペアワイズもつれゲート操作」とも呼ばれる）を実行できる。ＸＸゲート操作（ＸＸゲート）は、単一キュービット操作（Ｒゲート）とともに、所望の計算プロセスを実行するように構成された量子コンピュータを構築するために使用できるユニバーサルゲート｛Ｒ、ＸＸ｝のセットを形成する。ｍ番目とｎ番目のキュービット間でＸＸゲート操作を実行するには、条件χ^{（ｍ，ｎ）}（τ）＝θ^{（ｍ，ｎ）}（０＜θ^{（ｍ，ｎ）}≦π／２）を満たすパルス（すなわち、ペアワイズもつれ相互作用χ^{（ｍ，ｎ）}（τ）が所望の値θ^{（ｍ，ｎ）}を有することで、ペアワイズもつれ相互作用がゼロ以外である条件と呼ばれる）が構築され、ｍ番目とｎ番目のキュービットに適用される。上記のｍ番目とｎ番目のキュービットの結合状態の変換は、θ^{（ｍ，ｎ）}＝π／２のときに最大のもつれを伴うＸＸゲート操作に対応する。ｍ番目とｎ番目のキュービットに適用されるパルスの振幅と離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））と（Ω^（ｎ）（ｔ），μ^（ｎ）（ｔ））とは、ｍ番目とｎ番目のキュービットのゼロ以外の調整可能なもつれを保証して、ｍ番目とｎ番目のキュービットで所望のＸＸゲート操作を実行するように調整できる制御パラメータである。

【0031】

たとえば、ｍ番目とｎ番目のイオンのペア（単に（ｍ，ｎ）と呼ばれる）と，ｍ’番目とｎ’番目のイオンのペア（単に（ｍ’，ｎ’）と呼ばれる）など、イオンの２ペアに対して同時にペアワイズもつれゲート操作を実行するには、振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））、（Ω^（ｎ）（ｔ），μ^（ｎ）（ｔ））、（Ω^（ｍ’）（ｔ），μ^（ｍ’）（ｔ））、及び（Ω^（ｎ’）（ｔ），μ^（ｎ’）（ｔ））を有するパルスは、それぞれｍ番目、ｎ番目、ｍ’番目、及びｎ’番目のイオンに個別に適用される。パルスの振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））、（Ω^（ｎ）（ｔ），μ^（ｎ）（ｔ））、（Ω^（ｍ’）（ｔ），μ^（ｍ’）（ｔ））、及び（Ω^（ｎ’）（ｔ），μ^（ｎ’）（ｔ））が決定され、その結果、パルスの適用の終了時に、もつれるイオンのペア（すなわち、（ｍ，ｎ）、（ｍ’，ｎ’））が互いに結合され、もつれることのないイオンのペア（すなわち、（ｍ，ｍ’）、（ｍ，ｎ’）、（ｎ，ｍ’）、（ｎ，ｎ’））は、互いに分離される。すなわち、もつれるペアの場合、もつれ相互作用がゼロ以外である条件χ^{（ｌ，ｌ’）}（τ）＝θ^{（ｌ，ｌ’）}（（ｌ，ｌ’）＝（ｍ，ｎ）、（ｍ’，ｎ’））が満たされている必要があり、もつれることのないイオンのペアの場合、条件χ^{（ｌ，ｌ’）}（τ）＝０（（ｌ，ｌ’）≠（ｍ，ｎ）、（ｍ’，ｎ’））（つまり、分離された各ペアの全体的なペアワイズもつれ相互作用がゼロ）が満たされる必要がある（分離の条件と呼ばれる）。

【0032】

同時ペアワイズもつれゲート操作は、イオンのペアの数がより多くても実行できる。これらのペアのすべてのイオン（つまり、それぞれが別のイオンともつれるイオン）は、同時ペアワイズもつれゲート操作では「関与イオン」又は「関与キュービット」と呼ばれる。同時ペアワイズもつれゲート操作における関与イオンの数は、以下で、Ｎ_ＥＡＳＥとして記述する。Ｎ_ＥＡＳＥ個の関与イオンに対して同時ペアワイズもつれゲート操作を実行するために、それぞれ振幅Ω^（ｍ）（ｔ）及び離調周波数μ^（ｍ）（ｔ）（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）を有するパルスがｍ番目のイオンに個別に適用される。パルスの振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）は、もつれるイオンのペア（ｌ，ｌ’）については、ペアワイズもつれ相互作用がゼロ以外であるという条件χ^{（ｌ，ｌ’）}（τ）＝θ^{（ｌ，ｌ’）}（０＜θ^{（ｌ，ｌ’）}≦π／２）を満たすように、もつれることのないイオンのペア（ｌ，ｌ’）については、分離する条件χ^{（ｌ，ｌ’）}（τ）＝０を満たすように決定される。

【0033】

制御パラメータであるパルスの振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））は、パルスの照射によって運動モードが励起されたときに、それらの初期位置から変位している鎖１０２内のＮ個のトラップイオンのすべてが、初期位置に戻るという条件を満たさなければならない。重ね合わせ状態│０＞±│１＞のｍ番目のキュービットは、ゲート持続時間τの間にｐ番目の運動モードの励起により変位し、ｐ番目の運動モードの位相空間（位置と運動量）の軌道

【数7】

に従う。軌道

【数8】

は、パルスの振幅Ω^（ｍ）（ｔ）と累積位相

【数9】

によって決定される。ここで、ｇ^（ｍ）（ｔ）は、ｇ^（ｍ）（ｔ）＝Ω^（ｍ）（ｔ）ｃｏｓ（ψ^（ｍ）（ｔ））で定義されるパルス関数である。したがって、Ｎ個のトラップイオンの鎖１０２の場合、条件

【数10】

（つまり、軌道

【数11】

を閉じる必要がある）は、トラップイオンを元の位置と運動量値に戻すための条件（又は位相空間軌道の閉鎖）とも呼ばれ、ペアワイズもつれ相互作用がゼロ以外であるという条件χ^{（ｍ，ｎ）}（τ）＝θ^{（ｍ，ｎ）}（０＜θ^{（ｍ，ｎ）}≦π／２）に加えて、Ｐ個（ｐ＝１，２，…，Ｐ）のすべての運動モードに対して課されなければならない。

【0034】

パルスの制御パラメータである振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）は、結果として得られるパルスが最適な出力になるように調整され、必要なレーザ出力が最小化される（出力最小化の条件と呼ばれる）。必要なレーザ出力はゲート持続時間τに反比例するため、出力最適パルスは、ゲート持続時間τが固定されている場合は最小出力要件で、レーザ出力バジェットが固定されている場合は最短ゲート持続時間τでＸＸゲート操作を実装する。

【0035】

いくつかの実施形態では、振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））は、ゲート持続時間の中間点ｔ＝τ／２に対して時間的に対称又は反対称になるように選択され、すなわち、

【数12】

となる。以下に説明する例では、振幅と離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））は、簡単にするために、対称になるように選択されており（Ω^{（ｍ）（＋）}（ｔ）とμ^{（ｍ）（＋）}（ｔ））、上付き文字（＋）なしでΩ^（ｍ）（ｔ）及びμ^（ｍ）（ｔ）と呼ぶこともできる。対称離調周波数μ^（ｍ）（ｔ）の場合、累積位相ψ^（ｍ）（ｔ）は反対称、つまり

【数13】

となる。トラップイオンを元の位置と運動量値に戻すための条件は、パルス関数ｇ^{（ｍ）（ｔ）}の反対称成分ｇ^{（ｍ）（－）}（ｔ）（単に「パルス関数」と呼び、以下では上付き文字（－）なしで、ｇ^（ｍ）（ｔ）と記す）で、

【数14】

として書き直すことができる。

【0036】

パルスの振幅と離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））で決まるパルス関数ｇ^（ｍ）（ｔ）（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）は，これらの条件を満たすように，パルスの振幅と離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））を変化させ、こうして、パルス関数ｇ^（ｍ）（ｔ）を変化させることで，導出される。ここで、例えば、パルス関数ｇ^（ｍ）（ｔ）は、基底関数ｓｉｎ（２πｋｔ／τ）とパルス関数係数

【数15】

を使用して、ゲート持続時間τにわたってフーリエ正弦基底

【数16】

で展開される。ここで、Ｎ_Ａは、基底関数の数である。トラップイオンを元の位置と運動量値に戻すための条件は、

【数17】

のように書き直すことができる。ここで、Ｍ_ｐｋは

【数18】

として定義される。同様に、トラップイオンを元の位置と運動量値に戻すための条件（たとえば、位相空間軌道の閉鎖）は、行列形式で

【数19】

と書くことができる。ここで、ＭはＭ_ｐｋの２Ｐ×Ｎ_Ａ係数行列であり、

【数20】

は、

【数21】

のＮ_Ａのパルス関数係数ベクトルである。基底関数の数Ｎ_Ａは、運動モードの数Ｐより大きくなるように選択される。したがって、トラップイオンを元の位置と運動量値に戻すための条件を満たし、自明ではないパルス関数ベクトル

【数22】

は、Ｎ_０（＝Ｎ_Ａ－Ｐ）個ある。

【0037】

ペアワイズもつれ相互作用がゼロ以外であるという条件と、分離する条件は、以下：

【数23】

のように書き直すことができる。ここで、

【数24】

は、

【数25】

のように定義されるか、又は同等に、行列形式で、

【数26】

（ｍとｎとがもつれる場合）又は０（それ以外の場合）である。ここで、Ｄ^{（ｍ，ｎ）}は

【数27】

のＮ_Ａ×Ｎ_Ａ係数行列であり、

【数28】

は

【数29】

の転置ベクトルである。トラップイオンを元の位置と運動量値に戻すための条件、及びペアワイズもつれ相互作用がゼロ以外であり、分離するという条件は、原則として、パルス関数係数

【数30】

に関して、二次制約付き二次計画法（ＱＣＱＰ）として知られる形式で記述できることに留意すべきである。一般に、ＱＣＱＰは、非決定論的で困難な多項式時間（ＮＰ）問題であることが知られている（すなわち、少なくともＮＰ問題と同じくらい困難である）。しかしながら、ここで説明する実施形態では、これらの条件はＱＣＱＰの特殊なケースを形成するため、困難なＮＰ問題ではない形式に変換でき、その結果、Ｎ_ＥＡＳＥ×Ｎ_Ａ個の制御パラメータのセット

【数31】

は、関与イオンの数Ｎ_ＥＡＳＥに関して多項式的に増加するオーバーヘッドで決定される。Ｎ_ＥＡＳＥ×Ｎ_Ａ個の制御パラメータのセット

【数32】

を決定するこの方法は、効率的な任意同時もつれ（ＥＡＳＥ）プロトコルと呼ばれ、決定されたパルス関数係数

【数33】

を有するパルスによって実行されるゲート操作は、以下でＥＡＳＥゲートと呼ばれる。

【0038】

出力最小化の条件は、出力関数

【数34】

を最小化することに対応する。これは、ゲート持続時間τにわたって平均化されたパルス関数ｇ^（ｍ）（ｔ）の絶対二乗値である。したがって、出力最適パルスは、トラップイオンを元の位置と運動量値に戻すための条件を満たす、自明でないパルス関数ベクトル

【数35】

の線形結合

【数36】

を計算することによって構築できる。係数Λαは、もつれ相互作用がゼロでないという条件も、出力を最小化する条件も満たすように決定される。

【0039】

したがって、パルスの振幅と離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）は、ペアワイズもつれ相互作用がゼロ以外である条件と出力最小化の条件とを満たすパルス関数係数

【数37】

又は同等のパルス関数係数ベクトル

【数38】

を有するパルス関数ｇ^（ｍ）（ｔ）に基づいて計算できる。トラップイオンを元の位置と運動量値に戻すための条件と、ペアワイズもつれ相互作用がゼロ以外であるという条件は、パルス関数係数ベクトル

【数39】

に関して線形代数の形式であることに留意されたい。したがって、これらの条件を、出力最小化の条件とともに満たすパルス関数係数

【数40】

は、近似又は反復なしで既知の線形代数計算方法によって計算できる。

【0040】

パルス関数係数

【数41】

に関するパルス関数ｇ^（ｍ）（ｔ）の展開は、周波数領域（周波数２πｎ／τ）でのパルスの構築に対応し、したがって、パルス関数ｇ^（ｍ）（ｔ）によって構築されるパルスは、マルチトーンレーザ（すなわち、複数のトーンを有し、各トーンが別個の振幅及び対応する周波数を有するレーザビーム）によって直接実装され得る。すなわち、周波数が２πｋ／τで振幅が

【数42】

のＮ_Ａトーンレーザビーム（ｋ＝１，２，…，Ｎ_Ａ）は、レーザビームの位相が固定されており、ＥＡＳＥゲート操作を直接実装できる。パルス関数は、ゲート持続時間にわたって完全なセット又は不完全なセットを形成する任意の関数を使用して展開できる。しかしながら、パルス関数が不完全なセットで展開された場合、上記の方法で計算されたパルス関数ｇ^（ｍ）（ｔ）が出力最適であるという保証はない。

【0041】

上記の特定の例示的な実施形態は、本開示によるパルス関数の構築方法の単なるいくつかの可能な例であり、パルス関数の構築方法の可能な構成、仕様などを制限するものではないことに留意されたい。例えば、振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））の対称性は、システム１００の構成、仕様等に関連する利便性に基づいて、反対称である（負のパリティを有する）ように、あるいは混合対称を有する（混合パリティを有する）ように選択することができる。しかしながら、振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））に対称性を課すことは、特定の対称性を有する誤差を排除することにつながる可能性がある。

【0042】

（パルスの復調）
出力最適で誤差に強いパルスをｍ番目のキュービットに適用するために、出力最適パルスの振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））を決定されたパルス関数ｇ^（ｍ）＝Ω^（ｍ）（ｔ）ｓｉｎ（ψ^（ｍ）（ｔ））（ｍ＝１，２，…，Ｎ_ＥＡＳＥ）から復調する必要がある（すなわち、振幅Ω（ｔ）及び離調周波数μ（ｔ）振幅及び離調周波数（Ω^（ｍ）（ｔ），μ^（ｍ）（ｔ））が抽出され、パルス関数ｇ^（ｍ）（ｔ）は、単一のレーザビームの一連の時間依存パルスセグメントを有するパルスに変換される）。ここで、

【数43】

は、離調周波数μ^（ｍ）（ｔ）による累積位相である。この復調プロセスが固定離調周波数で実行される場合、すなわち、μ^（ｍ）（ｔ）＝μ_０の場合、結果として生じるパルスは振幅変調（ＡＭ）パルスであり、振幅Ω^（ｍ）（ｔ）は変調される。復調プロセスが固定振幅で、すなわち、Ω^（ｍ）（ｔ）＝Ω_０で実行される場合、結果として生じるパルスは位相変調（ＰＭ）パルスであり、位相ψ^（ｍ）（ｔ）は変調される。位相ψ^（ｍ）（ｔ）が離調周波数μ^（ｍ）（ｔ）を変調することによって実装される場合、結果として生じるパルスは周波数変調（ＦＭ）パルスになる。復調プロセスは、振幅Ω^（ｍ）（ｔ）、位相ψ^（ｍ）（ｔ）（これにより、離調周波数μ^（ｍ）（ｔ））、信号処理の分野で知られている従来の復調方法によって出力最適パルスを構築するための周波数の任意の組み合わせ変調で実行できる。

【0043】

例示的な復調プロセスの最初のステップは、ｔ＝ζ_ｊ（ｊ＝０，１，…，Ｎ_ｚ－１）でパルス関数ｇ^（ｍ）（ｔ）＝Ω^（ｍ）（ｔ）ｓｉｎ（ψ^（ｍ）（ｔ））の零点を見つけることである。ここで、Ｎ_ｚはパルス関数ｇ^（ｍ）（ｔ）の零点の総数である。振幅Ω^（ｍ）（ｔ）は、振幅Ω（ｔ）がゼロにならないように選択できる。したがって、ｓｉｎ（ψ^（ｍ）（ｔ））がゼロの場合（つまり、ｓｉｎ（ψ^（ｍ）（ζ_ｊ））＝０）、パルス関数ｇ^（ｍ）（ｔ）はゼロになる。正弦関数の性質により、ψ^（ｍ）（ζ_ｊ）＝ｊπ（ｊ＝０，１，…，Ｎ_ｚ－１）の場合、パルスのゲート持続時間τの開始時と終了時の零点（すなわち、ｔ＝ζ_０＝０及びｔ＝ζＮ_ｚ－１＝τ）も含め、ｓｉｎ（ψ^（ｍ）（ζ_ｊ））＝０となる。

【0044】

復調プロセスの第二のステップは、パルス関数ｇ^（ｍ）（ｔ）の零点に基づいて離調周波数μ^（ｍ）（ｔ）を計算することである。いくつかの実施形態では、離調周波数μ^（ｍ）（ｔ）は、パルス関数ｇ^（ｍ）（ｔ）の隣接する零点間の定数値として近似される（すなわち、ζ_ｊ－１＜ｔ＜ζ_ｊに対してμ^（ｍ）（ｔ）≒μ_ｊで，ｊ＝１，２，…，Ｎ_ｚ－１）。

【数44】

のように離調周波数μ^（ｍ）（ｔ）により位相ψ^（ｍ）（ｔ）が蓄積されるため、ｔ＝ζ_ｊとｔ＝ζ_ｊ－１の位相差は、

【数45】

となる。したがって、ｔ＝ζ_ｊ－１とｔ＝ζ_ｊの間の離調周波数μ_ｊは、μ_ｊ＝π／（ζ_ｊ－ζ_ｊ－１）して決定される。復調プロセスの第三のステップは、振幅Ω^（ｍ）（ｔ）を計算することである。ｔ＝ζ_ｊでのパルス関数ｇ^（ｍ）（ｔ）＝Ω^（ｍ）（ｔ）ｓｉｎ（ψ^（ｍ）（ｔ））の時間微分は、

【数46】

であり、ここで、ψ^（ｍ）（ζ_ｊ）＝ｊπ及び

【数47】

を使用する。したがって、ｔ＝ζ_ｊでの振幅Ω^（ｍ）（ｔ）は、計算されたパルス関数

【数48】

（つまり、

【数49】

）の時間微分を使用して

【数50】

として計算される。

【0045】

いくつかの実施形態では、計算された離調周波数μ_ｊ（ｊ＝１，２，…，Ｎ_ｚ－１）のセットがスプライン（例えば、１つ以上の多項式及び他の代数式によって区分的に定義された関数）で補間され、離調周波数μ^（ｍ）（ｔ）の補間された値は、振幅Ω^（ｍ）（ζ_ｊ）を計算するためのμ^（ｍ）（ζ_ｊ）に使用される。いくつかの実施形態では、μ^（ｍ）（ζ_ｊ）は、（ｉ）μ_ｊ、（ｉｉ）μ_ｊ＋１、及び（ｉｉｉ）（μ_ｊ＋μ_ｊ＋１）／２であり、振幅Ω^（ｍ）（ζ_ｊ）を計算するためにμ^（ｍ）（ζ_ｊ）として使用される。

【0046】

いくつかの実施形態では、計算された振幅Ω^（ｍ）（ζ_ｊ）のセットも、またスプラインで補間して、時間依存振幅Ω^（ｍ）（ｔ）を計算する。

【0047】

位相変調（ＰＭ）パルスの復調プロセスの場合、計算された位相ψ^（ｍ）（ζ_ｊ）のセットをスプラインで補間して、時間依存位相ψ^（ｍ）（ｔ）を計算することができる。

【0048】

（ＥＡＳＥプロトコル及びＥＡＳＥゲート）
図７は、一実施形態による、制御パラメータ

【数51】

のセット（ＥＡＳＥプロトコル）を決定するための方法７００を示すフローチャートを示す。この例では、Ｎ個のトラップイオンの鎖１０２は量子レジスタであり、その中に、Ｎ個のトラップイオンのそれぞれの２つの超微細状態がキュービットを形成する。

【0049】

本明細書で説明する例では、所望の量子アルゴリズムは、グラフィックスプロセッシングユニット（ＧＰＵ）（図示せず）などのユーザインターフェイスを使用することによって古典的コンピュータ（図示せず）によって選択され、古典的コンピュータ内のソフトウェアプログラムによる複数のキュービットのペアに対して、Ｒゲート操作（単一キュービット操作）と、ＸＸゲート操作（「もつれゲート操作」、「ペアワイズもつれゲート操作」、又は「２キュービット操作」とも呼ばれる）とに分解される。すべてのペアワイズもつれゲート操作の中で、選択されたキュービットのペア（合計でＮ_ＥＡＳＥ個の関与キュービット）に対するペアワイズもつれゲート操作が同時に実行されるように決定され（ＥＡＳＥゲート）、選択されたキュービットのペア間でもつれ相互作用を発生させ、ＥＡＳＥゲート操作を実行するために、Ｎ_ＥＡＳＥ個の関与キュービットに適用されるパルスシーケンスは、古典的コンピュータ内のソフトウェアプログラムによってさらに以下に説明されるように、方法７００（ＥＡＳＥプロトコル）を使用して決定される。ソフトウェアプログラムによって決定されたパルスは、量子レジスタ内のＮ_ＥＡＳＥ個の関与キュービット（Ｎ個のトラップイオンの鎖）に適用され、システムコントローラによって制御される選択されたキュービットのペアに対してＥＡＳＥゲート操作を実行する。

【0050】

ＥＡＳＥプロトコルの開始時に、Ｎ_ＥＡＳＥ個の関与キュービットのそれぞれは、最初に番号ｎ（ｎ＝１，２，…，Ｎ_ＥＡＳＥ）でラベル付けされ、例えば、Ｎ_ＥＡＳＥ個の関与キュービットがＮ個のトラップイオンの鎖１０２内で整列される順序でラベル付けされる。Ｎ_ＥＡＳＥ個の関与キュービットは、最初に、任意の他の順序でラベル付けされてもよい。番号ｎでラベル付けされたキュービットは、以下ではｎ番目のキュービットと呼ばれることもある。

【0051】

ブロック７０２（前処理）では、Ｎ_ＥＡＳＥ個の関与キュービットは、互いに素なキュービットのセットにグループ化される。Ｎ_ＥＡＳＥ個の関与キュービットは再びラベル付けされ、その結果、各互いに素なセット内のキュービットが連続番号でラベル付けされ、各互いに素なセット内で一番目と二番目に小さい番号でラベル付けされたキュービット（互いに素なセットの第一と第二のキュービットと呼ばれる）は、ペアワイズもつれゲート操作が実行される選択されたペアの１つに対応するようになる。例えば、鎖１０２は、１１個のトラップイオン（すなわち、最初に番号１～１１でラベル付けされたキュービット）を有し得、もつれゲートは、キュービット（１，２）、（１，４）、（１，５）、（３，６）、及び（３，８）のペアに対して同時に実行され得る。次に、関与キュービットは、１～６及び８でラベル付けされたキュービットである。第一の互いに素なセットには、キュービット１、２、４、及び５が含まれ、第二の互いに素なセットには、キュービット３、６、及び８が含まれることがある。第一の互いに素なセットのキュービットは１～４で再びラベル付けされ、第二の互い素なセットのキュービットは５～７で再びラベル付けされる。

【0052】

ｎ番目のキュービット（ｎ＝１，２，…，Ｎ_ＥＡＳＥ）に個別に適用されるパルス関数ｇ^（ｎ）（ｔ）を有するパルスを構築するために、最初に、第一のキュービット（１とラベル付け）に適用される第一のパルスのパルス関数ｇ^（１）（ｔ）と、第二のキュービット（２とラベル付け）に適用される第二のパルスのパルス関数ｇ^（２）（ｔ）を決定する。決定した第一のパルスのパルス関数ｇ^（１）（ｔ）と第二のパルスのパルス関数ｇ^（２）（ｔ）に基づいて、第三のキュービット（３とラベル付け）に適用される第三のパルスのパルス関数ｇ^（３）（ｔ）を決定する。次に、このプロセスは、Ｎ_ＥＡＳＥ番目のキュービット（Ｎ_ＥＡＳＥとラベル付け）に適用されるＮ_ＥＡＳＥ番目のパルスのパルス関数ｇ^{（ＮＥＡＳＥ）}（ｔ）を決定するまで続け、こうして、すべてのＮ_ＥＡＳＥ個の関与キュービットに適用されるパルスのパルス関数を決定する。

【0053】

ブロック７０４では、ｎ番目のキュービット（ｎ＝１，２，…，Ｎ_ＥＡＳＥ）に適用されるｎ番目のパルスのパルス関数ｇ^（ｎ）（ｔ）を決定するための最初のステップとして、ｍ（ｍ＝１，２，…，ｎ－１）でラベル付けされたキュービットの中でｎ番目のキュービットと、もつれることのないキュービットは、（「もつれを解消されたキュービット」と呼ばれ、ｓでラベル付けされる）。つまり、ｎ番目のキュービットとｓ番目のキュービット（ｓ≦ｎ－１）とは、ペアワイズもつれゲート操作が実行される選択されたペアのいずれにも対応していない。もつれを解消されたキュービットの数は、以下、Ｎ_ｓとして示す。

【0054】

ブロック７０６では、トラップイオンを元の位置と運動量値（

【数52】

）に戻すための条件を満たすＮ_０個の自明でないパルス関数

【数53】

の中で、ｎ番目のキュービットとｓ番目のもつれを解消されたキュービットの間で分離するための条件（

【数54】

）を満たす１つ以上のパルス関数

【数55】

は、キュービットに適用されるｓ番目のパルスのパルス関数

【数56】

に基づいて、決定される。パルス関数

【数57】

は、以前の反復で定義又は決定される。トラップイオンを元の位置と運動量値に戻すための条件を満たすパルス関数

【数58】

は、（Ｎ_０－（ｎ－１））個あり、これらのパルス関数を線形結合して、

【数59】

にすることができる。

【0055】

ブロック７０４及び７０６のｎ番目のキュービット（ｎ＝１，２，…，Ｎ_ＥＡＳＥ）に適用されるｎ番目のパルスのパルス関数ｇ^（ｎ）（ｔ）を決定するための最初のステップに続いて、すべてのｍ番目のキュービット（ｍ＝１，２，…，ｎ－１）がｎ番目のキュービットから分離されている場合、プロセスはブロック７０８に進む。ｍ番目のキュービットの一部（ｍ＝１，２，…，ｎ－１）（（ｎ－１）番目のキュービットを含む）がｎ番目のキュービットに結合されており、ｍ番目のキュービットの１つ（たとえば、（ｎ－１）番目のキュービット）に適用されるパルスがまだ決定されていない場合、プロセスはブロック７１０に進む。この場合は、すべてのばらばらのセットの第二のキュービットに対して発生する。ｍ番目のキュービット（ｍ＝１，２，…，ｎ－１）の一部がｎ番目のキュービットに結合され、ｍ番目のキュービットの全部に適用されるパルスが決定された場合、プロセスはブロック７１２に進む。

【0056】

ブロック７０８では、ｎ番目のキュービットのパルス関数

【数60】

を保存する。この場合、すべてのｍ番目のキュービット（ｍ＝１，２，…，ｎ－１）がｎ番目のキュービットから分離される。プロセスは、ブロック７０４に戻り、（ｎ＋１）番目のキュービットに適用される（ｎ＋１）番目のパルスのパルス関数ｇ（ｎ＋１）を決定する。

【0057】

ブロック７１０では、（ｎ－１）番目のキュービットに適用される（ｎ－１）番目のパルスのパルス関数ｇ^{（ｎ－１）}（ｔ）は決定されていない。したがって、ｎ番目の反復のブロック７０６で決定されたパルス関数

【数61】

の線形結合

【数62】

と、同様に、（ｎ－１）番目のキュービットの前の反復のブロック７０６で決定されたパルス関数

【数63】

の線形結合

【数64】

とを計算し（すなわち、係数

【数65】

及び

【数66】

を決定する）、その結果、ｎ番目と（ｎ－１）番目のキュービット間のもつれ相互作用がゼロ以外であるという条件（

【数67】

）が満たされ、結果として生じるｎ番目と（ｎ－１）番目のパルスのノルム（ｎｏｒｍｓ）は、ゲート持続時間τにわたって平均化され、

【数68】

及び

【数69】

は最小化される。上述のｎ番目と（ｎ－１）番目のパルスは、出力が最適になるように構成される（すなわち、必要なレーザ出力が最小化される）。プロセスは、ブロック７０４に戻り、（ｎ＋１）番目のキュービットに適用される（ｎ＋１）番目のパルスのパルス関数ｇ^{（ｎ＋１）}を決定する。

【0058】

ブロック７１２では、ブロック７０６で決定されたパルス関数

【数70】

の線形結合

【数71】

を計算して（すなわち、係数

【数72】

を決定する）、その結果、ｎ番目とｍ番目（ｍ＝１，２，…．ｎ－１）のキュービット間のペアワイズもつれ相互作用がゼロ以外であるという条件（ｎ番目とｍ番目のキュービットがもつれている場合、

【数73】

、その他の場合、０）が満たされ、結果として生じる、ゲート持続時間τにわたって平均化されたｎ番目のパルスのノルム

【数74】

は、最小化される。上述の（ｎ－１）番目のパルスは、ノルムが最適になるように構成される（すなわち、必要なレーザ出力が最小化される）。プロセスは、ブロック７０４に戻り、（ｎ＋１）番目のキュービットに適用される（ｎ＋１）番目のパルスのパルス関数ｇ^{（ｎ＋１）}を決定する。

【0059】

すべてのＮ_ＥＡＳＥ個の関与キュービットのパルスが構築されると、プロセスは終了する。

【0060】

上記のように構築されたパルスを関与キュービットに適用することは、選択された量子アルゴリズムが分解される一連のユニバーサルゲート｛Ｒ、ＸＸ｝操作の中で、関与キュービットのペアに、ペアワイズもつれゲート操作（ＸＸゲート操作）を実装する。一連のユニバーサルゲート｛Ｒ、ＸＸ｝操作のすべてのＸＸゲート操作（ＸＸゲート）は、選択された量子アルゴリズムを実行するために、単一キュービット操作（Ｒゲート）とともに上記の方法７００によって実装され、これは、古典的コンピュータによって定義され、実装される。選択された量子アルゴリズムの実行の最後に、量子レジスタ（トラップイオンの鎖１０２）内のキュービット状態（トラップイオン）の集団が、イメージング対物レンズ１０４で得られた測定によって決定（読み出し）され、ＰＭＴ１０６にマッピングされる。こうして、選択された量子アルゴリズム内の量子コンピューティングの結果を決定し、古典的コンピュータ（例えば、デジタルコンピュータ）への入力として提供できるようにする。次に、量子コンピューティングの結果を古典的コンピュータで処理することで、所望のアクティビティを実行したり、古典的コンピュータだけでは通常は解明できず、妥当な時間内に解明できない問題の解決策を取得したりすることができる。今日、従来のコンピュータ（すなわち、古典的コンピュータ）によって手に負えない、及び解明できないことが知られているが、実行された量子コンピューティングから得られた結果を使用することによって解決できる問題には、複雑な分子及び材料の内部化学構造のシミュレーションと、大きな整数の因数分解とが含まれるが、これらに限定されない。

【0061】

本明細書に記載のＥＡＳＥプロトコルは、効率的な方法で、誤差又は近似なしに、キュービットの複数ペアに対して、ペアワイズ同時もつれゲート操作を実行するためのパルスを決定することができる。すなわち、時間セグメントの数など、パルスを決定するためのオーバーヘッドは、関与キュービットの数に対して直線的に増加するだけである。これは、関与キュービットの数が増えるにつれて指数関数的に増加するオーバーヘッドを必要とする、以前に提案された非線形で近似的な方法とは対照的である。さらに、ＥＡＳＥプロトコルによって構築されたパルスは、単一のＸＸゲート操作に使用する場合に最適であり、振幅変調を使用して所望のもつれを実現する場合に、ゲート操作を実行するために必要なレーザ出力が最小になる。したがって、量子コンピュータは、パルスを実装するためのレーザ出力やパルスを決定するための計算リソースなどの利用可能なリソースを使用して、より高速な実行速度で、より大きなサイズにスケールアップできる。

【0062】

上記は特定の実施形態を対象としているが、他のさらなる実施形態は、その基本的な範囲から逸脱することなく考案することができ、その範囲は、以下の特許請求の範囲によって決定される。

【図1】

【図2】

【図3A】

【図3B】

【図3C】

【図4】

【図5】

【図6】

【図7】